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Abstract 



In this thesis wc focus on studying the physics of cosniological recombina- 
tion and how the details of recombination affect the Cosmic Microwave Back- 
ground (CMB) anisotropies. We present a detailed calculation of the spectral 
line distortions on the CMB spectrum arising from the Lya and two- photon 
transitions in the recombination of hydrogen (H), as well as the correspond- 
ing lines from helium (He). The peak of these distortions mainly comes from 
the Ly a transition and occurs at about 170 /im, which is the Wien part of the 
CMB. The detection of this distortion would provide the most direct supporting 
evidence that the Universe was indeed once a plasma. 

The major theoretical limitation for extracting cosmological parameters from 
the CMB sky lies in the precision with which we can calculate the cosmologi- 
cal recombination process. Uncertainty in the details of hydrogen and helium 
recombination could effectively increase the errors or bias the values of the cos- 
mological parameters derived from microwave anisotropy experiments. With 
this motivation, we perform a multi-level calculation of the recombination of 
H and He with the addition of the spin-forbidden transition for neutral helium 
(Hei), plus the higher order two-photon transitions for H and among singlet 
states of He I. Here, we relax the thermal equilibrium assumption among the 
higher excited states to investigate the effect of these extra forbidden transitions 
on the ionization fraction Xc and the CMB angular power spectrum C'e- We find 
that the inclusion of the spin-forbidden transition results in more than a percent 
change in , while the higher order non- resonance two-photon transitions give 
much smaller effects compared with previous studios. 

Lastly we modify the cosmological recombination code recfast by introduc- 
ing one more parameter to reprodiice recent mimca'ical results for the speed-up 
of helium rccoiribination. Together with the existing hydrogen 'fudge factor', we 
vary these two parameters to account for the remaining dominant uncertainties 
in cosmological recombination. By using a Markov Chain Monte Carlo method 
with Planck forecast data, we find that we need to determine the parameters to 
better than 10% for He I and 1% for H, in order to obtain negligible effects on 
the cosmological parameters. 
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The detection of the 2.725 K Cosmic Microwave Background (CMB) is one of 
the strongest pieces of supporting evidence for the Big Bang model, which is 
the widely accepted theory for the history of the Universe. Together with other 
observations, we know that the Universe is expanding implying that it was much 
denser and hotter in the past and used to be a plasma of ions and electrons. 
The CMB, which is the remnant of the early radiation, was last scattered when 
the atoms became neutral. This period is called cosmological recombination, 
and it happened when the Universe was a few hunderd thousand years old. 
In the decades following its discovery, the CMB was found to be remarkably 
homogeneous and isotropic, but its tiny temperature fluctuations give us the 
most distant image we have of the Universe. This carries important information 
about the geometry, the expansion rate and contents of the Universe, as well as 
clues about the origin of all the structure it contains (see, for example, |2^). 
Exploiting this information requires an extremely precise undertanding of the 
process of cosmological recombination, which is the main topic of this thesis. In 
order to explain why this is the case we should first review the physics of the 
standard cosmological model. 



1.1 A brief history of the Universe 

In the late 1920s, Hubble |TJ] discovered that the Universe is expanding. He 
found that atomic lines in the spectrum of nearly all distant galaxies are red- 
shifted (or shifted to longer wavelengths) compared with the laboratory values. 
This means that the galaxies are moving away from us due to the expansion of 
the Universe. The redshift z is defined as 

where Aobs and Aomit are the observed and emitted wavelengths, respectively. 
Here Bify is a time-dependent scale factor, which gives infinitestimal distances 
in space when multiplied by the comoving distance dr. This idea of a uniform 
scale factor for the expansion is consistent with Hubble finding that the velocity 
of galaxies v increases linearly with distance r, which is the famous Hubble's 
law: 

v = Hr. (1.2) 
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Here H is the Hubble constant and represents the rate of expansion so that 



and today we have HQ=H{to). Although the actual value of the constant 
determined by Hubble is far from our current estimates, the Hubble diagram 
nevertheless proves that the Universe is expanding, and the same principle is 
used for today's measurements: measure the redshifts and estimate the distances 
of distant objects to determine H. Rcdshift, can be easily estimated from the 
shifting of the spectral lines, but it is hard to determine the distances without 
any information of on the intrinsic brightness or the intrinsic size of an object, 
so that precision measurements of Hubble's constant have been elusive. 

The current value of the Hubble constant Hq (the subscript '0' represents 
the present value, that is at z = 0) was determined by the Hubble Key Project 
using 'standard candles', which basically have the same intrinsic brightness or 
have a correlation between some observables and the intrinsic brightness. For 
example, Cepheid variables and Type la supernovae are commonly used stan- 
dard candles. The measured value of Hq is equal to 72±8kms^^Mpc^^ [Q. We 
usually define a dimensionless constant for Hq, which is 



and therefore, /i = 0.72 ± 0.08. Assuming R(t) is constant, the age of the Uni- 
verse is then equal to I/Hq, which is about 13.7 Gyr. 

The Universe appears to be homogeneous and isotropic on large scales (dis- 
tances greater than about 300 Mpc) from observations of the distribution of 
galaxies ||2^ . This is the Cosmological Principle; based on that we can build 
a simple model of the expanding Universe within General Relativity. Here we 
temporarily ignore the density fluctuation on small scales, which are of small 
amplitude in the early Universe but important later for the formation of galax- 
ies and clusters (the structure formation). On large scales, the Universe can be 
described by the Friedmann-Robertson- Walker (FRW) metric and the geometry 
of the Universe depends on the total density (see, for example, H^, |2^). Given 
the expansion rate H, there is a critical density pcr that determines whether the 
Universe has flat geometry. This critcal density is 




(1.3) 



h = 



Ho 



(1.4) 



lOOkms-iMpc 



-1 



3i/2 



(1.5) 



Per = 



BttG ' 



and we usually define a density parameter 



= Pi/ p, 



(1.6) 



where i represents different components (e.g. matter, radiation and dark energy) 
in the Universe. The Universe is spatially closed if the total density of the 
Universe is larger than pcr, and spatially open if its density is lower than pcr- 
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Figure 1.1: A schematic picture of a brief history of the Universe. Boxes indicate 
the periods when radiation, matter or dark energy are dominant. 



From the combined results of recent observations of the CMB, acoustic sig- 
natures in galaxy clustering and Type la supernovae, the Universe is found to be 
very close to flat with the total density fitot = 1.0052± 0.0064 Q. At the present 
time the Universe consists of about 4% baryons (ilh), 20% cold dark matter (fie), 
76% dark energy (Ha) and a tiny portion of photons (riR = 4.17x lO^^h^^). 
Here 'baryon' means ordinary matter, for example atoms, nuclei and elec- 
trons. 'Dark matter' is some gravitationally interacting (weakly interacting 
with baryons) and non-luminous substance. The dark matter is considered to 
have velocity dispersion which is negligible for structure formation, meaning 
that it decoupled when it was non-relativistic (cold); its fluctuations are the 
seeds of structure formation. The concept of dark matter was first proposed 
by Zwicky in 1933 through observations at the rotational curves of stars 
in galaxies. This dark matter was introduced in order to explain the increas- 
ing rotational velocity of material with increasing distance from the centres of 
galaxies. As we will see later, since Big Bang Nucleosynthesis (as well as the 
CMB) gives a very low limit on the baryon density, some non-baryonic matter 
must exist in the Universe. In general, the density of matter is proportional 
to (1 -|- z)^, while that of the radiation is proportional to (1 -I- z)"^. The 
dark energy provides the negative pressure responsible for the recent acceler- 
ated expansion of the Universe and its density is constant over redshift (this is 
Einstein's cosmological constant, A) or very nearly so. 

Due to different scalings of the density of each species with redshift, the 
components dominate the Universe at different times. Figure |I.1| shows a brief 
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history of the Universe and indicates the important epochs using both time and 
redshift as coordinates. In cosmology, redshift z is usuaUy used instead of time, 
since it is (in principle at least) directly observable and so independent of the 
cosmological model. In the Big Bang picture, the very early times are still quite 
uncertain, but the physics of the thermal history of the Big Bang Nucleosyn- 
thesis (BBN) and recombination are well understood and firmly established. 

The earliest times were radiation dominated. The Universe was very hot (the 
background radiation temperature T ~ To{l + z), where Tq = 2.725 K) and dense. 
Due to the strong and highly energetic photon background, there were no bound 
nuclei until BBN occured at about 3 minutes after the Big Bang (z ~ 10^ — 10^). 
During BBN, the temperature decreased to about 100 keV/fce, which is lower 
than the typical binding energy of the nuclei. Therefore, nuclei of deuterium (D), 
helium (^He, ''He) and lithium (^Li) were able to form without being destroyed 
by the photons. Given the baryon density fib, the theoretical calculation of 
standard BBN can predict the abundance of different species of nuclei with very 
small uncertainties due to nuclear and weak-interaction rates (see Figure 1 in 

or Figure 5 in ]35|] ). In particular, the abundance of D is very sensitive to 
r^b . By measuring the primodial abundance of D through the absorption lines 
in the hydrogen clouds at redshift z ~ 3 — 4, we can put tight constraints on 
flh using the theoretical BBN prediction (see ||] and references therein). BBN 
gives a limit that the baryons can contribute at most 5% of the critical density, 
and therefore the rest of the matter must be non-baryonic. 

At about 3 X 10^ years (z ~ 1100) after the Big Bang, the radiation temper- 
ature dropped to around leV/fce, which is lower than the ionization energy of 
typical atoms. This period is called cosmological recombination. During this 
time, the ions and electrons were able to bind together without being ionized 
by the background photons. After the Universe became neutral, the photons 
were no longer scattered by the electrons and could basically travel freely to the 
present, being redshifted in the expanding Universe. These are the CMB pho- 
tons that we detect today. The CMB has been found to be remarkably smooth, 
the amplitude of the temperature deviations AT/T is only about 10~^, which is 
a strong contrast to the non-linear structure formed by the galaxies and clusters 
we observe today. Therefore this fluctuation amplitude of temperatures in the 
CMB gives us an idea about the strength of the matter density fluctuations 
at the time of recombination, which evolved into the large scale structures we 
observe now. 

After recombination, the Universe remained dark and neutral (20 < z < 900) 
until the first stars formed. There has not been any detection of informtion from 
this 'dark age' and we are still not sure how and when exactly the first stars 
formed. Up until now, the most distant quasar that has been observed is at 
about z = 6.5 ||l5|, From the hydrogen absorption line spectra from such 
high-z quasars M we know that the Universe was fully ionized by ultraviolet 
radiation from hot stars at z < 6. Moreover the CMB provides a constraint 
on the optical depth Trcion during this reionization epoch through the Thomson 



Chapter 1 . Introduction 



5 



scattering effect on the photons. The integrated optical depth is 

caTnc{z)^dz , (1.7) 

where ctt is the Thomson scattering cross-section, is the number density of 
free electrons and 2;rcion is the redshift at which the Universe became ionized. 
From the latest CMB measurement and assuming that the Universe became 
fully ionized instantaneously, the current estimate is Zroion — H Stars and 
galaxies are created basically due to the gravitation collapse of dense regions, 
but the process is non-linear and also involves the pressure of the gas. Therefore, 
although the current matter inhomogeneites in the Universe and the temperture 
fluctuations of the CMB originated from the same source, they appear very 
different today. 

In inflationary models, the primodial perturbations are generated by quan- 
tum fluctuations (see Q for a general review). For the simplest model, by 
assuming the matter is adiabatic and its fluctuations are Gaussian, the initial 
conditions for density perturbations can be described by only two parameters: 
the scalar amplitude As and the spectrum index rig (the slope of the power 
spectrum; the subscript 's' distinguishes these scalar perturbations from pos- 
sible tensor, or gravity wave, contributions). The variance of the comoving 
curvature perturbations is usually defined as 

where = A^(fco), k is the wavenumber and fcp — 0.05 Mpc^^. 

Since the CMB photons come from the time before stars formed, the anisotro- 
pics in the CMB provide us with information about density perturbations at the 
recombination time and in combination with measurements made today, they 
are a powerful tool for constraining the parameters of the cosmological model. 
From the above discussion, and assuming a flat Universe, the standard cosmolog- 
ical model (the A Cold Dark Matter model, ACDM) consists of six parameters: 
rib, ^^m, Trcion, and Tig. There could of course be more parameters in the 
cosmological model (see for a review), for example, including the tensor 
mode of the primodial perturbations or allowing the Universe to deviate from 
flatness (f^tot 7^ !)■ 

Since the CMB photons were mostly last scattered during the epoch of cos- 
mological recombination, we need to understand in detail how the photons de- 
coupled from the matter during that period in order to obtain the correct CMB 
anisotropy power spectrum for constraining the cosmological parameters using 
the observations. In this thesis, we focus on the physics of recombination and 
how the details of the recombination process affects the CMB. We now therefore 
present an introduction to the physics of cosmological recombination (the last 
scattering surface of the CMB photons), and also the basic principles of the 
formation of the CMB anisotropics. 
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1.2 Cosmological recombination 

Recombination in an expanding Universe is not an instantaneous process. It 
is basically controlled by the recombination time and by the Hubble expansion 
time. If the recombination time is much shorter than the expansion time, then 
the electrons and ions follow an equilibrium distribution. For the ionization of 
a plasma, the equilibrium situation is described by the Saha equation. Taking 
hydrogen as an example (see Equation (13) in |29j and references therein), 

l2 \ 3/2 



Here is the number density of electrons in the ith energy level of the H atom, 
rip is the number density of free protons, TOc is the mass of the electron, fee is 
the Boltzmann constant, hp is Planck's constant, gi is the degeneracy of the 
energy level i and Ei is the ionization energy of level i. Due to the higher 
ionization energy, helium recombined at higher redshifts, first by forming He"*" 
(Hen) and then neutral He (Hei). Hydrogen started to recombine shortly after. 



Figure 1.2 shows the full ionization history of recombination by plotting the 
ionization fraction (xc^nc/nn, where njj is the number density of H nuclei) 
versus z. Based on standard BEN, about 8% (by number) of the atomic nuclei 
are helium. And since the ionization fraction is normalized to the total 
number density of hydrogen, is equal to about 1.16 when the Universe is 
fully ionized. 

Peebles (1968) [|2| and Zeldovich (1968) Q first calculated the Hi recom- 
bination evolution in detail and found that the recombination process is much 
slower than Saha equilibrium (for example, see Figure 6.8 in|Q). The Saha 
equation is good for describing the initial departure from full ionization, but the 
equilibrium situation breaks down shortly after recombination starts. When the 
temperature of the Universe reached about CSeV/fce at zcil700, there were 
not enough photons in the Wien tail to keep ionizing the H atoms. Due to the 
high photon to baryon ratio n^/rib — 10^, direct recombinations to the ground 
state were highly prohibited. The 'spectral distortion' photons emitted from di- 
rect recombination are highly energetic and easily re-ionize the nearby neutral 
atoms. This is very similar to the 'Case B' recombination familiar in other areas 
of astrophysics (see e.g. ^^), in which the electrons mostly cascade down to the 
ground state through the first excited state n — 2. However, in cosmological Hi 
recombination, the resonant 2p-ls Lya transition is also strongly suppressed, 
because the line is optically thick. These line photons can only escape reab- 
sorption through redshifting out of the line and the probability for this is very 
low. The other way for the electrons to move from the first excited state to the 
ground state is through the 2s-ls two-photon forbidden transition. Almost half 
of the electrons cascade down from the n = 2 state through this process (see 
Chapter 2 & 3 for details). Overall, the net recombination rate to ground state 
from n = 2 state is lower than the recombination rate into the n — 2 state, and 
this causes a 'bottleneck', which is responsible for making the net recombination 
rate much smaller than the one given by Saha equilibrium. 
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Figure 1.2: The ionization history for cosmological recombination generated by 
the current version of recfast. The dashed line shows the visibility function 
g{z) as a function of redshift (multipied by 100 for better illustration). The 
cosmological ACDM model used here has: r2b = 0.04; 17^ = 0.24; f2A = 0.76; 
/i = 0.70; yp = 0.25; and Tq = 2.725 K. 
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In the next chapter, we will discuss details of the radiative processes during 
recombination and also recent development in performing the numerical calcu- 
lations. However, all the updates are based on the basic picture of the standard 
recombination given here. We have already discussed how H I recombination is 
not an equilibrium process. The situation is similar for helium recombination. 
He I recombination is also slower than Saha equilibrium due to the 'bottleneck' 
at the first excited state, but Hell deviates from the Saha value at only the 0.2% 
level due to the relatively fast two-photon rate to the ground state ^ . 
The ionization fraction Xc affects the CMB anisotropics (see Equation 



( |1.13D for the definition of Ci) through the shape of the last scattering surface 



which is given by the visibility function g{z), 

, , ^dr , ^ 

9{z)^e--—, (1.10) 

where t is the Thomson optical depth during recombination (excluding the 
effects of reionization if we are only considering primary anisotropics). Here r is 
defined the same as in Equation ( |1.7| ) , but with different integration limits (say, 
from z = c» to 100). One can consider g{z) as the probability that a photon 
last scattered at redshift z. In Figure |1.2| , the function g{z) is plotted on top of 
the ionization history of cosmological recombination. Since t changes rapidly 
with z, g{z) is sharply peaked, and its width gives us the thickness of the last 
scattering surface (which means that the CMB photons we see last scattered in 
the specific range of redshift 600 < z < 1500). It is usual to define the location of 
the peak of g{z) as the redshift of the recombination epoch, when the radiation 
effectively decoupled from the matter z^cc- This is approximately equal to 1100 
in the current cosmological ACDM model. From the profile of g{z), we can see 
that Hi recombination affects the Ce much more than He. The later stages of 
He I recombination can also change the high-z tail of g{z) (see Chapter 6 for 
more details) , but He ii recombination occurs too early to bring any significant 
effects on d. 



1.3 Cosmic microwave background 

From many measurements, particularly those of the Far-InfraRed Absolute 
Spectrophotometer (FIRAS) on board with the Cosmic Background Explorer 
{COBE) § § [19|, the CMB was found to be very close to a pure blackbody 
spectrum, which is described by the Planck function J: 

J- nil) 

Figure |l.3| shows the data points from FIRAS (|, D , with error bars multiplied by 
100 and compared with the theoretical blackbody spectrum with Tr = 2.725 K. 
We can see that the data points match the blackbody shape incredibly well 
within the frequency v range from 2 to 20cm~^ (i.e. 60 to 600 GHz). The 
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deviation is less than 5 x 10^^ at the peak of the CMB spectrum Q. The back- 
ground photons originate from an epoch much earher than that of recombi- 
nation, coming from the electron-positron annihilations before BBN and from 
when the energy of the photons was so high that bremsstrahlung and double 
Compton scattering could create and destroy photons so that they were rapidly 
thermalized into a blackbody spectrum p6[ . Hence spectral distortion constrain 
any energy injection later than that epoch. The FIRAS data put strong limits 
on the chemical potential < 9 x 10~^ and the Compton-scattering distortion 
parameter |?/| < 1.5 x 10~^ ||, These strong constraints eliminated many 

earlier competing cosmological models and provide strong evidence that the ra- 
diation temperature Tr scales accurately as (1 + ^) (see, for example, [|2^, ^ 
for more details). The small value of y shows that the hydrogen remained neu- 
tral for quite a long time, otherwise distortions of the blackbody spectrum due 



to Compton scattering by the hot electrons would be observed (see |37 and 
references therein). 

The other main feature of the CMB is the dipole variation of the temperature 
across the sky, with an amplitude equal to 3.358 mK (see ||2^ for a review). This 
anisotropy is determined by the Doppler shift from the solar system's motion 
relative to the 'rest frame' of the radiation, which is supported by measurements 
of the radial velocities of relatively local galaxies. When we talk about the 
temperature anisotropics of the CMB, this contribution from our relative motion 
is usually removed. 

The first detection of the CMB temperature anisotropics was made by the 
COBE Differential Microwave Radiometer (DMR; ) . The variations in tem- 
perature, AT/T, were found to be of the order of 10^'^. We usually decompose 
maps of the CMB temperature fluctuations using the spherical harmonic expan- 
sion: 

^ - Zl^l^ . Y: a,,MO, </>) . (1.12) 

i,m 

If the fluctuations are Gaussian and the sky is statistically isotropic (indepen- 
dent of m), then the temperature field is fully charaterized by the amplitudes 

{^*l,m^i' ,m') = Su'd.mm'Cl . (1-13) 

We usually plot £{£ + 1)Ci/2tt, since this is the contribution to the variance of 



the power spectrum per logarithmic interval in £ (see, for example, [|13, 27, 4l|). 
The radiation temperature itself corresponds to the monopole £ ^ 0, while the 
dipole variation corresponds to ^ = 1. 

Temperature fluctuations in the CMB are essentially a projection of the mat- 
ter density perturbations at the recombination time. There are many reviews 
covering details of the formation of the CMB anisotropics (see |l^, |27| and ref- 
erences therein) and we just briefly recount the basic mechanism here. Photons 
from high density regions were redshifted when they climbed out of the poten- 
tial wells (the Sachs- Wolfe effect). And the adiabaticity between matter and 
photons also gives a higher temperature in higher density regions. The other 
primary source is the oscillating density and velocity of the photon fluid itself. 
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Figure 1.3: Intensity of cosmic microwave background radiation as a function of 
frequency. The crosses are the data points from FIRAS ||^, ^ and the sohd Hne 
is the expected intensity from a pure blackbody spectrum with Tr = 2.725 K. 
Note that the plotted one-sigma error bars have been magnified by 100. Other 
experiments extend the frequency range, but typicaUy with much larger errors, 
and add nothing substantially new to the constraints on the spectral shape. 
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Before the epoch of recombination, the baryons and the radiation are tightly 
coupled as a single photon-baryon fluid, through Thomson and Compton scat- 
terings. The structure seen in the anisotropy power spectrum is mainly due to 
the acoustic oscillations in this photon-baryon fluid, driven by the evolving per- 
turbations in the gravitational potential. One can think of these oscillations as 
standing waves in a harmonic series, with the fundamental mode being the scale 
which has reached maximal compression at the time of last scattering. After 
recombination, when the Universe became neutral, the photons decoupled from 
the atoms and could propagate freely to us (although there are some secondary 
anisotropics formed when the photons travel along the line of sight). Therefore, 
the correct interpretation of the relationship between the underlying matter fluc- 
tuation spectrum and the photon distribution depends strongly on the angular 
diameter distance between us and the last scattering surface. This distance 
depends on the expansion and curvature of the Universe or equivalently, the 
energy content of the Universe. Therefore, the CMB temperature anisotropics 
can provide precise constraints on the cosmological expansion model, as well as 
the scale dependence of the primodial fluctuations. 

In addition, the Thomson scattering between electrons and photons also 
leaves a characteristic signature in the polarization of the CMB photons. The 
quadrupole temperature anisotropy in the photon field generates a net linear 
polarization pattern through Thomson scattering. It has became conventional 
to decompose the polarization pattern into two modes: a part that comes from a 
divergence ('B-mode'); and another part from a curl ('B-mode'). Scalar pertur- 
bations (i.e. spatial variations in density) coming from the inflation epoch only 
give an E-mode signal, while tensor perturbations (i.e. gravity waves) produce 
both E and B-modes. Much current activity in CMB experiments is focussed 
on trying to measure these B-modes, in order to probe the physics of inflation. 
In fact, there are 6 possible cross power spectra from the full temperature and 
polarization anisotropy data set. Cross-correlation between the B-mode and 
either the T or E-mode is zero due to having opposite parity. This leaves us 
with 4 possible observables: Cf^, Cf ^, Cj^ and Cf ^. 

Figure L4 shows the anisotropies Cf^ and Cj^ with £>2 from recent result 
based on the Wilkinson Microwave Anisotropy Probe ( WMAP; jo) ) 5-year data. 
The points show the WMAP data, while the sohd fine is the best-fit ACDM 
model. We can see that the first two acoustic peaks of the temperature spectrum 
are well measured and there is clearly a rise for the third peak. Together with 
other ground based experiments (see and references therein), perhaps the 
first five acoustic peaks have now been localized. 

The amplitude of the polarization signal is about 2 orders of magnitude 
smaller than the temperature one and so it is much harder to detect. The 
DASI jl^ experiment first demonstrated the existence of CMB polarization in 
2002 and the WMAP experiment has measured the TE power spectrum to 



high precision . Figure 1.4 shows the recent measurements of Cf^ from the 
WMAP 5 year results. 
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Figure 1.4: The temperature auto-correlation (TT) and temperature- 
polarization cross-correlation (TE) power spectra with 2 < £ < 1000. The 
points are from the 5 year WMAP data and the error bars are the noise er- 
rors only. The solid line is the best-fit 6 parameter ACDM model, fit to the 
WMAP data only j|] . The grey shaded area shows the 1 a error band due to 
cosmic variance (i.e. the fact that our realization of the CMB sky can vary 
from the underlying expectation value). This figure is taken from Hinshaw et 
al. (2008) §. 
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1.4 Why are we interested in recombination? 

1.4.1 Distortion photons from recombination 

From the previous section, we know that the photons in the radiation back- 
ground were thermahzed to a nearly perfect blackbody spectrum by bremsstrah- 
lung and double Compton scattering processes before recombination. As well 
as the photons from this blackbody background, there were some extra distor- 
tion photons produced during the epoch of cosmological recombination. When 
an electron combined with an ionized atom and cascaded down to the ground 
state, there was at least one distortion photon emitted for each recombination. 
These recombination photons give a distinct series of spectral line distortions 
on the nearly perfect blackbody CMB spectrum. The main contribution to the 
distortion comes from the Hi Lya transition at about z ~ 1500, and this line 
will be observed in the Wien tail (~ 100 /xm) of the CMB spectrum today (see 



Figure 3.3 in Chapter 3). Since these distortion photons are produced directly 
from each recombination of the atoms, the overall shape and amplitude of the 
line are very sensitive to the details of the recombination process. Therefore 
the detection of this spectral distortion would provide direct contraints on the 
physics of recombination and also provide incontrovertible evidence that the 
Universe was once a hot, dense plasma which recombined. 

FIRAS showed that the CMB spectrum around the peak is well-modelled by 
a 2.725 K Planck spectrum. It was found that there is also a Cosmic Infrared 
Background (CIB; see |§, |[ |6)), which peaks at about 150 /zm, right above the 
recombination distortion on the CMB spectrum (see Figure |3.3|) . This back- 
ground is mainly due to luminous infrared galaxies at fairly recent epochs and 
it makes the detection of the recombination distortion even more challenging. 
The first calculations of the line distortion on the CMB tail were presented by 
Peebles (1968) H and by Zeldovich et al. (1968) @. However, they provided 
no details about the line shape, and since then there have been no explicit cal- 
culations showing different contributions to the line shape. Today we have a 
better understanding of the cosmological model as well as improved detection 
techniques, and so it is time to calculate these spectral distortion lines to much 
higher accuracy, in order to investigate whether they could be detected and 
whether such a detection would be cosmologically interesting. A detailed study 
of this line distortion on the CMB spectrum coming from the recombination 
time will be presented in Chapter 3. 



1.4.2 Precision cosmology 

The CMB anisotropics have been well studied theoretically, and the calculations 
are robust, because they can be based on linear perturbation theory (see p3[ 
and references therein), given that the primordial fluctuations are of small am- 
plitude. CMBFAST 1^ is one of the most widely used numerical Boltzmann 
codes for calculating the C^. It has been tested over a large set of cosmological 
models and is consistent with other codes with an accuracy at better than the 1% 
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level [|2| . We have already entered the era of precision cosmology [T^, ^ ^ . 
With the release of the WMAP 5 year data, we can constrain the cosmological 
parameters extermely well from the shape of the anisotropy power spectrum Q . 
The next generation of CMB satellites, Planck |^ , which will be launched in 
early 2009, has been designed to sensitively measure the Cg of the TT- and 
T&modes up to £ = 2500 and the EE^mode for £ < 2000. In order to extract 
the correct cosmological parameters from the experimental data, theoretical cal- 
culations with consequently higher accuracy are required. It now seems clear 
that we need to obtain the theoretical CiS to better than the 1% level. And the 
main theoreical uncertainty comes from details of the ionization history during 
recombinatio n |]3^ . 

RECFAST [^8[ is the most common numerical code for calculating the evolu- 
tion of the ionization fraction Xc during recombination; it is embedded into all 
of the widely-used Boltzmann codes. It is written to be a short and quick pro- 
gram for reproducing the results from a multi- level atom calculation [p9[ , which 
follows the evolution of the number density of electrons at each of more than 100 
atomic levels for each species of atom. The accuracy of the Xc obtained from 
RECFAST is at the percent level, which is sufficient for WMAP, but may not 
be good enough for Planck. This fact has recently motivated many researchers 
to investigate several detailed physical processes during recombination which 
may cause roughly percent level changes on Xe- Although the basic physical 
picture for standard cosmological recombination is quite well established, the 
non-equilibrium details of recombination are unexpectedly complicated to solve. 
That is because it must be done consistently with the interaction between the 
matter and radiation field, in order to reach the required sub- 1% accuracy in 
Xc (see Chapter 2 for a review). In one specific example (Chapter 4), we in- 
vestigated the effect of inclusion of the higher order non-resonant two-photon 
transitions and the semi-forbidden transitions in a multi-level atom calculation, 
which was first suggested by Dubrovich & Grachev (2005) Q] using a three-level 
atom model. 

There have recently been comprehensive studies of calculations of the He i 
reombination, with all relevent radiative processes to the 0.1% accuracy level [^o[ 
^ I39II . However there still has not been a single numerical calculation which 
includes all the improvements in H i recombination (which of course has greater 
effect on the than for He). From another point of view, given the precision 
of the experimental Ci measurement, we may want to ask how accurate the 
theoretical model needs to be in order not to bias the determination of the 
cosmological parameters. In another of our projects (Chapter 6), we investigated 
how the remaining uncertainties in recombination affects the constraints on the 
cosmological parameters using Planck forecast data. We do this through use 
of the COSMOMC code [|8| , which is a numerical code for exploring the multi- 
dimensional cosmological parameter space with the Markov Chain Monte Carlo 
method. 
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1.5 Outline of the thesis 

This thesis focuses on the study of cosmological recombination and its effects on 
the CMB. Here we have briefly reviewed the standard model for the evolution 
of Universe, including the basic picture of cosmological recombination and the 
formation of the CMB. Chapter 2 provides an overview of progress in the the- 
oretical calculation of recombination, and the recent updates for obtaining the 
ionization fraction to better than 0.1% accuracy. In Chapter 3 we present 
a calculation of the spectral distortions in the CMB due to Hi Lya and the 
lowest 2s-ls line transitions, as well as the corresponding lines of Hei and Hell, 
during the epoch of recombination. Next, in Chapter 4, we investigate the ef- 
fects of including non-resonant two-photon transitions and the semi-forbidden 
transitions in the process of H i and He i recombination. Chapter 5 is a brief 
study to clarify that the previously claimed effect of the reheating of matter due 
to the distortion photons emitted during recombination is neligible. In Chat- 
per 6 we investigate how uncertainties in the recombination calculation affects 
the determination of the cosmological parameters in future CMB experiments. 
Finally we present our conclusion and ideas for future directions in Chapter 7. 
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Chapter 2 

Progress in recombination 
calculationsn 



In this chapter, we wiU give a review of progress in controUing the accuracy 
of the recombination calculation, starting from a traditional three-level atom 
model and going up to the recent multi-level atom models including interactions 
between matter and radiation. We will also describe the remaining uncertainties 
which will need to be tackled in the numerical codes in order to obtain the 
ionization fraction to better than 1%. 



2.1 Standard picture of recombination 

Cosmological recombination calculations were first performed forty years ago 
by Peebles (1968) |^ and Zeldovich, Kurt and Sunyaev (1968) Q using a 3- 
level atom model in hydrogen (H"*" -I- c^ ^ H + 7). In this simplified model, one 
only follows the detailed rates of change of electrons in the continuum, the 
first excited state and also the ground state of the atom. The higher excited 
states are assumed to be in thermal equilibrium with the first excited state. 
The cosmological recombination of H is slower than that via the Saha equation; 
it is 'Case B' recombination, since direct recombination to the ground state is 
highly prohibited and the Ly a line is optically thick. Due to the short mean free 
time of the ionizing photons compared to the expansion time of the Universe 
(by a factor of ~ 10~^), the ionizing photons emitted from direct recombination 
to the ground state easily photoionize the surrounding neutral atoms. There- 
fore, the electrons recombine mainly through the first excited state (n — 2) and 
cascade down to the ground state by the Ly a or the 2s-ls two-photon transi- 
tion. The two-photon transition plays an important role in recombination and 



the net rate is comparable to the net Lya rate (see Fig. p.5[), because only a 
tiny amount of the Ly a photons redshift out of the line and escape to infinity 
without getting absorbed or scattered. To account for the redshifting of the 
Lya resonance photons, Peebles (1968) |6^ approximated the intensity distri- 
bution as a step and scaled the Ly a rate by multiplying by the ratio of the rate 
of redshifting of photons through the line to the expansion rate of the Universe. 
The radiation field and the matter are strongly coupled through Compton scat- 
tering, and therefore the matter temperature Tm can be well approximated as 



version of this chapter will be submitted for publication; Wong W. Y. 'Progress in 
recombination calculations'. 
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the radiation temperature Tr. These two temperatures start to depart only 
in the very late stages of recombination (at z ~ 200) |2| , when most 

of the electrons have already recombined. After that the matter temperature 
decreases adiabatically, Tm oc (1 + z)^, while Tr decays as {1 + z). 

The above description gives us the standard picture for the H i recombina- 
tion. It has been argued that we should also include stimulated recombination in 
the three level atom model but the effect is quite negligible. An analogous 
physical situation was proposed for He i recombination (He"'' + e~ — > He + 7) by 



Matsuda et al. (1969, 1971) I^J, 55 , and a slower recombination than Saha 
equilibrium was then found. However, it was later argued that Hei recombina- 
tion should be well approximated by the Saha equation by taking into account 
the tiny amount of neutral hydrogen formed at the same time. Since these 
Hi atoms can capture the Hei 2^P-1^S resonant line photons as well as the 
photons from direct recombination to the l^S ground state this speeds up 
Hei recombination. This issue was not entirely cleared up until some recent 
calculations included the continuum opacity of Hi in the Hei recombination 
evolution ^ , as will be discussed in Section 2.3.4 . 

He II recombination (He^"*" -I- e~ — > He"*" -I- 7) was found to remain very close 
to Saha equilibrium |7^, |8^ due to the fast radiative rates. This, together with 
the fact that Hen recombination occurs too early to have any effects on the 
CMB anisotropics, means that we do not discuss He 11 recombination in detail 
in this chapter. 



2.2 Multi-level atom model 

Thirty years later, after the first calculations there was an increased demand 
for an accurate ionization history for modeling the CMB power spectrum for 
new experiments, for example, WMAP. Seager et al. (1999, 2000) [||, Q set a 
benchmark precision for the numerical recombination calculation by following 
the evolution of the occupation numbers of 300 atomic energy levels in H i and 
200 levels in He i without any thermal equilibrium assumption between each 
state. This multi-level Hi atom consisted of maximum 300 separated quan- 
tum number energy levels (n-states) , while the He i atom included the first 
four angular momentum states (^-states) up to n = 22 and just the separated 
n-states above that. The rate equation for each level was written down using 
the photoionization and photorecombination (bound-free) rates, the photoex- 
citation (bound-bound) rates and the collision rates. The bound-bound rates 
included all the resonant transitions but only one forbidden transition, the low- 
est spontaneous two-photon transition (2s-ls for Hi and 2^S-1^S for Hei). The 
bound-bound rate of the Lyman-series transitions were scaled with the Sobolev 
escape probability Ps |^ to account for the redshifting and trapping of the dis- 
tortion photons in the radiation field; this reduces to Peebles' step method when 
Ps oc 1/t (see Equation (3.14) for an explicit expression for ps) and t :i> 1 (r is 
the optical depth of the line). All the rates were obtained with the radiation 
background approximated as a perfect blackbody spectrum. 
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In this multi-level atom model, Seager et al. (1999, 2000) @ |7|l fomid 
a speed-up in H i recombination at low redshift compared with the standard 
Case B recombination Q at low redshift, and a delayed Hei recombina- 
tion in contrast to that from the Saha equation. The Hi recombination was 
faster than previously estimated because of the non-zero bound-bound rates 
among higher excited (n > 2) states. In the three-level atom model, the higher 
excited states are assumed to be in thermal equilibrium and the bound-bound 
rates between these states are negligible. However, these bound-bound rates are 
actually dominated by spontaneous de-excitations due to the strong but cool 
radiation field, which means that the electrons prefer to cascade down to the 
lower energy states rather than staying at the higher excited states. This results 
in faster H i recombination. 

On the other hand in this study, the He i was found to follow a standard 
hydrogen-like Case B recombination, agreeing with the earlier study of Mat- 
suda et al. (1969, 1971) |5|, which adopted a three-level atom model for 
He I by considering the singlets only. There are two sets of states in the He i 
atom; the singlets and the triplets. In this multi- level atom model calculation, 
the triplet states were found to be highly unpopulated because the collisional 
transitions between singlets and triplets are weak, and therefore the electrons 
mainly cascade down to the ground state via the singlet states. Concerning any 
mechanisms which might bring the He i recombination into Saha equilibrium, 
Seager et al. (2000) [Q found that the photoionization rate of H i was much 
lower than the Hei PS-2^P photoexcitation rate, and concluded that the Hi 
atoms have negligible effect on stealing the He i resonance line photons in order 
to speed up the He i recombination. 

Evolution of the matter temperature with all the relevant cooling processes 
(specifically, Compton, adiabatic cooling, free-free, photorecombination and line 
cooling) and the formation of hydrogen molecules (for example, H2) was also 
considered, but these effects, along with the collisional transitions, were found to 
be negligible for the ionization fraction Xc (see Table 2.1 for the magnitude of the 
change in Xc)- Seager et al. (2000) j7|] also discussed the effects of the secondary 
distortions due to photons emitted from the H i Ly a and 2s-ls transitions, 
and also from the corresponding transitions in Hei and Hen. These distortion 
photons can be redshifted into a frequency range where they could photoionize 
the electrons in the first excited (n = 2) state and the ground state of H i during 
the time of Hi recombination. Again, the effect, was found to be very small. 

In order to reproduce the accurate numerical results without going through 
the full multi-level calculation, the authors used a so-called 'effective three- 
level model' [72I by multiplying by a 'fudge factor' Fji the recombination and 
ionization rates in the standard Case B recombination calculation to reproduce 
the speed-up of H i recombination. For He i, the result can be well approximated 
by considering the standard Case B recombination situation in a three-level 
atom model (singlets only), consisting of the continuum, the first excited singlet 
state and the ground state (see Section 3.2 for details), recfast is the 
publicly available computer code which calculates the ionization fraction Xe 
(the detailed profile of the last scattering surface) using the effective three-level 
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model discussed above. It is currently adopted in most of the commonly used 
Boltzmann codes, for example, CMBFAStQ, CAMB Q and CMBEASy[|i^, to 
numericalUy evolve the accurate CMB anisotropy spectrum for different sets of 
cosmological parameters. 



2.3 Recent improvements and suggested 
modifications 

Recently, driven mainly by the up-coming high-£ CMB experiments ^ and 
also the possibility of detecting spectral distortions ^ |6| , there have 
been many suggested updates and improvements to the the multi-level atom 
model suggested by Seager et al. (2000) Q. In this section, we discuss these 
new physical processes included in recombination, concentrated on those which 
may lead to more than 1% level change in the ionization fraction x^- 



2.3.1 Energy levels 

In the Seager et al. (2000) ||7^ calculation, they only considered separated n- 
states for Hi, while the /-states are assumed to be in thermal equilibrium within 
each rt-shell. Rubino-Marti'n et al. (2006) first tried to relax this assumption 
by resolving all the /-states up to n — 30, and later the same authors]^ even 
pushed the maximum level to n = 100. They found that the total population 
of any shell is smaller than the value obtained from the Saha equation during 
H I recombination. The deviation of the populations from Saha was claimed to 
increase from ^ 0.1% at z = 1300 to ~ 10% at z = 800 and this in general led to 
a slower recombination at lower redshift compared with previous studies (see 



Table 2.1). There seems to be no need to further consider the separate states 
with different spin orientations (for example, hyperfine splitting), since the rates 
of the resonant transitions connecting individual hyperfine splitting states (or 
only one /-state to the other states) are the same, even if these splitting states 
are not in equilibrium. 

Note that there is a serious problem in the Rubiho-Martfn et al. (2006) 
model: the ionization fraction Xc does not converge when increasing numbers 
of 71-shells are included. Due to computational limitations, the most intensive 
calculation involved 5050 separate /-states with maximum n=100. Although 
the /-changing and n-changing coUisional transitions were additionally consid- 
ered in their model, these authors found that these transitions were not strong 
enough to bring the higher excited states back into thermal equilibrium and the 
divergence problem remained. Although there is undoubtedly still some physics 
missing in this model (which will be discussed later in this chapter), it is still 
worth asking how many levels we need to consider to solve for the recombination 
of H I atom. Using the thermal equilibrium assumption in each n-shell, Seager 
et al. (2000) |Q found that the ionization fraction Xc converges well when con- 
sidering a maximum of 300 energy levels, and claimed that this should be the 
maximum number of levels which needs to be considered by arguing that for 
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such a large n state the thermal broadening width of the level is larger than the 
gap between that level and continuum. 

For He i recombination, no such convergence problem exists. Switzer & 
Hirata (2008) Q performed a similar multi-level atom model calculation 
including the interaction of matter with the radiation field by resolving all the 
^-states with n < 10. The number of resolved /-states are limited by the avail- 
ability of the atomic data of He i and so this calculation is the best that can be 
carried out for now (the limited availablility of the atomic data will be discussed 



in Section 2.3.5). Switzer & Hirata (2008) IgSj found that the change of Xc was 
smaller than 0.004% when reducing the maximum principal number n of the 
levels from 100 to 45. In their model, the effect was not significant because the 
feedback of the spectral distortion from these highly excited states su ppres sed 
the net recombination to the ground state via these states (See Section ^.3. 4 for 
more details). 

2.3.2 Bound-bound transitions 

In the formerly standard recombination model, the lowest two-photon transition 
is the only forbidden transition considered. Dubrovich & Grachev (2005) |^ 
first suggested that it might be important to include more intercombination (i.e. 
transitions connecting triplets and singlets), the Hci 2'^Pi-l^So spin-forbidden 
transition specifically and the non-resonant two-photon transitions from the 
higher excited states (ns, nd Is for Hi and n^S, n^D— > l^Sg for Hei). They 
demonstrated that these additional transitions significantly speed up the re- 
combination in an effective three-level atom model calculation. We first focus 
on the effect of including the intercombination 2'^Pi-l^So transition. With 
the addition of this transition in a standard multi-level atom model, Wong & 
Scott (2007) (see also Chapter 4) found that more than 40% of the photons 
from n = 2 state cascaded down to the ground state through the triplet 2'^Pi 



state (see Table 4T). This is almost the same as the amount of electrons going 



from 2"'^Pi through the resonant transition, and the net rates of the 2^Pi-l^So 



and 2^Pi-l^So transitions are comparable (see Figure 4.8 and Section p~3| for de- 
tails) under the Sobolev photon escape approximation. The He i recombination 
speeds up due to this extra channel through the triplets to the ground state and 
the change on Xc is about 1.1% at z ~ 1750. Switzer & Hirata (2008) also 
found that the 2'^Pi-l^So transition is important for Hei recombination in their 
improved multi-level model calculation including the evolution of the radiation 
field; the effect of the radiative feedback between the 2^Pi-l^So and 2^Pi-l^So 
transitions was found to bring a 1.5% change in Xc (the details of the feedback 



effect will be discussed in Section 2.3.4). 



For the higher order two-photon transitions, Dubrovich & Grachev (2005) |2^ 
attempted to include the corresponding rates in an analogous way to the lowest 
2s-ls two-photon transition. However, they found that these two-photon tran- 
sitions from high n states are more complicated than the 2s-ls one. That is 
because the matrix elements for these transition rates have poles when the inter- 
mediate states are not virtual (i.e. ns, nd — > mp — > Is with 1 < m < n) and this 
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is in-distinguishable from the resonant one-photon transitions themselves. If we 
include all the poles in calculating the two-photon decay rate, then we obtain 
a very fast rate since the process is dominated by the resonant Lyman-series 
transitions. Additionally we double count the number of electrons recombining 
through those one-photon resonance transitions. In order to avoid these prob- 
lems due to the resonance poles, Dubrovich & Grachev (2005) Q approximated 
the non-resonant two-photon rate by considering only one pole (the np state) 
as the intermediate state in the matrix element. Their estimated rate was very 
fast (scaling as n for large n) and this dramatically sped up the recombination 
process, with Axc equal to a few percent. 

Wong & Scott (2007) proposed an improved, net non-rcsonant two- 
photon rate for Hi from rt = 3[^3[ and this was significantly lower. This 
calculation included all the non-resonant poles (i.e. all the n > 3 intermedi- 
ate states). By comparing with the Dubrovich & Grachev (2005) estimate 
at n = 3, the rate obtained is an order of magnitude smaller, due to the de- 
structive interference of some matrix elements, which was ignored in Dubrovich 
& Grachev (2005) (since they only considered one pole). Using this rate 
with the same n scaling given by Dubrovich & Grachev (2005) for the higher 
n two-photon rates, Wong & Scott (2007) |85| found that the maximum change 
in Xc was only 0.4%. Chluba & Sunyaev (2007) performed a more detailed 
calculation of the high n two-photon rates for Hi by studying the frequency 
distribution profile of the photons from these transitions. They estimated the 
effective two-photon rates by subtracting Lorentz profiles of the possible reso- 
nant transitions directly from the full two-photon profile in order to avoid the 
double-counting from the one-photon resonant transitions. The rates they found 
were lower than the ones given by Dubrovich & Grachev (2005) |Q, also due 
to destructive interference of the matrix elements. With their effective rates, 
Chluba & Sunyaev (2007) obtained essentially the same value for the maxi- 
mum change in at similar redshift range as found by Wong & Scott (2007) |^ . 

Hirata & Switzer (2008) Q and Hirata (2008) @ further studied the role 
of these high n two-photon transitions in He i and H i recombination, respec- 
tively, by including the related two-photon scattering (Raman scattering) and 
the possibility of re-absorption of the photons from the resonant intermediate 
states. In their model, they separated the spectrum of the photons into non- 
resonant (photons emitted through a virtual intermediate state) and resonant 
regions. They added an additional rate due to these non-resonant photons in 
analogy to the lowest 2s-ls two-photon rate. The higher order non-resonant 
two-photon rates were also found to be much lower than those estimated by 
Dubrovich & Grachev (2005) again because of the destructive interference 
of the matrix elements, and the rates scale as n~^. The resonant transitions were 
considered as being photons from the corresponding one-photon resonant tran- 
sitions, but with a modified line profile; these photons were highly probable to 
be scattered or absorbed by other atoms. For He i, Hirata & Switzer (2007) |Q 
found that inclusion of these higher order two-photon transitions brought no 
more than a 0.04% change in Xc- But for Hi, with the additional consideration 
of the feedback between the Lya line and the two-photon transitions p3] (see 
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Section 2.3.4 for details), the change in Xc was found to be more than a percent 
around the peak of the visibihty function. 

Some other forbidden transitions were also included in He i recombination, 
specifically, the magnetic dipole 2'^S-l^So transition ^ the electric 
dipole transitions with n < 10 and I < 7 |p5[, the intercombination n^Pi-l-'^So 
transitions, the electric quadrupole n^D-PSo ^ 4) transitions, the magnetic 
quadrupole 2'^P2-l^So transition and the electric octupole n'^F-l^So {n > 4) 
transitions iQ. However, the effect of the inclusion of all the above transitions 
is very small (Axe < 0.001%) and can therefore be neglected. One may ask 
whether we should include the one-photon 2s-ls magnetic dipole transition for 
H I. The rate of this transition is equal to 2.49 x 10~^ s^^|§, which is about 6 
orders of magnitude smaller than the 2s-ls two-photon transition. Therefore we 
expect that the effect of the inclusion of this magnetic dipole transition should 
be neghgible. 

It is worth remembering that there are two electrons bound in each He i 
atom. In the standard He I recombination calculation, we usually consider the 
inner electron to be in the ground state. One may wonder whether these other 
electrons might sometimes leave the ground state by stealing photons and get- 
ting excited to higher levels. However, He ii recombination occurs much earlier 
than He i recombination, and therefore almost all of the inner electrons were 
already in the ground state based on the Boltzmann distribution at the time 
when He i recombination began. In order to excite the inner electrons from 
the ground state to the first excited state, the energy of the incident photons 
would need to be about 40 eV, which is almost double the ionization energy 
of He I. This means that the abundance of such 40 eV photons is 10^ of the 
He I ionization photons at z — 2500, based on the blackbody spectrum, implying 
that the inner electrons have almost no chance to get excited from the ground 
state during the He i recombination. Hence we can completely neglect all such 
transitions. 



2.3.3 Bound-free transitions 

One of the approximations adopted in the standard recombination calcula- 
tion is that there are no direct recombinations to the ground state. This is 
because the photons emitted in this transition immediately reionize another 
neutral atom (the same situation applies to both Hi and Hei). Chluba and 
Sunyaev (2007) jTo| revisited this approximation by calculating the net rate of 
direct recombinations to the ground state for H i through detailed consideration 
of photon escape. Although the escape probability of a photon emitted from 
the continuum to the ground state is about 10-100 times larger than that of the 
Lya photons, the inclusion of these direct recombinations only brings about a 
0.0006% change in Xc- Hu et al. (1995) had earher argued that the direct re- 
combination of He I should be possible, due to the absorption of the continuum 
photons by the tiny amount of H i atoms in the later stages of He i recombi- 
nation. However, Switzer & Hirata (2008) showed that this effect on the 
speed-up of He i recombination is neghgible (Axe — 0.02%) by calculating the 
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effective cross-section of tlie bound-free transition to the ground state due to 
the presence of H i. From the above dissension, we can therefore safely neglect 
direct recombinations to the ground state for both H i and He i. 

2.3.4 Radiative transfer 

In the standard multi-level calculation of recombination, the radiation back- 
ground field is approximated as a perfect blackbody spectrum. For the inter- 
action between atoms and the radiation field, the Sobolev approximation is 
adopted to account for the escape probability of the photons redshifting out of 
the line. But in order to calculate to better than the 1% level, the above 
approximations are not sufhcient, and fundamentally we need to solve for the 
evolution of the number densities of the atomic levels and the radiation field, 
with the distortion photons from recombination process solved consistently in an 
expanding environment. Several recent studies M, nn, pi 
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have suggested that additional radiative transfer processes (for example, the 
feedback between lines) might cause significant effects on recombination. In par- 
ticular, Switzer & Hirata (2008) |3^, |8^, |8^ have performed the most complete 
and systematic multi-level He i atom model calculation, with the consideration 
of both coherent and incoherent scattering process between atoms and photons. 
They specifically included the feedback between lines, absorption due to the 
continuum opacity of Hi, stimulated and induced two-photon transitions, the 
collisional transitions and Thomson scattering (all examples of incoherent scat- 
tering) , together with partial redistribution of the line profile due to coherent 
scattering. We will discuss each of these processes in turn. 

Feedback from spectral distortions 

In the standard multi-level atom calculation, no feedback between resonant lines 
is considered. However, in practice distortion photons escaping from the higher 
order resonance transitions will redshift to a lower line frequency and excite 
electrons in the corresponding state. For example, photons emitted from Ly7 
transitions can excite electrons in the ground state after redshifting to Ly (3 or 
Ly a line frequencies. In general, this feedback process will suppress the net 
recombination rate to the ground state thereby slowing down recombination. 
Switzer & Hirata (2008) |Q used an iterative method to include the feedback 
between transitions connecting the excited states and the ground state during 
He I recombination. They only considered the radiation being transported from 
the next higher transition [{i + l)th state to l^So] to the ith transition (to the 
ground state in the same species) . They found that the most significant change 
to the ionization fraction {Axe = 1.5%) is due to the feedback between the 2^^?- 



l^So and 2^Pi-l^So transitions. Chluba & SunyaevjlO also studied the same 
feedback effects among the Lyman-series transitions during Hi recombination. 
They found that feedback from the Ly /3 transition on the Ly a line accounts for 
most of the contribution, and the maximum change in x^ is about 0.35%, this 
appearing to be a convergent result when including Lyman-series transitions up 
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to 71 = 30. For H i recombination, we also need to consider the distortion photons 
from He i recombination feeding back to the H i line transitions (especially the 
Lyman series), which brings about a 0.1% change in Xe jSOj. In the Seager 
et al. (1999) recombination model, only the photons from Hei 2^Pi-l^So 
and 2^So-l^So transitions were considered as secondary distortions on the Hi 
recombination. This should clearly be extended by calculating a detailed Hei 
line spectrum, including all the released photons. 

Stimulated and induced two-photon transitions 

The standard recombination model only includes the spontaneous 2s-ls two- 
photon emission rate and the corresponding two-photon absorption rate coming 
from detailed balance. Taking Hi as an example, the spontaneous two-photon 
decay is 

H(2s) ^ H(ls) -I- 7spon + 7spo„ , (2.1) 
and the two-photon excitation is 

H(ls) + 7bb + 7bb H(2s) . (2.2) 

Here 7spon represents a spontaneously emitted photon and 7bb represents a pho- 
ton taken from a blackbody radiation spectrum. Chluba & Sunyaev (2005) 
suggested that one should include the stimulated Hi 2s-ls two-photon emis- 
sion due mainly to the low frequency background photons. The two stimulated 
decays are 

H(2s) ^ H(ls) + (2.3) 

and 

H(2S) H(ls) + 7stim + 7stim , (2.4) 

where 7stim refers to a photon from stimulated emission. The recombination 
is found to speed up, and these authors claimed that the effect can be more 
than 1% in x^- Later, Kholupenko & Ivanchik (2006) ^l| pointed out that 
the induced Hi 2s-ls two-photon absorption of a thermal background photon 
and a redshifted distortion photon from the H i Ly a transition should also be 
considered, i.e. 

H(ls) + 7bb + 7dist ^ H(2s) , (2.5) 

where 7dist represents a spectral distortion photon. By including this absorption 
process, recombination is actually delayed overall, and the maximum change in 
a;c is about 0.6%[|4[ |l|. Hirata (2008) |||| extended the above ideas further 
to include the higher order two-photon transitions (H i nd, ns-ls) using the 
steady-state approximation. Instead of adopting an effective rate, he performed 
a radiative transfer calculation to account for the emitted line photons, whether 
they are being re-absorbed or scattered later. The result showed that the recom- 
bination speeds up after inclusion of the stimulated and induced higher order 
two-photon transitions, the maximum change being 1.7% in at 2; ~ 1250, 
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which is bigger than the resuh of using only the effective rates in the previ- 
ous studies Hirata (2008) [|4| also investigated the effect of two other 
relevant two-photon process: Raman scattering 

H(nO+7^H(ls)+7', (2.6) 

where 7' is a photon with higher energy compared with 7; and direct two-photon 
recombination to the ground state 

H++e- ->H(ls)+7 + 7'. (2.7) 

The direct two-photon recombination process was found to be negligible, but on 
the other hand, the Raman scattering brought about dramatic effects on Hi re- 
combination. Raman scattering is dominant in the 2s-f s transition, since the 2s 
state is the most populated among all the ns and nd states with n>2. Through 
Raman scattering, the CMB photons can excite atoms in the 2s state and the 
atoms will decay down to the ground state by emitting photons with frequencies 
between the Ly/3 and Lya lines. Therefore, Raman scattering provides another 
channel for the electrons to get down to the ground state and this initially speeds 
up recombination. However, the photons emitted from the Raman scattering 
process having energy larger than Ly a will redshift and feed back on the Ly a 
and 2s-ls transitions. This additional feedback delays recombination and Xe 
increases by about 1% at z ~ 900. 

For He i recombination, a similar study was performed by Hirata & Switzer 
(2008) and a much smaller effect was found on Xc (< 0.01%). The reason is 
that the abundance of H is much greater than for He (about a factor of 12 in 
number) which leads to a lower optical thickness in the case of the He i 2^Pi-l^So 
line than the Hi Lya line for the resonant photons from two-photon transitions. 
The other reason comes from the different shapes of the frequency spectra of 
the lowest two-photon transition at low frequencies. The frequency spectrum 
for Hei 2-'^Pi-l^So is proportional to i'^, while that at Hi 2s-ls is proportional 
to I/. This is because the H i 2p and 2s states are essentially degenerate (actually 
the 2p state is slightly lower than 2s due to the Lamb shift [Q, but the shift 
is only 4.372 x 10~^ eV) and so there is a pole in the matrix element at zero 
frequency when 2p is the intermediate state | |3^ . The stimulated and induced 
two-photon transitions dominate at low frequencies and therefore there is a 
larger probability in the H i two-photon spectrum at both ends (where one of 
the photons has a small frequency) . As a result the effect is more significant in 
H I recombination. 



Photon absorption due to continuum opacity of H i 



The other important improvement in He i recombination is inclusion of the con- 
tinuum opacity of H 1 1^ . In the later stages of He i recombination a tiny but sig- 
nificant amount of neutral hydrogen H i is formed (nni/j^H < 10"'* at z ~ 2000), 
and these H i atoms can absorb (through photoionization) the distortion photons 
emitted during Hei recombination. In Section 2.3.3, wc have already discussed 
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how the effect on the direct recombination of He i due to this continuum opacity 
of Hi is neghgible. However, the presence of the Hi continuum opacity signifi- 
cantly affects the transitions connecting the excited states and the ground state, 
particularly the 2-'^Pi-l-'^So transition. The 2-'^Pi-l^So transition, which is the 
lowest He i resonance transition, is also one of the main paths for the electrons 
to cascade down to the ground state. In the standard multi-level atom model, 
about 60% of the electrons in the n = 2 state reach the ground state through 



this transition (see Table 4.1). The energy of the photons emitted from the 
2^Pi-l^So transition is equal to 21.2 eV, which is much larger than the ioniza- 
tion energy of Hi. Therefore, the Hi atoms can absorb: (1) the Hei 2^Pi-l^So 
line photons directly; or (2) the redshifted line photons from the next higher 
transitions before they redshift down to the 2^Pi-l^So line and excite another 
atom. This process removes these distortion photons and prevents them from 
re-exciting other He i atoms: 

He(2ip) ^ He(liS)+7 
H(ls)+7 ^ H++C-. (2.8) 

For process (1), the usual Sobolev escape probability can be modified due to 
the direct line photon absorption by the H i atoms instead of He i |4l| [s^ . The 
modified escape probability, which is in general larger than the Sobolev value, 
has been applied to the Hei resonant n^P-l^So, intercombination n^D-l^So 
and quadrupole n-'^F-l^So lines Recombination is significantly sped up 

mainly due to the extra continuum opacity of Hi within the 2^Pi-l^So line. 
This effect gives more than a 2% change in Xo, while the opacity in other lines 
only contributes about 0.05%. For process (2), the absorption of the redshifted 
line photons suppresses feedback between the lines. For example, there are some 
distortion photons from Hei 2'^P-l^So which are absorbed by Hi before they 
can redshift down to the Hei 2^Pi-l^So line frequency to excite electrons in 
the ground state of He i atoms. Therefore, the number of redshifted distortion 
photons available for the feedback between Hei lines is smaller, and hence the 
Hei recombination speeds up a little. Overall, the continuum opacity of Hi 
modified to include these feedback process brings about a 0.5% change in a^e- 



Coherent scattering 

In the Sobolev escape probability method, a Voigt profile is assumed for both 
the frequency spectra of the emitted and absorbed photons in the line transi- 
tions. However, this will only be true when the system is very close to thermal 
equilibrium. For an optically thick line (for example, Hei 2^Pi-l^So) without 
the continuum opacity of other species of atoms, the radiation field in the region 
of the line frequency is in thermal equilibrium with the population ratio of the 
corresponding two levels relevant for this transition. However, in the presence 
of the continuum opacity of H i, the H I and He i atoms complete for the the 
distortion photons from the line transitions and so no such thermal equilibrium 
exists. The emission and absorption line profiles may not be the same as each 
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other or equal to a Voigt function, since there is no complete redistribution in 
the line profile. In such a non-equilibrium situation, we need to consider all the 
possible paths for an electron at each state to go after it is excited by a resonant 
photon from a lower state. Therefore, besides the incoherent scattering pro- 
cesses, we also need to consider coherent scattering (relative to the atom's rest 
frame). An electron excited by a resonant photon to a higher state can decay to 
the original lower state by emitting a photon with the same energy without any 
intermediate interaction. The emitted and absorbed photons have no energy 
difference in the atom's rest frame, but there is a small fractional change in 
the photon's frequency (on the order of v/c, where v is the atomic velocity) in 
the comoving frame. If the effects of coherent scattering are significant, the line 
profile is only partially redistributed and the frequency spectrum of the emission 
line photons depends on the radiation background. Switzer & Hirata (2008) |^ 
performed a Monte Carlo simulation for the partial redistribution of the profile 
in the Hei rt^Pi-l^So resonance line due to coherent scattering. The effect they 
found was about 0.02% in Xc, compared with the model having feedback and 
continuous opacity of H i, as discussed above. 

Thomson scattering and collisional transitions 

Thomson scattering and collisional transitions were also considered by Switzer 
& Hirata (2008) ||8^ in the He i recombination calculation, but both of these 
processes were found to be negligible. During Hei recombination, Thomson 
scattering may be significant, since a large fraction of electrons have not yet 
recombined. The photons can gain energy after multiple electron scatterings 
and the photons which had previously redshifted out of the line can be scattered 
back into the line. This reduces the escape probability of the line and hence 
delays the recombination. However, in the presence of feedback between lines 
and the continuum opacity of H i, the distortion photons are more likely to get 
re-absorbed instead, and therefore Thomson scattering is strongly suppressed. 
The net effect becomes only 0.03% in Xc- During Hi recombination, Thomson 
scattering should also be considered, because the optical depth of the Lyman- 
series lines is very high {r-^ 10^ for Lya, which is 10'^ times that of the He i 2^Pi- 
l^So hue). Due to this high optical depth, a similar calculation is necessary 
for studying the partial redistribution of the Lyman-series line profiles with all 
the possible coherent and incoherent scattering processes. However, no such 
systematic calculation (similar to the He i one) has been performed yet. Since 
the rate of H i recombination is mainly controlled by the trapping of the Ly a 
photons, there are several studies concerning only the line profile of the Lya 
transition. Rybicki & Dell'Antonio |Q studied the time-dependent spectral 
profile of the Lya transition in an expanding environment using the Fokker- 
Planck equation and, found that the quasi-static assumption is an adequate 
approximation for this transition. Several other works 
included the effect of the frequency shift of Ly a due to the recoil of the H atoms, 
with suggestions that the effect on Xe may be at the level of 1%. 

Collisional transitions, caused by the collisions between atoms and ions, are 



, 31, Hgl have also 
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usually neglected in recombination, because of the high photon to baryon ratio 
(~ 10^). Such coUisional processes tend to bring the species into equilibrium and 
maintain statistical balance between the energy levels. The bound-free transi- 
tions, bound-bound transitions, and charge exchange (He+ + H ^ He +H+ + 7) 
between H i and He i due to the collisions are found to be too slow to have any 
effect on Hei recombination js^j . In the later stage of Hi recombination, the 
separated Z-states fall out of equilibrium and the collisional transitions become 
very important for redistributing the electrons within each n shell, at least for 
the higher excited states (n > 50) For the lower excited states, radiative 
transitions are dominant and the effect of collisional processes between the H I 
2s and 2p states was found to be negligible . 

The electrons, ions and neutral hydrogen are well approximated as a single 
tightly coupled component in the standard recombination picture ||3^ . They are 
considered as a single 'baryon' fluid and described by a single temperature, the 
matter temperature Tm- The matter temperature is very close to the radiation 
temperature Tr, during recombination, due to the strong effects of Compton 
scattering. During Hei recombination (2:>1600), for example, the fractional 
temperature difference (Tr — Tm)/Tr is smaller than 10^^ The effects 

of the Compton scattering become weaker during H i recombination, since most 
of the electrons are captured to form neutral neutral atoms. Adiabatic cooling 
starts to become important for matter when Compton scattering effects become 
slow compared with expansion time. The matter temperature then starts to 
depart from the radiation temperature, because the matter cools faster. But 
actually, even during Hi recombination (700 < z < 1500), the fractional differ- 
ence between these two temperatures is no more than 1% (see Fig. 2 in ||8^). 
This summarizes the general picture for the evolution of matter temperature. 
Several authors ^ |2| have performed detailed calculations of the evolu- 
tion of Tm by including all the relevant heating and cooling processes between 
the matter and radiation fields, in addition to Compton and adiabatic cooling. 
The results found are basically the same as in previous studies Q , with the 
additional processes bringing negligible change on the matter temperature. One 
study suggested that one should include the heating of matter due to the dis- 
tortion photons emitted during H i recombination, and that this effect delayed 
recombination |]5l| . However, the coupling between matter and these distortion 
photons is very weak. Almost all of these photons go into the radiation field 
and form the spectral distortion lines on the CMB blackbody spectrum (see 
Chapter^ for details). This additional suggested effect is therefore negligible. 

2.3.5 Atomic data 

In the multi-level atom calculation of cosmological recombination, the non- 
equilibrium situation existing between states is important, since radiative pro- 
cesses are much stronger than collisional ones |3^, . Based on recent studies, 
it is necessary to include energy levels with principal quantum number n < 50 
for Hei and n < 300 for Hi in the multi-level atom model. Therefore, for the 
numerical recombination calculation, detailed and accurate atomic data are re- 
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quired for the energies of the states, and for the bound-free and bound-bound 
transition rates, not only for the lower states, but also for the higher excited 
states. 

For the H i atom, there is an exact solution for the non-relativistic Schrodinger 
equation, and the energies of each (n, I) state are given by En = —Rn/n^, where 
i?H is the hydrogen Rydberg constant and hpcRn = 13.5984 eV Q. With the 
exact wavefunctions, the rates of the bound-bound resonant (electric dipole) 
transitions between resolved Z-states can also be determined to very high ac- 
curacy . This is also true for the lowest two-photon 2s-ls transition, and 
there are many papers in the literature determining the theoretical value of this 
spontaneous rate, |, 1^, |9[ |o[ 0, 1?!. The latest value of A" 

is 8.2206 s~^, given by Labzowsky et al. (2005) |46 , and this agrees with other 
calculations to about the 0.1% level of accuracy. This small uncertainty has 
negligible effect on recombination. For two-photon transitions from the higher 
excited states (n > 2) to the ground state, we need to have the detailed spectra 
of the emitted photons in order to avoid double counting the photons in the res- 
onant transitions. By direct summation of the matrix elements or by by using 
the Green functions method, the spectra can be calculated to 0.1% accurarcy 



(see 1 12, ^ and references therein). 

For n < 10, TOPbase provides spectra for the photoionization cross 
sections a{iy) for each (n, I) level. And we can use the Gaunt factor approx- 
imation |5^, ^ to calculate the photoionization cross section for the states 
with n > 10. The Gaunt factor is the ratio of the photoionization cross- 
section from a quantum-mecahnical calculation to the value obtained from the 
semi-classical electromagnetism formalism (see, for example. Chapter 6 i n p^ ). 
Rubiho-Martfn et al. (2006) p7| compared three numer ical methods |, |, |gfor 
obtaining these cross-sections and found that the results agree to the percent 
level. 

The atomic physics of He i is more complicated than H i because it is a two- 
electron system. Morton et al. (2006) |Q have provided the largest and most 
recent set of ionization energies of resolved I states for n < 10 and I < 7, with 
accurarcy better than 10~^, combined with both experimental and theoretical 
results. For the other states, it is usual to adopt re-scaled hydrogenic values; 
it should be a good approximation to consider an electron orbiting a point- 
like He+ ion for l>2. For the bound-bound transition rates, Drake & Mor- 
ton (2007) have also presented the most up-to-date data-set of the emission 
oscillator strengths / for the electric dipole transitions, and also the intercombi- 
nation (spin-forbidden) transitions between the singlets and triplets with n < 10 
and l<7. Bauman et al. (2005) |^ ^ developed a computer code for generat- 
ing fs and the Einstein coefficients A for the bound-bound transitions for even 
higher excited states (n<13 and l<ll) by combining different data sources 
and approximations. The accuracy of these two approaches is about 5% to 
10%, respectively, which is estimated by comparing the results with experimen- 
tal data using the adopted approximations ^ . For the higher order reso- 
nant transitions, n> 12, the rescaled hydrogenic values are used for the bound- 
bound resonant rates and the uncertainty should be at least at the 10% level. 
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Since the resonant lines are optically thick, the intercombination transitions 
(n^Pi-l-'^So with n > 2) play an important role in recombination, especially 
the 2^Pi-l^So transition. The theoretical value of the spontaneous 
transition rate ^23Pi-iiSi ranges from 171 s~^ to 233 s^^ in different calcula- 
tions (l|, 1^, |9[||. The latest Azsp^.iigi value is ITTs"! given by Lach & 
Pachucki (2001) p^. Although variations of among estimates of this rate are 
about 30%, the effect on Xc is only at the 0.1% level . The most important 
forbidden transition for Hei is the lowest two-photon 2^o~l^So transition. The 
latest value of the spontaneous rate A^i''g^_^ig^ is 51.02 s~-^ which agrees 
with other theoretical values |l| at the 1% level. For the higher order two- 
photon transitions (n^Sg, n^D2-l^So), Hirata & Switzer (2008) |^ have tried to 
estimate the corresponding rates and also the freqency spectrum of the emitted 
photons by direct summation of the matrix elements. The accuracy of their 
method is at about the 10% level. But this uncertainty brings almost no change 
on recombination, since the effect of the inclusion of the higher order two-photon 
transitions was found to be insignificant for He i recombination JM . 

For the bound-free cross-sections, Hummer & Storey (1998) |38| provided the 
largest set of data for the spectrum of the cross-section with n < 25 and 
I < 4, while Topbase |l^ only contains a{h') with n < 10 and I < 2. Bauman et 
al. (2005) ^ have combined these two results with other approximations in 
a computer code which can generate cr(i^) up to n = 27 and I = 26. These three 
sets of data (although not entirely independent) agree at the few percent level. 
For higher excited states {n > 10), the re-scaled hydrogenic cross-section [j79[ 
can be used. This is a reasonable approximation, giving accuracy at about the 
10% level [|8j. Overall there is about a 10% error in the atomic data of Hei, but 
the effect on Xc should be no more than the 0.1% level, helped considerably by 
the low abundance of Hei (about 8% of the total number of H and He atoms). 

2.3.6 Fundamental constants, cosmological parameters 
and other uncertainties 

The accuracy of the available fundamental physical constants is of course im- 
portant for the numerical recombination calculation. The biggest uncertainty 
comes from the gravitational constant G , due to the inconsistency among 
different experimental measurements (see Chapter 10 in for details). The 
latest recommended value by the Committee on Data for Science and Technol- 
ogy (CODATA) is G = 6.67428(67) x 10"" m3kg-ls"^ with a fractional uncer- 
tainty equal to 10^"* [Q. The gravitational constant mainly affects the overall 
time scale of the expanding Universe. However, this uncertainty brings almost 
no effect on x^ {Ax^ <C 10^'^) All other relevant physical constants are 

measured to much higher accuracy and their effects on recombination can be 
ignored. 

The CMB monopole temperature Tcmb is one of the few cosmological param- 
eters that can be measured directly by experiments, and is usually considered 
as one of the fundamental 'input' parameters in the standard six parameter 
ACDM cosmological model for calculating the CMB temperature and polariza- 
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tion anisotropies. Given TcmB: we can determine the radiation density or the 
photon background field of the Universe, and this strongly affects the speed of 
recombination. The latest value of Tcmb is 2.725 ± 0.001 K, which is the fi- 
nal assessment, including calibration and other systematic effects, coming from 
measurements made with FIRAS instrument (on the COBE satellite) |2^ . Al- 
though the relative uncertainty in AT/T is only 0.04% , it leads to a 0.5% 
change in at z~900. But the corresponding effect on the Ci is only at 

the 0.1% level for €~ 2500 ]ll|, ||. 

The other uncertainty among the input cosmological parameters is the pri- 
mordial helium abundance Yp (defined to be the mass fraction of helium) [|ll| . 
In the standard Big Bang Nucleosynthesis (BBN) calculation, the derived value 
of Yp only depends on the baryon to photon ratio tjiq = 10^° {n-g/n^), and can be 



numerically calculated to about the 0.2% level of accuracy |78|. Note that the 
number of neutrino species N^, is assumed to be 3 in standard BBN (although 
it is not quite correct if there is mixing between different kinds of neutrinos). 
The number of neutrino species affects the He abundance because of the change 
in the expansion rate of the Universe (AFp w 0.013 AA'y) j?^. Based on stan- 
dard BBN and the WMAP five-year results, Yp is determined to be equal to 
0.2486 ±0.0005 [|2^, which is a little larger than the value estimated from the di- 
rect observational results Yp = 0.240 ± 0.006 1|7||. After the BBN epoch, helium 
can be produced in all H-burning stars, while some other heavier elements, such 
as oxygen O, are produced only in short-lived massive stars. In low-metallicity 
regions, the measured He abundance should be close to Yp if the oxygen to 
hydrogen ratio 0/H is very low. Therefore, the observed value of Yp is usually 
determined by studying line emission from the recombination of ionized H and 
He in low-metallicity extragalactic Hii regions. However, the observed value of 
Yp is still quite uncertain, due to the sysmatic errors and the lack of evidence 
for the correlation between helium and oxygen abundances (see Section 3.3 in 
[ [78| for details). Due to discrepancies between the theoretical and observational 
results, the uncertainty of Yp should be considered to be about 5% and this 
brings a change in at about the 1% level at redshifts around the peak of the 
visibility function. 

In most recombination codes, only the masses of the constituents of atomic 
hydrogen and helium are taken into account for converting the baryon density 
r^B to the number of hydrogen atoms nu, i.e. 

nH^^i^. (2.9) 
SttG Tor 

It has been argued that we should also consider the binding energy in each 
atom fr\ as well as the abundance of lithium in the above formula. However, 
the binding enerygy is about 10"^ of the mass of a proton and the mass fraction 
of lithium is only 10~^. Therefore, the effects on recombination should be very 
small. 

When we calculate the ionization history of cosmological recombination, we 
mainly talk about the hydrogen and helium because these two elements com- 
prise more than 99% of the total number of atoms in the Universe, particularly 
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in the primodial abundance. However, from the standard BBN, there are also 
tmy amount of deuterium (D) and hthium (Li) produced. Since Li^+ and Li"*" 
have higher ionization energies (122.4 and 75.6 eV respectively), they actually 
recombined before helium pO{ . On the other hand, neutral Li recombined at 
a much later time (z < 300) than hydrogen recombination ]8l[| . However, the 
lithium recombination brings negligible effect on Xg, because of its low abun- 
dance. Deuterium recombined at the same time as the rest of the hydrogen, due 
to having almost the same atomic structure, but with a heavier nucleus. Similar 
to Li, the abundace of D is also low (~ 10~^) and therefore, the recombination 
of D brings a negligible effect on ■ 

There is a tiny fraction of free electrons left {ric/nii ~ 10~^) after hydrogen 
recombination, and this allows for the formation of molecules in the later stages 
of evolution (see [|6[ |o[ and reference therein). Due to the high photon 
to baryon ratio, then at early times there are huge numbers of photons about 
the dissociation energy of H2 and hence coUisional processes (e.g. three-body 
reactions) are inefficient in molecule formation. The molecules are only pro- 
duced through radiative association [Q. For example, H2 (H~ -|-H^H2 +e~) 
is produced via the formation of H^ (H -t- e^ — > H~ + 7; see |g5|, |5^ for the latest 
calculations). Since the process of radiative association requires the existence of 
H I, significant production of molecules occurs only after recombination. These 
primordial molecules are important coolants in the star formation process and 
hence are crucial for understanding the formation of the first stars and galaxies, 
but they again have negligible effect on Xc- Due to the very low fraction of free 
electrons available for molecule formation, the abundance of these molecules is 
very low and they are also produced too late (2 < 300) to significantly affect 
the CMB photons. 

In all of this discussion we have focussed on the standard picture of recom- 
bination. Of course it is possible that we are still missing important pieces 
of the big picture. Some other non-standard physics could also easily alter 
the ionization fraction Xe at more than the percent level. Examples include a 
non-negligible interacting cross-section of dark matter ^1 , strong primor- 
dial magnetic fields ||2^, ^ , strong s pat ial inhomogeneities ||37| , 59 , extra Ly a 
emission from primordial black holes |64[| and a time- varying fine structure 



con- 
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2.4 Discussions 

In this chapter, we have briefly reviewed the recent updates and remaining un- 
certainties in the numerical recombination calculation. In order to obtain the 
ionization fraction to better than the percent level, then complicated details 
of the non-equilibrium situation need to be included. Most of the significant 
improvements have been mainly from additional radiative processes controlling 
the population of the n = 2 states. This is because there is no direct recombina- 
tion to the ground state and cascading down through n — 2 states is the main 
path for electrons to reach the ground state. If the existing studies have already 
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considered all the relevant physical processes in Hei recombination, then we 

currently have the corresponding numerical calculation to an accuracy better 
than 1%. However, for Hi recombination, there is still no single computational 
code which includes all of the suggested improvements. Hydrogen recombina- 
tion is even more important for calculating the CMB anisotropics Ce, because 
it dominates the detailed profile of the visibility function. A comprehensive 
numerical calculation of H i recombination, including at least all the suggested 
processes here, is neccessary and urgent in order to obtain high accuracy for 
future experiments. 
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Effect 



References 



Energy level 

Separate ^-states in H i atom 

Bound-bound transitions 

Inclusion of He I 23Pi-liSo 

Inclusion of He I n^Pi-liSo (n > 3) 
Inclusion of Hi ns, nd-ls (n > 3): 

- effective rate only 

- with feedback 

- with feedback and Raman scattering 
Inclusion of He I n^S, niD-fiSo {n > 3): 

- effective rate only 

- with feedback and Raman scattering 
Bound-free transitions 

Direct recombination for H i 
Direct recombination for He i 
Radiative transfer 
Continuum opacity of H I 

in He I 2iPi-liSo 
Feedback between Hei 23p-liSo 

and 2iPi-liSo 

Stimulated and induced H i 2s-ls 
Diffusion of Lya line profile 

(with recoil of H atoms) 
Continuum opacity of H i modified 

to feedback in He i lines 
Continuum opacity of H i in He i 

niP-fiSo, n^P-fiSo {n > 3), niD-fiSo 
Coherent scattering in n^P-l^So 
Evolution of Tm 
Secondary distortions from He i 

& H I in H I recombination 
Other 

He I 2'^Pi-l^So spontaneous rate 
CMB monopole uncertainty Tcmb ±lmK 
Primordial He abundance Yp ±1% 
Formation of hydrogen molecules 



-0.7% 
+1% 

-1.1% 
-0.3%* 
-0.004%* 

-0.4% 
-1.2% 
+ 1.3% 

-0.5% 
-0.05% 

-0.0006% 
-0.02% 

-2.5%* 

+1.5%* 

+0.6% 
~ -1% 

-0.5%* 

-0.05%* 

-0.02%* 
±0.001% 
+0.1% 



±0.1% 
±0.5% 

±1% 
-1% 



1090 
< 900 

1750 
1900 
2000 

1200 
1250 
900 

1800 
2000 

1280 
1900 



1800 
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900 
900 
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1900 
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45 


1 
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Table 2.1: Summary of the improvements and uncertainties in the numerical 
recombination calculation. Here Axe/xc is the maximum ratio difference of 
the ionization fraction Xe from the value given by recfast Version 1.3 and 
ZmaK is the approximate redshift at which this occurs. *Note: This is the relative 
change compared with the full radiative model in Switzer & Hirata (2008) |p3|. 
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Chapter 3 

Spectral Distortions^ 



3.1 Introduction 

Physical processes in the plasma of the hot early Universe thermalize the ra- 
diation content, and this redshifts to become the observed Cosmic Microwave 
Background (CMB; see and references therein). Besides the photons from 
the radiation background, there were some extra photons produced from the 
transitions when the electrons cascaded down to the ground state after they 
recombined with the ionized atoms. The transition from a plasma to mainly 
neutral gas occurred because as the Universe expanded the background tem- 
perature dropped, allowing the ions to hold onto their electrons. The photons 
created in this process give a distortion to the nearly perfect blackbody CMB 
spectrum. Since recombination happens at redshift z~1000, then Lya is ob- 
served at ~ 100/im today. There are ^ 1 of these photons per baryon, which 
should be compared with the ~ 10^ photons per baryon in the entire CMB. 
However, the recombination photons arc superimposed on the Wien part of the 
CMB spectrum, and so make a potentially measurable distortion. 

From the Far-Infrared Absolute Spectrophotometer (FIRAS) measurements, 
Fixsen et al. (1996) ^ and Mather et al. (1999) [|4| showed that the CMB is 
well modelled by a 2.725 ±0.001 K blackbody, and that any deviations from this 
spectrum around the peak are less than 50 parts per million of the peak bright- 
ness. Constraints on smooth functions, such as fi- or y-distortions are similarly 
very stringent. However, there are much weaker constraints on narrower features 
in the CMB spectrum. Moreover, within the last decade it has been discovered 
that there is a Cosmic Infrared Background (CIB; see and references 
therein), which peaks at 100-200/im and is mainly comprised of luminous in- 
frared galaxies at moderate redshifts. The existence of this background makes 
it more challenging to measure the recombination distortions than would have 
been the case if one imagined them only as being distortions to Wien tail of the 
CMB. However, as we shall see, the shape of the recombination line distortion is 
expected to be much narrower than that of the CIB, and hence the signal may 
be detectable in a future experiment designed to measure the CIB spectrum in 
detail. 

The first published calculations of the line distortions occur in the semi- 



version of this chapter (except Section 3.6) has been published: Wong W. Y., Seager S. 
and Scott D. (2006) 'Spectral distortions to the cosmic microwave background from the re- 
combination of hydrogen and helium', Monthly Notices of the Royal Astronomical Society, 
367, 1666-1676. 
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nal papers on the cosmological recombination process by Peebles (1968) p9[ 
and Zel'dovich et al. (1968) j5^. One of the main motivations for studying the 
recombination process was to answer the question: 'Where are the Ly a line 
photons from the recombination in the Universe?' (as reported in |Q). In fact 
these studies found that for hydrogen recombination (in a cosmology which is 
somewhat different than the model favoured today) there are more photons cre- 
ated through the two-photon 2s-ls transition than from the Ly a transition. 
Both Peebles (1968) @ and Zel'dovich et al. (1968) ||§ plot the distortion 
of the CMB tail caused by these line photons, but give no detail about the 
line shapes. Other authors have included some calculation or discussion of the 
line distortions as part of other recombination related studies, e.g. Boschan 
& Biltzinger (1998) ||l|, and most recently Switzer & Hirata (2005) H]. How- 
ever, the explicit line shapes have never before been presented, and the helium 
lines have also been neglected so far. The only numerical study to show the 
hydrogen lines in any detail is a short conference report by Dell'Antonio & Ry- 
bicki (1993) meant as a preliminary version of a more full study which never 
appeared. Although their calculation appears to have been substantially cor- 
rect, unfortunately in the one plot they show of the distortions (their figure 2) it 
is difficult to tell precisely which effects are real and which might be numerical. 

Some of the recombination line distortions from higher energy levels, n > 2, 
have also been calculated |, |, |, |l|, |l|, |l|, |l|, |2|, |3|. However, these high n 
lines lie near the peak of the CMB and therefore are extremely weak compared 
with the CMB (below the 10"'' level), while the Lya hue is well above the CMB 
in the Wien region of the spectrum. 

As trumpeted by many authors, we are now entering into the era of precision 
cosmology. Hence one might imagine that future delicate experiments may be 
able to measure these line distortions. Since the lines are formed by the photons 
emitted in each transitions of the electrons, they are strongly dependent on the 
rate of recombination of the atoms. The distortion lines may thus be a more 
sensitive probe of recombination era physics than the ionization fraction Xc, 
and the related visibility function which affects the CMB anisotropics. This is 
because a lot of energy must be injected in order for any physical process to 
change Xe substantially (for example, |^o|| ). In general that energy will go into 
spectral distortions, including boosting the recombination lines. 

This also means that a detailed understanding of the physics of recom- 
bination is crucial for calculating the distortion. The basic physical picture 
for cosmological recombination has not changed since the early work of Pee- 
bles (1968) and Zel'dovich et al. (1968) l^^l- However, there have been sev- 
eral refinements introduced since then, motivated by the increased emphasis on 
obtaining an accurate recombination history as part of the calculation of CMB 
anisotropics. Seager et al. (1999,2000) ||5^, |l) presented a detailed calculation 
of the whole recombination process, with no assumption of equilibrium among 
the energy levels. This multi-level computation involves 300 levels for both 
hydrogen and helium, and gives us the currently most accurate picture of the 
recombination history. In the context of the Seager et al. (2000) [|l| recombi- 
nation calculation, and with the well-developed set of cosmological parameters 
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provided by Wilkinson Microwave Anisotropy Probe {WMAP; [^) and other 
CMB experiments, it seems an appropriate time to calculate the distortion lines 
to higher accuracy in order to investigate whether they could be detected and 
whether their detection might be cosmologically useful. 

In this Chapter we calculate the line distortions on the CMB from the 2p-ls 
and 2s-ls transitions of Hi and the corresponding lines of He (i.e. the 2^p-l^s 
and 2^s-l^s transitions of He I, and the 2p-ls and 2s-ls transitions of Hell) 
during recombination, using the standard cosmological parameters and recom- 
bination history. In Section 3.2 we will describe the model we used in the nu- 
merical cal cula tion and give the equations used to calculate the spectral lines. 
In Section 3.3 we will present our results and discuss the detailed physics of 
the locations and shapes of the spectral lines. An approximate formula for the 
magnitude of the distortion in different cosmologies will also be given. Other 
possible modifications of the spectral lines and their potential detectability will 
be discussed in Section 3.4. And finally, we present our conclusions in the last 
section. 



3.2 Basic theory 
3.2.1 Model 

Instead of adopting a full multi-level code, we use a simple 3-level model atom 
here. For single-electron atoms (i.e. Hi and Heii), we consider only the ground 
state, the first excited state and the continuum. For the 2-electron atom (Hei), 
we consider the corresponding levels among singlet states. In general, the upper 
level states are considered to be in thermal equilibrium with the first excited 
state. Case B recombination is adopted here, which means that we ignore 
recombinations and photo-ionizations directly to ground state. This is because 
the photons emitted from direct recombinations to the ground state will almost 
immediately reionize a nearby neutral H atom |3^, ^ . We also include the 
two-photon rate from 2s to the ground state for all three atoms, with rates: 
A^s-is = 8.229063 s-i 0; ^fi^^ii^ = 51.02s-i§, although it makes no 
noticeable difference to the calculation if one uses the older value of 51.3 s~^ 
from Drake, Victor & Dalgarno (1969) §; and A|<=i\3 = 526.532 s"^ |l|, |l|. 
This 3-level atom model is similar to the one used in the program recfast, 
with the main difference being that here we do not assume that the rate of 
change of the first excited state n2 is zero. 

The rate equations for the 3 atoms are similar, and so we will just state the 
hydrogen case as an example: 

+ z)^^ = -^[Ai?H + Ai?|_ J + 3nf; (3.1) 
dn^(z) 1 

~ARg^,,]+3nf; (3.2) 



Chapter 3. Spectral distortion,s 



45 



, . dne(z) 1 r H /I To 



(1 + ^) 



dnp(z) 
dz 



1 



Hiz) 



(3.3) 
(3.4) 



Here the values of Ui are the number density of the zth state, where rie and rip 
are the number density of electrons and protons respectively. Ai?^^- is the net 
bound-bound rate between state i and j and the detailed form of Ai?^p_j^g and 
Ai?^g_]^g will be discussed in the next subsection. H{z) is the Hubble factor, 



H{zr = Hi 



1 + Zcq 



-(1 + Z)'* + an(l + zf + f^A'(l + Z)^ + 



(3.5) 



Here Zoq is the redshift of matter-radiation equality [ pT| , 

1/2 2 

where C/ is radiation energy density U = ctRlR, an is the radiation constant, fi, is 
the neutrino contribution to the energy density in relativistic species. Finally 
an is the Case B recombination coefficient from Hummer (1994) [ p3[ , 

at'' 

an = 10"^^- ^ m^s-\ (3.7) 

which is fitted by Pequignot et al. (1991) with a = 4.309, b = -0.6166, 
c = 0.6703, d = 0.5300 and t = Tm/IO^K, while /3h is the photo-ionization 
coefficient: 

3 

( 2TTmckBTM \ " _ J hpV2s.,c 



*--(^^J -n-iSfl- 

where Tm is the matter temperature and 1^25,0 is the frequency of the energy 
difference between state 2s and the continuum. For the rate of change of Tm, 
we only include the Compton and adiabatic cooling terms pT| , i.e. 

+ = :777n 7 1 (7m -Tr) -|-2rM, (3.9) 



where c is the speed of light and ctt is the Thompson scattering cross-section. 

We use the Bader-Deuflhard semi-implicit numerical integration scheme (see 
Section 16.6 in [Q) to solve the above rate equations. All the numerical results 
in this chapter are made using the ACDM model with parameters: flh = 0.046; 
n,^ = 0.3; f^A = 0.7; f^K ^ 0; Yp = 0.24; Tq = 2.725 K and h = 0.7 (see for 
examples, |5^). 
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3.2.2 Spectral distortions 

We want to calculate the specific line intensity 1^^ (z — 0) (i.e. energy per 
unit time per unit area per unit frequency per unit solid angle, measured in 
Wm~^Hz~^sr~^) observed at the present epoch, z = 0. The detailed calculation 
of Ii'aiz — 0) for the Lya transition and the two-photon transition in hydrogen 
are presented as examples (the notation follows Section 2.5 in |^). A similar 
derivation holds for the corresponding transitions in helium. To perform this 
calculation we first consider the emissivity ju{z) (energy per unit time per unit 
volume per unit frequency, measured in Wm~'^Hz~^) of photons due to the 
transition of electrons between the 2p and Is states at redshift z: 

3,{z) = h^y^R}^^_^^{z)cp[y{z% (3.10) 

where '^{v) is the frequency distribution of the emitted photons from the emis- 
sion process and AB^^_i^ is the net rate of photon production between the 2p 
and Is levels, i.e. 

Ai??p_is = Pi2 «i?2i ~ n'^Ri2) ■ (3.11) 

Here n^^ is the number density of hydrogen atoms having electrons in state i, 
the upward and downward transition rates are 

Ri2 = B,2J, (3.12) 
and R21 = {A21+B21J), (3.13) 

with A21, B12 and B21 being the Einstein coefficients and pi2 the Sobolev escape 
probability (see ||5l|), which accounts for the redshifting of the Lya photons due 
to the expansion of the Universe. As nY ^ n^p, P12 can be expressed in the 
following form: 

1 - e^'^' 

Pi2^ ,with (3.14) 



^2iAip43 (g2p/ffi)ni 
87riJ(z) 



(3.15) 



We approximate the background radiation field J as a perfect blackbody spec- 
trum by ignoring the line profile of the emission (see ||5l| ) . We also neglect sec- 



ondary distortions to the radiation field (but see the discussion in Section 3.4.1). 
These secondary distortions come from photons emitted earlier in time, during 
recombination of H or He, primarily the line transitions described in this paper. 
Assuming a blackbody we have 



(3.16) 



where Va = c/121.5682 nm= 2.466 x 10^^ Hz and corresponds to the energy 
difference between states 2p and Is, while the frequency of the emitted photons 
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is equal to v^- Therefore, we can set — 6[i'{z) — i.e. a delta function 

centred on i/a, so that 

j^y^iz) = hp:.ARl_,,{z)5Mz) - v^]. (3.17) 

The increment to the intensity coming from time interval dt at redshift z is 

dh{z) = ^j^dt, (3.18) 
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which redshifts to give 



d/.„(z = 0) = ^^^dt. (3.19) 



We assume that the emitted photons propagate freely until the present time. 
Integration over frequency then gives 

/^(-O)^^/^.. (3.20) 



_ chp ARfp_j^^{za) 
An if(z„)(l + z„)3^ 

with 

^ ~r Zq^ — , 

using 

1 

i/(z) = 1^0 (1 + and 



(3.21) 



dz H{z){l + z) 



Equation ( 3.21 ) is the basic equation for determining the Lya line distortion, 
using AR^p_iJz) from the 3-level atom calculation. 

For the two-photon emission between the 2s and Is levels, the emissivity at 
each redshift is 

Mz) = hpiyARl_,,iz)^[iyiz)i (3.22) 

and the calculation is slightly more complicated, since for 4>{h') we need the 
frequency spectrum of the emission photons of the 2s-ls transition of H ^ 
as shown in Fig. 3.1. Here A_R^g_]^g is the net rate of photon production for the 
2s-ls transition, i.e. 

^Rfs-is = A?s-is {nl - nfe-^-^-/'^-^^) . (3.23) 
Therefore, using equation ( 3.20 ), we have 

r27/ c^P'^o r Ai?|_i,(z)0K(l + z)] ^ 

We u se th e simple trapezoidal rule (see Section 4.1 in |42|) to integrate equa- 
tion ( [3.24 ) numerically from z = to the time when AR is sufficiently small 
that the integrand can be neglected. 
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Figure 3.1: The normalized emission spectrum for the two-photon process (2s- 
Is) of hydrogen The top panel shows 4>{X) vs A, while the bottom panel 

shows (p{y) vs y, where ly = yi^a- Note that the spectrum is symmetric in 
about 1^0/2, but the A spectrum is very asymmetric, being zero below Ac, and 
having a tail extending to high A. 
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3.3 Results 



Each of the hne distortions is shown separately in Fig. 3.2 and summed for each 



species in Fig. 3.3. The shape of the hnes from Hi, Hei and Hen are fairly 
similar. There are two distinct peaks to the 2p-ls emission lines. We refer to 
the one located at longer wavelength as the 'pre-recombination peak', since the 
corresponding atoms had hardly started to recombine during that time. The 



physics of the formation of this peak will be discussed in detail in section 3.3.1 
The second (shorter wavelength) peak is the main recombination peak, which 
was formed when the atoms recombined. While the longer wavelength peak 
actually contains almost an order of magnitude more flux, it makes a much 
lower relative distortion to the CMB. The ratio of the total distortion to the 



CMB intensity is shown in Fig. |3.4 It is ~ 1 for the main recombination peak, 
but - 10"'' for the pre-recombination peak. 



In Fig. 3.3, we plot the lines from H i and He i together with the CMB and an 
estimate of the CIB. We can see that the lines which make the most significant 
distortion to the CMB are the Lya line and the 2^p-l-'^s line of Hei, and that 
these lines form a non-trivial shape for the overall distortion. The sum of all the 



spectral lines and the CMB is shown in Fig. 3.2. Note that these lines will also 
exist in the presence of the CIB - but the shape of this background is currently 
quite poorly determined ^ . 

We now discuss details of the physics behind the shapes of each of the main 
recombination lines. 

3.3.1 Lines from the recombination of hydrogen 

During recombination, the Lyman lines are optically thick, which means that 
nearly all photons emitted from the transition to ?i = 1 are instantly reabsorbed. 
However, some of the emitted photons redshift out of the line due to the expan- 
sion of the Universe and this makes the Ly a transition one of the possible ways 
for electrons to cascade down to the ground state. The other path for electrons 
going from n = 2 to n = 1 is the two-photon transition between 2s and Is. 



Fig. 3.5 shows the net photon emission rate of the Lya and two-photon transi- 
tions as a function of redshift for the standard ACDM model. The two-photon 
rate dominates at low redshift, where the bulk of the recombinations occur. 
This means that there are more photons emitted through the two-photon emis- 
sion process (54% of the total number of photons created during recombination 
of H) than through the Ly a rcdshifting process. This conclusion agrees with 
Zeldovich et al. (1968) - although of course the balance depends on the 
cosmological parameters (see Q) and for today's best fit cosmology the two 
processes are almost equal. Despite this fact, the overall strength of the two- 
photon emission lines are weaker because the photons are not produced with a 
single frequency, but with a wide spectrum ranging from to fo,. The location 
of the two-photon peak (see Fig. ^) is also somewhat unexpected, since it is 
almost at the same wavelength as the Ly a recombination peak, rather than 
at twice the wavelength. The reason for this will be discussed in the following 
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100 1000 
(nm) 



Figure 3.2: The line intensity Aq/aq from the net Ly a emission of H (thick sohd), 
the two-photon emission (2s-ls) of Hi with the spectrum ^(i/) (thin sohd), the 
2^p-l^s emission of He I (thick dashed), the 2^s-l^s two-photon emission of Hei 
(thin dashed), the 2p-ls emission of Hell (thick dotted) and the 2s-ls two- 
photon emission of Hen (thin dotted). 
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Figure 3.3: The line intensity Aq/aq from the sum of the net Ly a emission and 
two-photon emission (ls-2s) of H (thick sohd), the sum of the 2ip-l^s emission 
and 2is-lis two-photon emission of Hel (thick dashed), and the sum of the 2p- 
Is emission and 2s-ls two-photon emission of Hell (thick dotted) , together with 
the background spectra: CMB (long-dashed); and estimated CIB (dot-dashed; 
[0) The sum of all the above emission lines of H and He plus the CMB is also 
shown (thin solid). 
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subsection. 

We should also note that the tiny dip in our curves for the long-wavelength 



tail of the pre-recombination peak (see Fig. 3.2) is due to a numerical error, 
when the number density of the ground state is very small. This can also be 
seen in the pre-recombination peak for Heii. 

The pre-recombination emission peak 

The highest Lya peak (shown in Fig. is formed before the recombination 
of H has already started, approximately at z > 2000. During that time the 
emission of Ly a photons is controlled by the bound-bound Ly a rate from n = 2 



(i.e. the ?t.2^21 term in equation (3.11)) and the photo-ionization rate (?T.2aH)- 
From Fig. we can see that at early times the bound-bound Lya rate is 
larger than the photo-ionization rate. This indicates that when an electron 
recombines to the n = 2 state, it is more likely to go down to the ground state 
by emission of a Ly a photon than to get ionized. The excess Ly alpha photons 
are not reabsorbed by ground state H, but are redshifted out of the absorption 
frequency due to the expansion of the Universe; they escape freely and form the 
pre-recombination emission line. Note that there is very little net recombination 
of H, since the huge reservoir of > 13.6 eV CMB photons keeps photo-ionizing 



the ground state H atoms (see Fig. p.l2[) 



We now turn to a more detailed explanation of the pre-recombination emis- 
sion peak. The bound-bound Ly a rate from n = 2 is initially approximately 
constant, as it is domi nated by the spontaneous de-excitation rate (the A21 



term in equation (3.23)). At the same time the photo-ionization rate is al- 
ways decreasing as redshift decreases, since the number of high energy photons 
keeps decreasing with the expansion of the Universe. Therefore, with a constant 
bound-bound Lya rate and the decreasing photo-ionization rate, the emission 
of Ly a photons rises. The peak of this pre-recombination line of H occurs at 
around z = 3000, by which time only a very tiny amount of ground state H atoms 
have formed {ni/nu < 10^'', see Fig. 3.6). These ground state H atoms build up 



until they can reabsorb the Ly a photons and this lowers the bound-bound Ly a 
rate. The decrease of the bound-bound Lya rate is represented in the Sobolev 



escape probability pi2 in equation (3.14). At high redshift, pi2 is 1 and there is 



no trapping of Lya photons. When H starts to recombine, the optical depth Tg 
increases and the Ly a photons can be reabsorbed by even very small amounts 
of neutral H. For Tg S> 1, we can approximate pi2 ~ l/rg and pi2 oc H{z)/ni. 
Because of the increase of the number density of the ground state and the de- 
crease of H{z), the pre-recombination line decreases. One can therefore think 
of the 'pre-recombination peak' as arising from direct Lya transitions, before 
enough neutral H has built up to make the Universe optically thick for Lyman 
photons. This process occurs because the spontaneous emission rate {A21 term) 
is faster than the photo-ionization rate for n = 2; it increases as the Universe 
expands, due to the weakening CMB blackbody radiation, and is quenched as 
the fraction of atoms in the n = l level grows. The shorter wavelength peak, 
on the other hand, comes from the process of redshifting out of the Ly a line 
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Figure 3.4: The ratio of the total Une distortion to the CMB intensity is plotted. 

The ratio is larger than 1 (i.e. the intensity of the distortion line is larger than 
that of the CMB) when Aq ~ 170 /um which is just where the main Lya line 
peaks. 
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Figure 3.5: Comparison of the net 2p-ls (solid) and 2s-ls (dashed) transition 
rates of H. The Ly a rcdshifting process dominates during the start of recom- 
bination, while the 2-photon process is higher during most of the time that 
recombination is occurring. It turns out that in the standard ACDM model 
about equal numbers of hydrogen atoms rccombinc through each process, with 
slightly over half the hydrogen in the Universe recombining through the 2-photon 
process. 




Figure 3.6: The top panel shows the bound-bound Lya rate n2-R2i and the 
photo-ionizing rate n2aH for n=2. The lower panel shows the fraction of ground 
state H atoms ni/nu, and also the ionization fraction Xe- 
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during the bulk of the recombination epoch. 

By using the recfast program , we can generate the main Ly a recom- 
bination peak and also the two-photon emission spectrum, by simply adding a 
few lines into the code. However, the pre-recombination peak cannot be gen- 
erated from RECFAST, since there the rate of change of the number density of 
the first excited state 712 is assumed to be negligible and is related to ni via 
thermal equilibrium. Moreover, in the effective 3-level formalism, the Ly a line 
is assumed to be optically thick throughout the whole recombination process of 
H (in order to reduce the calculation into a single ODE), which is not valid at 
the beginning of the recombination process. Hence, one needs to follow the rate 
equations of both states (i.e. n — 1 and n = 2) to generate the full Ly a emission 
spectrum. The pre-recombination peak of H was mentioned and plotted in the 
earlier work of Dell'Antonio & Rybicki (1993) Q as well, although they did not 
describe it in any detail. 

Another way to understand the line formation mechanism is to ask how many 
photons are made in each process per atom. We find that for the main Ly a 
peak there are approximately 0.47 photons per hydrogen atom (in the standard 
cosmology). During the recombination epoch, net photons for the n = 2 to 
n — I transitions are only made when atoms terminate at the ground state. 
Hence we expect exactly one n = 2 to n = 1 photon for each atom, split between 
the Lya redshifting and 2-photon processes (and the latter splits the energy 
into two photons, so there are 1.06 of these photons per atom). For the 'pre- 
recombination peak', on the other hand, the atoms give a Lya photon when 
they reach n = 1, but they then absorb a CMB continuum photon to get back 
to higher n or become ionized. The number of times an atom cycles through 
this process depends on the ratio of the relevant rates. If we take the rate 



per unit volume from Fig. 3.5 and divide by the number density of hydrogen 
atoms at z ~ 3000 then we get a rate which is about an order of magnitude 
larger than the Hubble parameter at that time. Hence we expect about 10 'pre- 
recombination peak' photons per hydrogen atom. A numerical calculation gives 
the more precise value of 8.11. 

The tvifo-photon emission lines 

Surprisingly, the location of the peak of the line intensity of the 2s-ls transition 



is almost the same as that of the Lya transition, as shown in Fig. 3.2, while 



one might have expected it to differ by a factor of 2. In order to understand 



this effect, we rewrite the equation (3.24) in the following way: 



/2;'(z = 0) = / <^'(z')< = 0; z'] dz\ (3.25) 
where 4>'{z') = and 



(3.26) 
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with 

1 + z' = —. 

Equation (3.26) gives the redshifted flux (measured now at z = 0) of a single fre- 
quency v' coming fr'om redshift z' and corresponding to the redshifted frequency 

We first calculate the line intensity of the two-photon emission with a sim- 
ple approximation: a delta function spectrum 5{v — i/q,/2), where 1^0/2 is the 
frequency corresponding to the peak of the two-photon emission spectrum 4>{i^)- 
Fig. 3.7 shows the intensity spectrum of two-photon emission using a delta fre- 
quency spectrum 5{i' — 1^0/2) compared with the two-photon emission using the 
correct spectrum (f){v) . We can see that there is a significant shift in the line 
centre compared with the J-function case. Where does this shift come from? 

We know that the frequencies of emitted photons are within the range of 
to Va at the time of emission. For a fixed redshifted frequency vq now, we can 
calculate the range of emission redshifts contributing to vq (referred to as the 
'contribution period' from now on), which is represented by (f>'{z') or (j){v') . 
In Fig. 3.8, we show the spectral distribution 0[:/'(z')] as a function of redshift 
z' for specific values of vq. For example, if we take vq — 5 x 10^^ Hz, then 
photons emitted between 1 -I- 2; = 1 (i.e. v = vq) and ~500 {v = Va) will 
give contributions to vq. The smaller the redshifted frequency vq, the wider 
the contribution period. We might expect that the line intensity of this two- 
photon emission will be larger if the contribution period is longer, as there 
are more redshifted photons propagating from earlier times. However, this is 
not the case, because the rate of two-photon emission i?27 also varies with 
time, and is sharply peaked at z ~ 1300-1400. Hence I^^^ [z = 0; z'] is also 
sharply peaked at z ~ 1300-1400. In Fig. |3.8| , the redshifted fiux integrand 
I^g{z = 0, z) and the emission spectrum (/)[:/(z)] are plotted on the same redshift 
scale. For j/q = 5 x 10^^ Hz (lowest panel), we can see that the contribution 
period covers a redshift range when /^^ (z = 0, z) and i?27 are small in value. The 
contribution period widens with decreasing i^q and covers more of the redshift 
range when two-photon emission was high. Therefore, the flux 1^^ (z = 0) is 
expected to increase with decreasing vq until the contribution period extends to 
the redshifts at which the two-photon emission peaks. As fo gets even smaller 
(e.g. fQ = lO^^Hz), then the contribution period becomes larger than the redshift 
range for two-photon emission and hence only lower energy photons can be 
redshifted to that redshifted frequency. As a result, the fiux I,^„{z — 0) starts 
to decrease, and so we have a peak. The flux peaks at i^q ~ 10^^ Hz when we 
use the (S-function approximation. However, from Fig. 3.5, we can see that the 
contribution period for vq ~ 10^^ Hz is much greater than that of the two-photon 
emission period, and therefore this is not the location of peak. Based on the 
argument presented above, we expect the peak to be at around 1.6 x lO^'^Hz, 
or 200 fim. 

The basic mathematical point is that (j){y) is extremely poorly represented 
by a (5-function. Since the spectrum (f>{h') is quite broad, it can be better ap- 
proximated as a uniform distribution than as a (5-function. Another crude ap- 
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Figure 3.7: The line intensity of the 2s-ls transition (two-photon emission) 
Ivaiz = 0) as a function of redshifted frequency vq for three different assump- 
tions: the correct frequency spectrum of two-photon emission (solid); the delta 
function approximation 5{i' — Va/^) (dashed); and the flat spectrum approxi- 
mation (dotted). 
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Figure 3.8: The top panel shows the redshifted flux from single emission fre- 
quency lt,„{z = 0;z') plotted against the redshift of emission, 1 + z'. The 
bottom panel shows the frequency spectrum of two-photon emission 
plotted against z' for three redshifted frequencies: fo = 10^^ Hz; 1.6 x 10^^ Hz; 
and 5 x lO^^ jjz. 
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proximation would be to assume a flat spectrum for in Fig. |3.l| . Fig. \i.7\ 
compares the intensity 1^,^ [z — 0) found using the correct form for ipi^) with 
the (5- function and flat spectrum approximations. This shows that the flat spec- 
trum gives qualitatively the same results as the correct form of the spectrum, 
and that the peak occurs fairly close to that of Ly a, but is much broader. The 
same general arguments apply to the two-photon lines of He i and He ii (as we 
discuss in Section 3.3.2). 



Dependence of and fib 

The largest distortion on the CMB is from the shorter wavelength recombination 



peak of the hydrogen Ly a line (see Fig. 3^ ) . It may therefore be useful estimate 



the peak of this line's intensity as a function of the cosmological parameters. The 
relevant parameters are the matter density (oc rim^i^) and the baryon density 
(oc flhh^)- This is because ft^h^ affects the expansion rate, while il\^h^ is related 
to the number density of hydrogen. No other combinations of cosmological 
parameters have a significant impact on the physics of recombination. 

We can crudely understand the scalings of these parameters through the fol- 
lowing argument. Regardless of the escape probability pi2, the remaining part 
of the rate (n^pi?2i — "1^-^12) is roughly proportional to nf oc flhh'^{l — Xe)- 
The escape probability pi2 can be approximated as 1 at the beginning of recom- 
bination (ts <C 1) and 1/ts during the bulk of the recombination process (with 
Ts » 1). Note that oc H{z)/nf oc {n^h'^)'^/'^[nbh'^{l - Xe)]~^. Therefore, 

r {n,,h^nny,h^i - X,)] for n « 1 

oc I (f^^;,2)i/2[j^^^2(i _ ^^)]0 for » 1, ^^'^^^ 



and thus 



°' H{z) ^ 1 {n,^h^f[n^h\l - xe)]0 for T, » 1. ^''•^^^ 

From this rough scaling argument, we may expect that the dependence of 
the peak of the Ly a line is an approximate power law with index between —1/2 
and 0, while for f^b the corresponding power-law index is expected to lie between 
and 1. The dependence of is actually more complicated when one allows 
for a wider range of values (see Q). The above estimation just gives a rough 
physical idea of the power of the dependence. 

A more complete numerical estimate of the peak of the recombination Ly a 
distortion is: 

0.15 



ST-\ (3.29) 



where we have normalized to the parameters of the currently favoured cosmo- 
logical model. The peak occurs at 

Ao-lTO^im (3.30) 
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for all reasonable variants of the standard cosmology. 

3.3.2 Lines from the recombination of helium (He I and 
He II) 

We compute the recombination of He ii and He i in the same way as for hydrogen. 
For the two-electron atom He i, we ignore all the forbidden transitions between 
singlet and triplet states due to the low population of the triplet states (see 
[|o[ |5^). The 2^p-l^s transitions of He I are optically thick, the same situation 
as for H. This makes the electrons take longer to reach the ground state and 
causes the recombination of He i to be slower than Saha equilibrium. However, 
unlike for H, and despite the optically thick 2^p-l^s transition line, the 2^p-l-'^s 



rate dominates, as shown in Fig. 3.9. For Hen, due to the fast two-photon 
transition rate (see Fig. 3.1C), there is no 'bottleneck' at the n = 2 level in the 
recombination process. Hence Hen recombination can be well approximated by 
using the Saha equilibrium formula [|T] . 

We can see the effect of the above differences in recombination history on 
the lines: the width of the recombination peak of both H and He i is larger than 
that of Hen. Overall, the spectral lines of Hen are of much lower amplitude 



than those of H (see Fig. 3.2) with the distortion to the CMB about an order 
of magnitude smaller. 

The peaks of the line distortions from H and He n arc located at nearly the 
same wavelengths. For hydrogenic ions the ls-2p energy (and all the others) 
scales as Z^, where Z is the atomic number. Hence for Hen recombination 
takes place at z ~ 6000 rather than the z ~ 1500 for hydrogen. Hence the line 
distortion from the 2p-ls transition of He n redshifts down to about 200 /im, 
just like Lya. 

The two-photon frequency spectrum of He n is the same as for H, since they 
are both single-electron atoms However, it is complicated to calculate the 
two-photon frequency spectrum of He i very accurately, since there is no exact 
wave-function for the state of the atom. Drake et al. (1969) |^ used a variational 
method to calculate the two-photon frequency spectrum of He i with values given 
up to 3 significant figures. Drake (1986) |^ presented another calculation, giving 
one more digit of precision, and making the spectrum smoother, as shown in 



Fig. 3.11. These two calculations differ by only about 1 % , which makes negligible 
change to the two-photon He i spectral line. 

All of the H and He lines (for n = 2 to n=l) are presented in Fig. |3.2| and 



the sum is shown as a fractional distortion to the CMB spectrum in Fig. 3.4. We 
find that in the standard cosmological model, for He i recombination, there are 
about 0.67 photons created per helium atom in the 'main' 2^p-l^s peak, 0.70 
per helium atom in the 'pre-recombination peak', and 0.66 in the two-photon 
process. The numbers for Hen recombination are 0.62, 0.76 and 6.85 for these 
three processes, respectively. 
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Figure 3.9: Comparison of the net 2^p-l^s (solid) and 2^s-l^s two-photon 
(dashed) transition rates of Hei. The two-photon rate is sub-dominant through 
most of the Hei recombination epoch, and hence, unhke for hydrogen, most 
heUum atoms did not recombine through the two-photon process. 
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Figure 3.10: Comparison of the net 2p Is (solid) and 2s-ls two-photon (dashed) 
transition rates of Hell as a function of redshift. The two-photon process is 
greater through most of the recombination epoch, so that most of the cosmo- 
logical Helll — > Hell process happens through the two-photon transition. 
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Figure 3.11: The normalized emission spectrum for the two-photon emission pro- 
cess (2^s-l^s) in He I. Here y = i^/t^2s-is, where j^2s-1s = 4.9849 x 10^^ Hz. The 
crosses are the calculated points from Drake et al. (1969) and Drake (1986) 
while the line is a cubic spline fit. 
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3.4 Discussion 

3.4.1 Modifications in the recombination calculation 

There are several possible improvements that we could make to the line distor- 
tion calculation. However, as we will discuss below, we do not believe that any 
of them will make a substantial difference to the amplitudes of the lines. 

In order to calculate the distortion lines to higher accuracy, we should use 
the multi-level model without any thermal equilibrium assumption among the 
bound states. And we also need to take into account the secondary spectral 
distortion in the radiation field, i.e. we cannot approximate the background 
radiation field J as a perfect blackbody spectrum. This means, for example, that 
the extra photons from the recombination of He i may redshift into an energy 
range that can photo-ionize H(n=l) g We can assess how significant this 
effect might be by considering the ratio of the number of CMB background 
photons with energy larger than E^, n^{> E-y), to the number of baryons, tib, 



at different redshifts (see Fig. 3.12) 



Roughly speaking, the recombination of H occurs at the redshift when n^(> 
hpi'a)/nB is about equal to 1. This is because at lower redshifts there are not 
enough high energy background photons to photo-ionize or excite electrons from 
the ground state to the upper states (even to n = 2), while at higher redshift, 
when such transitions are possible, there are huge numbers of photons able to 
ionize the n = 2 level. The solid line in Fig. 3.12| shows the effect of the helium 



line distortions on the number of high energy photons (above Ly a) per baryon. 
The amount of extra distortion photons with redshifted energy larger than h^Va 
coming from the recombination of Hei is only about 1 per cent of the number 
of hydrogen atoms. Their effect is therefore expected to be negligibly small for 
Xe- We neglect the effect of the helium recombination photons on the hydrogen 
line distortion, since it is clearly going to make a small correction (at much less 
than the 10 per cent level). 

As well as this particular approximation, there have been some other recent 
studies which have suggested that it may be necessary to make minor mod- 
ifications to the recombination calculations presented in Seager et al. (1999, 
2000) |5^. Although these proposed modifications would give only small 
changes to the recombination calculation, it is possible that they could have 
much more significant effects on the line amplitudes and shapes. Recent papers 
have described 3 separate potential effects. 

In the effective three-level model, Leung et al. (2004) argued that the 
adiabatic index of the matter should change during the recombination process, 
as the ionized gas becomes neutral, giving slight differences in the recombination 
history. Dubrovich & Grachev (2005) |13| have claimed that the two-photon rate 
between the lowest triplet state and the ground state and that between the upper 
singlet states and the ground state should not be ignored in the recombination 
of Hei. And Chluba & Sunyaev (2005) Q suggested that one should also include 
stimulated emission from the 2s state of H, due to the low frequency photons in 
the CMB blackbody spectrum. Even if all of these effects are entirely completely 
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correct, wc find that the change to the ampUtude of the main spectral distortion 
is much less than 10%. We therefore leave the detailed discussion of these and 
other possible modifications to a future work. 



3.4.2 Possibility of detection 

There is no avoiding the fact that detecting these CMB spectral distortions will 
be difficult. There are three main challenges to overcome: (1) achieving the 
required raw sensitivity; (2) removing the Galactic foreground emission; and 
(3) distinguishing the signal from the GIB. 

Let us start with the first point. We can estimate the raw sensitivity achiev- 
able in existing or planned experiments (even although these instruments have 
not been designed for measuring the line distortion). Since the relevant wave- 
length range is essentially impossible to observe from the ground, it will be 
necessary to go into space, or at least to a balloon-based mission. One ex- 
isting experiment with sensitivity at relevant wavelengths is BLAST which 
has an array of bolometers operating at 250 /xm on a balloon payload. The 
estimated sensitivity is 236 mJy in 1 s, for a 30 arcsec FWHM beam, which 



corresponds to XI\ = 1.2 x lO^^Wm^'^sr"^. Comparing with equation (3.29) 
for the peak intensity, it would take ^ 10^ such detectors running for a year 
to detect the line distortion. The SPIRE instrument on Herschel will have a 
similar bolometer array, but with better beamsize. The estimated sensitivity 
of 2.5 mJy at 5cr in 1 hour for a 17.4 arcsec FWHM beamsize corresponds 
to A/a — 4.4 X 10~®Wm~^sr~^ per detector for the la sensitivity in 1 second. 
So detection of the line would still require ^ 10^ such detectors operating for a 
year. 

These experiments are limited by thermal emission from the instrument it- 
self, and so a significant advance would come from cooling the telescope. This 
is one of the main design goals of the proposed SAFIR and SPICA |3^ mis- 
sions. One can imagine improvements of a factor ~ 100 for far-IR observations 
with a cooled mirror. This would put us in the regime where arrays of - 10^ de- 
tectors (of a size currently being manufactured for sub-mm instruments) could 
achieve the desired sensitivity. 

One could imagine an experiment designed to have enough spectroscopic 
resolution to track the shape of the expected line distortion. The minimum 
requirement here is rather modest, with only X/6X 10 in at least 3 bands. 
An important issue will be calibration among the different wavelengths, so that 
the non-thermal shape can be confidently measured. To overcome this, one 
might consider the use of direct spectroscopic techniques rather than filtered or 
frequency-sensitive bolometers. 

Another way of quoting the required sensitivity is to say that any experi- 
ment which measures the recombination line distortion would have to measure 
the GIB spectrum with a precision of about 1 part in lO'^, which is obviously 
a significant improvement over what has been currently achieved. A detection 
of the line distortion might therefore naturally come out of an extremely pre- 
cise measurement of the GIB spectrum, which would also constrain other high 
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Figure 3.12: The ratio of number of CMB photons with energy larger than 
{nj{> Ej)) to number of baryons (ne) is plotted against redshift z. The solid 
line includes the extra distortion photons from the recombination of Hel. Prom 
the graph, we can sec that the recombination of H occurs approximately at the 
redshift when the ratio of photons with energy > hua to baryons is about unity. 
By the time the helium recombination photons are a significant distortion to 
the CMB tail above Ly a the density of the relevant photons has already fallen 
by 2 orders of magnitude, and so the effects can make only a small correction. 
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frequency distortions to the CMB spectrum. 

Some of the design issues involved in such an experiment are discussed by 
Fixsen & Mather (2002) jl^. They describe a future experiment for measur- 
ing deviations of the CMB spectrum from a perfect blackbody form, with an 
accuracy and precision of 1 part in 10^. This could provide upper limits on 
Bose-Einstein distortion n and Compton distortion y parameters at the ^ 10"'' 
level (the current upper limits for y and fi are 15 x 10~^ and 9 x 10"^, re- 
spectively; ||l^). The frequency coverage they discuss is 2-120 cm~^ (about 
80-5,000 yitm), which extends to much longer wavelengths than necessary for 
measuring the line distortion. The beam-size would be large, similar to FIRAS, 
but the sensitivity achieved could easily be 100 times better. An experiment 
meant for detecting the line distortion would have to be another couple of orders 
of magnitude more sensitive still. 

Turning to the second of the major challenges, it will be necessary to detect 
this line in the presence of the strong emission from our Galaxy. At 100 /xm 
the Galactic Plane can be as bright as ^ lO^MJysr"^ which is about a billion 
times brighter than the signal we are looking for! Of course the brightness falls 
dramatically as one moves away from the Plane, but the only way to confidently 
avoid the Galactic foreground is to measure it and remove it. So any experiment 
designed to detect the line distortion will need to cover some significant part of 
the sky, so that it will be possible to extrapolate to the cosmological background 
signal. The spectrum of the foreground emission is likely to be smoother than 
that of the line distortion, and it may be possible to use this fact to effectively 
remove it. However, it seems reasonable to imagine that the most efficient sep- 
aration of the signals will involve a mixture of spatial and spectral information, 
as is done for CMB data (see, for example, Q). 

In the language of spherical harmonics, the signal we are searching for is 
a monopole, with a dipole at the ^ 10"'^ level and smaller angular scale fluc- 
tuations of even lower amplitude. Hence we would expect to be extrapolating 
the Galactic foreground signals so that we can measure the overall DC level of 
the sky. This is made much more difficult by the presence of the GIB, which 
is also basically a monopole signal. Hence spatial information cannot be used 
to separate the line distortion from the GIB. The measurement of the line dis- 
tortion is therefore made much more difficult by the unfortunate fact that the 
GIB is several orders of magnitude brighter - this is the third of the challenges 
in measuring the recombination lines. 

The shape of the GIB spectrum is currently not very well characterised. 
It was detected using data from the DIRBE and FIRAS experiments on the 
COBE satellite. Estimates for the background {XI\) are: 9nWm~^sr~^ at 
60/xm||; 23nWm-2sr-i at 100Aim|7|; 15nWm-2sr-i at 140 ^im |l[ |§; 
and llnWm~^sr~^ at 240 /im [|l[ Hi; In each case the detections are only at 
the 3-5a level, and the precise values vary between different prescriptions for 
data analysis (see also |l5], ^ @)- The short wavelength distortion of the 
CMB, interpretted as a measurement of the GIB can be fit with a modified 
blackbody with temperature 18.5 K and cmissivity index 0.64 (although there 



is degeneracy between these parameters), which we plotted in Fig. 3.3 
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The CIB is thus believed to peak somewhere around 100 ^m, which is just 
about where we are expecting the recombination hne distortion. The accuracy 
with which the CIB spectrum is known will have to improve by about 5 orders 
of magnitude before the distortion will be detectable. Fortunately the spectral 
shape is expected to be significantly narrower than that of the CIB - the line 
widths are similar to the dz/z ^ 0.1 for the last scattering surface thickness, as 
opposed to (5A/A 1 for a modified blackbody shape (potentially even wider 
than this, given that the sources of the CIB come from a range of redshift 
Az ~ 1). 

One issue, however, is how smooth the CIB will be at the level of detail 
with which it will need to be probed. It may be that emission lines, absorption 
features, etc. could result in sufficiently narrow structure to obscure the recom- 
bination features. We are saved by 2 effects here: firstly the CIB averaged over a 
large solid angle patch is the sum of countless galaxies, and hence the individual 
spectral features will be smeared out; and secondly, the far-IR spectral energy 
distributions of known galaxies do not seem to contain strong features of the 
sort which might mimic the recombination distortion (see, for example, ||2^ ). 
As we learn more about the detailed far-IR spectra of individual galaxies we 
will have a better idea of whether this places a fundamental limit on our ability 
to detect the recombination lines. 

Overall it would appear that the line distortion should be detectable in 
principle, but will be quite challenging in practice. 

3.5 Conclusion 

We have studied the spectral distortion to the CMB due to the Lya and 2s- 
Is two-photon transition of Hi and the corresponding lines of Hei and Hen. 
Together these lines give a quite non-trivial shape to the overall distortion. 
The strength and shape of the line distortions are very sensitive to the details 
of the recombination processes in the atoms. Although the amplitude of the 
spectral line is much smaller than the Cosmic Infrared Background, the raw 
precision required is within the grasp of current technology, and one can imagine 
designing an experiment to measure the non-trivial line shape which we have 
calculated. The basic detection of the existence of this spectral distortion would 
provide incontrovertible proof that the Universe was once a hot plasma and its 
amplitude would give direct constraints on physics at the recombination epoch. 
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3.6 Remarks 

Since this work was published, there have been other studies calculating the 
same spectral distortions with different approach in a different independent 
numerical code Rubino-Martm et al. (2006) |^ pointed out a correction 

in the treatment of the two-photon spectrum, and found no pre-recombination 
peak in the H Ly a line distortion, in contradiction to the results of this chapter. 
As an addition to our published study, we now discuss these two issues. 



Normalization of the two-photon spectrum 

In our calculation of the two-photon line distortion, the emission spectrum 



is normalized to 1 (see Figure B.l and Equation (3.22)). However, the two- 
photon spectrum </)(i/) should be normahzed to 2 (as pointed out by |Q) because 
there are two photons emitted in each electron transition from the 2s state to 
the ground state. Due to this correction, the intensity of the two-photon line 
distortion presented before should be doubled. For H i, since the amplitude of 
the distortion from Lya emission is about 10 times larger than the two-photon 
contribution, the overall shape and the peak location of the line spectrum remain 
almost the same as before. The same correction should be made for the helium 
line distortion spectrum as well, and again the effect or the overall distortion 
from He is small. 



The pre-recombination peak 

Rubiho-Martfn et al. (2006) performed an independent calculation of the 
spectral line distortions from H i recombination with a multi-level atom model. 
The authors adopted the same procedure described in Seager et al. (2000) |^ 
but considered separate Z-states within each n-shell of H i with no thermal equi- 
librium assumption. In this Chapter, we obtained a pre-recombination peak 
using a 3-level atom model also based on the recombination model given by 
Seager et al. (2000) In contrast, Rubino-Martm et al. (2006) ||| found no 
pre-recombination peak in their calculation. 

As with earlier work , our pre-recombination peak was only found from 
numerical calculation and no explicit theoretical argument for the formation of 
this peak was given. We can consider the calculation a different way in order to 
understand the underlying physics. Since the population of the hydrogen atom 
states is well described by the Boltzmann equations before recombination (say 
z < 1700 for Hi; see for example [Q), we now present an analytical estimate of 
the H I pre-recombination peak under the local thermal equilibrium assumption 
in a 3-level atom model. 

The pre-recombination peak was previously found in the calculation of the 
Hi Lya emission line, and also for the corresponding He I and Heii emission 
lines. Here we only discuss the case of Hi, since the physics is basically the 
same for other species within the standard recombination model |Q . Since the 
spectrum of the photon emission in this transition is narrowly peaked at the Ly a 
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frequency, the distortion shape is mainly controUed by the net Ly a emission 
rate Ai?2p-is (see Equation ( |3.21 )). From Equation (3.11), the net Lya rate 
can be rewritten as 



AR 



•2p-ls 



pi2nisA2i 



n2p 



92p ^-hpu^/kBTn 



(3.31) 



Here we approximate Tm — Tr. We can also write the net 2s-ls two-photon rate 
as 

AR2S-U = A|_i,ni3^ - 52p^_,,.„/fc3TH\ (3 33) 

52p \nis gis ) 

by assuming that the 2p and 2s states are in thermal equilibrium. From the 
above equations, we can see that these two rates are controlled by the same 
difference, i.e. the difference between the ratio ri2p/nis and its local thermal 
equilibrium value from the Boltzmann factor. 

Rubiho-Martfn et al. (2006) argued that A_R2p-is is equal to zero at 
z > 2000 because the states are in thermal equilibrium. This is not entirely 
true, since the expanding Universe is a fundamentally out-of-equilibrium system; 
we will show that Ai?2p_is is non-zero (although the rate is very low) even 
if the population of the states in H i is well approximated by the Boltzmann 
distribution at each instant of time during the pre-recombination period. From 
the Saha equation, we have 



where Ui^c is the frequency of the energy difference between the zth state and 



the continuum. We can take Equation (3.33) for i = l (Is, the ground state), 



differentiate with respect to z and substitute into Equation (3.1), giving 

Ai?2^™, + AR^J\ ^ n,,Hiz) - \ - • (3.34) 

The right-hand side of the above equation is dominated by the first two terms, 
since 

/ipi/i^c 5.792 X 10^ 



with Tr = 2.725(1 + z) K, (3.35) 

z dxn 



^bTr 1 

and 

< 0.1 (3.36) 



before the recombination of H i {z~ 1800) . This makes the sum of the two rates 
larger than zero and implies that there are net recombinations to the ground 
state even in the case that the number density of each state closely follows 
the thermal equilibrium distribution. Physically, the non-zero net recombina- 
tion rate of H i is due to the decreasing number of high-energy photons in the 



Chapter 3. Spectral distortion,s 



72 



expanding Universe. And as we know from Equation (3.31) and (3.32), 

Ai?2p-ls 3^12^21 



Ai?2s-ls 



AH 



10'' at z > 1700 . 



(3.37) 



Before the Hi recombination (2:> 1800 say), the net Lya rate dominates, be- 
cause the neutral hydrogen abundance is very low and the escape probability 
Pi2 is very close to 1. Therefore, we can ignore Ai?2™is in Equation (3.34) and 
we have 



Ai?: 



2p-ls 



(3.38) 



In Figure 3.13, the approximate rate Ai?2J^'is is plotted along with the previ- 



ous result from the numerical recombination code. We can see that Ai^jJ-is 
matches the numerical rate Ai?2p~is very well when the hydrogen is about to 
recombine at z =1800-2000. These two rates are expected to depart at z ~ 1750 
when the ground state goes out of thermal equilibrium with the higher excited 
states due to the bottleneck at the first excited state. On the other hand, we 
expect the thermal equilibrium assumption to be valid at even higher redshifts 
(z>2000). Under this assumption, we find no significant emission before re- 
combination and therefore, there is no pre-recombination peak. Why does this 
result contradict with what we found in the numerical calculation? 

In fact we found that the pre-recombination peak that we presented before 
arose due to a systematic error in the ODE (ordinary differential equation) 
solver. Any ODE solver allows us to find a numerical approximation to an 
exact (or real) solution of the equations to within some error. We usually want 
the relative error to be small, and this is controlled by setting the required 
accuracy as an input parameter in the solver. In our case, the relative error of 



the number density is Arii/rii — (n" 
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Now consider the effect of this error on the net rates. The net rates are 
strongly controlled by the deviation of the ratio n2p/?iis from its Boltzmann 
value. Note that 
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where ei2 — Ari2p/n2p — Ariig/nis, should be the same order of magnitude as 
the uncertainty in Ui (i.e. Ani/rii). In the above equation, the first bracket 
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2000 4000 6000 

z 

Figure 3.13: The net Lya transition rate Ai?2p-is (solid line) and the net 2s-ls 
two-photon transition rate Ai?2s-is (dashed line) of H i as a function of redshift 
z. These two curves are generated from the multi-level numerical recombination 
code. The dotted line (red) is the approximate analytical Lya transition rate 
Ai?^™i, from Equation (|^. 
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accounts for how much the first excited state and the ground state are out of 
equihbrium and this gives us the actual net Lya rate. The second term is 
the error in the Ly a rate due to the numerical errors in the number densities. 
Somewhat surprisingly, it is directly proportional to the actual value of n2p/nis, 
which increases with z. In the pre-recombination epoch, we can approximate 
{n2p/nisy°'^^ using Boltzmann equations in order to calculate the error of the 
rate. For comparison, we use Equation (3.31) to obtain the estimate 



A 



/ \ real A pLTE 



V "Is / Pi2n 



ls^21 V / 



In Fig. 3.14, we separately plot the two terms in Equation (3.39) as well as 
^("■2p/?^is)""™ from the numerical code. The estimated numerical error dom- 
inates at z > 2000 and it matches well with the A(n2p/nis)""™ curve if we 
take ei2 = lO"'''®, which is even smaller than the required accuracy in the ODE 
solver. This error term explains why there is an anomalous increasing trend 
of A(n2p/rtis)""™ at high z, while we expect the difference in the ratio to get 
smaller with increasing z due to the tight thermal equilibrium relation between 
the states. This estimated numerical error is directly proportional to rt2p/"-is 
and decreases with decreasing z. On the other hand, A(n2p/7iis)'''^'*' is getting 
larger and larger as the recombination of hydrogen begins. So at z ~ 2000, 
A(n2p/nis)''°''' takes over. This explains why the Ai^^^^^ and AR^"^^^^ val- 
ues agree with each other only in the range of z ~ 1600-2000. Overall, the 
pre-recombination peak that we found earlier seems to have beem caused by a 
systematic error. This should serve as a warning for blindly accepting numerical 
result. 

Physically, any possible pre-recombination peak or extra emission through 
the H I Ly a transition would require channels for electrons in the ground state 
to get back to higher excited states or the continuum, since there is almost no 
net neutral hydrogen H i formed due to these processes. In the three-level atom 
model, no such path exists since all the transitions to and from the ground 
state are connected with the first excited state {n — 2). We have also investi- 
gated this problem in a multi-level atom model based on the paper of Seager et 
al. (2000) and still find no such net excitation to the higher excited states. 
From this physical reasoning and the previous analysis of numerical errors, the 
pre-recombination peak we found in the standard recombination calculation 
therefore appears to have been a false signal. 

On the other hand, additional radiative processes not considered in the stan- 
dard recombination model of Seager et al. (2000) provide the channels nec- 
essary for a pre-recombination peak. In the recent studies which include ad- 



dditional continuum opacity of Hi in the Hei recombination 1 25 , p5[ evolution, 
the extra high energy distortion photons from He i recombination can excite the 
electrons in the ground state for Hi atom before the recombination of Hi at 
z = 1600-2200. This allows for direct ionization from the ground state, and a 
narrow 'pre-recombination peak' in the Hi Lya line is formed at z~ 1870 (see 
Figure 9 in |4^). This effect is, however, much smaller than the previous false 
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Figure 3.14: The difference between the ratio n2p/n\s and its Bohzmann value 
as a function of redshift z. The sohd hue is A(r^2p/f^ls)""™ from the numerical 
recombination code, while the dashed line i s the approximate value of the actual 
difference A(n2p/nis)''''^' from Equation ( 3.40 ). The dotted (red) line is the 



estimated numerical error ei2(f^2p/"-is)'^"''^ with ei2 = 10 
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signal. The amplitude of this pre-recombination peak is about an order of 
magnitude smaller than the main peak of the H i Ly a line formed during H i 
recombination. 

To conclude, the pre-recombination peaks found in our previous studies 
seem to have been false signals coming from a systematic error in the nu- 
merical code (for both H and He). By including the direct ionization from the 
ground state of H i due to the distortion photons from He i recombination, we 
can find a pre-recombination peak in the H i Ly a line, but with a much reduced 
amplitude. 
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Chapter 4 

Forbidden transitions^ 



4.1 Introduction 

The release of the third year data from the Wilkinson Microwave Anisotropy 
Probe (WMAP) has further improved the precision with which we can con- 
strain the cosmological parameters from the shape of the Cosmic Microwave 
Background (CMB) anisotropics C£[^. The Planck satellite, scheduled for 
launch in 2008 , will provide even higher precision Ci values and data down 
to smaller angular scales {£ < 2500). Higher precision in the observations re- 
quires increased accurarcy from the theoretical calculations, in order for the 
correct cosmological parameters to be extracted. It now seems crucial to obtain 
the CfS down to at least the 1 percent level over a wide range of £. 

CMBFAST ||2^ is the most commonly used Boltzmann code for calculating the 
CiS, and it gives consistent results with other independent codes (see and 
references therein). The dominant uncertainty in obtaining accurate ds comes 
from details in the physics of recombination, for example, the 'fudge factor' in 
the RECFAST routine [E3, M . Calculations of cosmological recombination were 



first published by Peebles (1968) ^ and Zeldovich et al. (1968) |35|. Seager 
et al. (2000) |^ presented the most detailed multi-level calculation and intro- 
duced a fudge factor to reproduce the results within an effective three-level atom 
model. Although the multi-level calculation already gives reasonable accuracy, 
the required level of accuracy continues to increase, so that today any effect 
which is ^ 1 per cent over a range of multipoles is potentially significant. Sev- 
eral modifications have been recently suggested to give per cent level changes in 
the ionization fraction and/or the C^s (see Section 4.4 for details). Most of these 
modifications have been calculated only with an effective three-level code, and 
so the results may be different in the multi-level calculation, since there is no 
thermal equilibrium assumed between the upper states. Here we want to focus 
on one of these modifications, namely the extra forbidden transitions proposed 
by Dubrovich & Grachev (2005) 0, which we study using a multi-level code. 

In the standard calculations of recombination, one considers all the reso- 
nant transitions, but only one forbidden transition, which is the 2S-1S two- 
photon transition, and this can be included for both H and He. Dubrovich & 
Grachev (2005) Q suggested that one should also include the two-photon transi- 
tions from higher excited S and D states to the ground state for H and He i, and 



•^A version of this chapter has been pubhshed: Wong W. Y. and Scott D. (2007) 'The 
effect of forbidden transitions on cosmological hydrogen and helium recombination', Monthly 
Notices of the Royal Astronomical Society, 375, 1441-1448. 
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also the spin-forbidden transition between the triplet and singlet ground 

state l-'^So for Hei. They showed that the recombination of both Hi and He I 
sped up in the three-level atom model. The suggested level of change is large 
enough to bias the determination of the cosmological parameters ||l^ . 

In this chapter we try to investigate the effect of the extra forbidden tran- 
sitions suggested by Dubrovich & Grachev (2005) [Q in the multi-level atom 
model without assuming thermal equilibrium among the higher excited states. 
The outline of this chatper is as follows. In Section|4.2|i 



we will describe details of 



the rate equations in our numerical model. In Section 4.3 we will present results 
on the ionization fraction and the anisotropies Cg, and assess the importance 
of the addition of the forbidden transitions. Other possi ble i mprovements of 
the rec omb ination calculation will be discussed in Section 4.4. And finally in 
Section 4.5 we will present our conclusions. 



4.2 Model 



Here we follow the formalism of the multi-level calculation performed by Seager 
et al. (2000) ||. We consider 100 levels for Hi, 103 levels for Hei, 10 levels for 
Hell, 1 level for Hem, 1 level for the electrons and 1 level for the protons. For 
Hi, we only consider discrete n levels and assume that the angular sub- levels 
(^-states) are in statistical equilibrium within a given shell. For both He i and 
Hell, we consider all the Z-states separately. The multi-level Hei atom includes 
all states with n < 10 and I <7. Here we give a summary of the rate equations 
for the number density of each energy level i, and the equation for the change 
of matter temperature Tm- The rate equation for each state with respect to 
redshift z is 



1 



N 

{ncHcRci - UiRic) + ^ Ai?j_ 



(4.1) 



where is the number density of the ith excited atomic state, ric is the number 
density of free electrons, and ric is the number density of continuum particles 
such as a proton, Hc+, or Hc^+. Additionally Rd is the photo-recombination 
rate, Ric is the photo-ionization rate, ARj^i is the net bound-bound rate for 
each line transition, and H{z) is the Hubble parameter. We do not include the 
coUisional rates, as they have been shown to be negligible ||25||. 

For Hei, we update the atomic data for the energy levels |18||, the oscillator 
strength for resonant transitions ^and the photo-ionization cross-section spec- 
trum. We use the photo-ionization cross-section given by Hummer & Storey 
(1998) [|l^ for n < 10 and I < 4, and adopt the hydrogenic approximation for 
states with I > 5 . It is hard to find published accurate and complete data 
for the photo-ionization cross-section of Hei with large n and I. For example, a 
recent paper by Bauman et al. (2005) ||l| claimed that they had calculated the 
photo-ionization cross-section up to n = 27 and I = 26, although, no numerical 
values were provided. 
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For the matter temperture Tm, we only include the adiabatic and Comp- 
ton cooling terms in the rate and it is given by Equation (^^). Seager et 
al. (2000) 1^ considered all the resonant transitions and only one forbidden 
transition, namely the 2S-1S two-photon transition, in the calculation of each 
atom, (for Hei, 2S= 2^80 and 1S= l^So). The 2S-1S two-photon transition 
rate is given by 

Ai?2s-.is - A2S-1S (n2S - nis^e-''-''^^"^^/'^-^"') , (4.2) 

where A2s_is is the spontaneous rate of the corresponding two-photon transi- 
tion, j^2S-is is the frequency between levels 2S and IS and gt is the degeneracy 
of the energy level i. Here we include the following extra forbidden transitions, 
which were first suggested by Dubrovich & Grachev (2005) 0. The first ones are 
the two-photon transitions from nS and nD to IS for H, plus n^So and n^D2 to 
l^So for He I. For example, for Hi, we can group together the nS and nD states 
coming from the same level, so that we can write the two-photon transition rate 
as 

^-"nS+riD^lS " ^^nS+riD "nS+riD "-IS 6 ■ 

V 5lS J 



Here n (without a subscript) is the principle quantum number of the state, 
nnS+nU is the total number density of the excited atoms in either the nS or nD 
states, and A^g_|_„Q is the effective spontaneous rate of the two-photon transition 
from nS -|- nD to IS, which is approximated by the following formula 



where A^g_.^g is equal to 8.2290 s^^ The latest value of A^g_j^g is equal 

to 8.2206 s~^ jl^ and does not bring any noticeable change to the result. Here 
g-nS+nD is cqual to 1 for n = 2, and 6 for n > 3. This spontaneous rate is es- 
timated by considering only the non-resonant two-photon transitions through 
one intermediate state nP. Dubrovich & Grachev (2005) ignored the resonant 
two-photon transition contributions, since the escape probability of these emit- 
ted photons is very low. The above formula for A^g^^p is valid up to n~40, 
due to the dipole approximation used, although it is not trivial to check how 
good this approximate rate is. Besides the 2S--1S two-photon rate, only the non- 
resonant two-photon rates from 3S to IS and 3D to IS are calculated accurately 
and available in the literature. Cresser et al. (1986) Q evaluated A^g and A^j-, 
by including the non-resonant transitions through the higher-lying intermediate 
nP states (n>4), which are equal to 8.2197s^^ and 0.13171 s""'^, respectively. 
These values were confirmed by Florescu (1988) H and agreed to three signif- 
icant figures. Using these values, we find that A^g_|_„D is equal to 1.484s~^, 
which is an order of magnitude smaller than the value from the approximated 
rate coming from equation (4.4). The approximation given by Dubrovich & 
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Grachev (2005) Q therefore seems to be an overesthxiate. This leads us mstead 
to consider a scaled rate ^nS+nD' which is equal to A^g_|_^j-, multiplied by a 
factor to bring the approximated two-photon rates of Hi (equation (4.4)) with 
n = 3 into agreement with the numerical value given above, i.e. 

~^ns+nu - 0.0664 A«s+„D- (4.5) 

Note that the use of the non-resonant rates is an approximation. The res- 
onant contributions are suppressed in practice because of optical depth effects, 
and in a sense some of these contributions are already included in our multi- 
level calculation. Nevertheless, the correct way to treat these effects would be 
in a full radiative transfer calculation, which we leave for a future study. For 
He I, we treat n^So and n^D2 separately and use equation (4.3) for calculating 
the transition rates. The spontaneous rate A^g^^j^ is estimated by Dubrovich 
& Grachev (2005) § by assuming a similar form to that used for A^g_|_„j-,: 

Ks/nu = — — T > (4-6) 

where A^^^ is a fitting parameter (which is still uncertain both theoretically and 
experimentally). According to Dubrovich & Grachev (2005) 0, resonable values 
of A range from 10 to 12 s~^, and we take A = 11 s~^ here. In our calculation, 
we include these extra two-photon rates up to n = 40 for H and up to n = 10 
for He I. 

The other additional channel included is the spin-forbidden transition be- 
tween the triplet 2'^Pi and singlet I^Sq states in Hei. This is an intercombina- 
tion/ semi- forbidden electric-dipole transition which emits a single photon and 
therefore we can calculate the corresponding net rate by using the bound-bound 
resonant rate expression, i.e. 

Ai?23Pi_iiS(, = P23Pi,liSo ("-23Pi-R23Pi,liSo - '^2iSo-RliSo,23Pi) , (4-7) 



where 



-R23Pi,liSo = ^23Pi-liSo + -S23Pi-liSo'^7 

1 - e-^" 



P23Pi,liSo 



^23Pi-liSoA23p^^liSo 

^ttH(z) 



with 

523 Pi 



5iiSo 



(4.8) 
(4.9) 

(4.10) 
(4.11) 



Here yl23pj_iig^, -B23pi-iiSo a-nd -BiiSo-23pi are the Einstein coefficients, p23Pi,iiSo 
is the Sobolev escape probability, Tg is the Sobolev optical depth (see |2^ and 
references therein), A23pj^ jigj, is the wavelength of the energy difference between 
states 2'^Pi and l^So, and J is the blackbody intensity with temperature Tr. 
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This 2'^Pi~l^So transition is not the lowest transition between the singlet and 
the triplet states. The lowest one is the magnetic-dipole transition between 
and l^So, with Einstein coefficient ^2='Si-iiSo ~ ^■'^•^ ^ 10^^ s^^ However, 
this is much smaller than ^23pi-iiSo ^ 177. 58s^^ |l4|, so this transition can 
be neglected. Note that Dubrovich & Grachev (2005)1^ used an older value of 
^23Pi-iiSo = 233 s~^ in their calculation. 

We use the Bader-Deuflhard semi-implicit numerical integration scheme (see 
Section 16.6 in [^) to solve the above rate equations. All the numerical results 
are carried out using the ACDM model with cosmological parameters: flh = 0.04; 
Oc=0.2; f7A = 0.76; 17k = 0; lp = 0.24; To = 2.725 K and (consistent 
with those in ||||). 



4.3 Result 

The recombination histories calculated using the previous multi-level code [E5[ 



and the code in this paper are shown in Fig. 4.1, where Xc = ne/nn is the 
ionization fraction relative to hydrogen. As we have included more transitions 
in our model, and these give electrons more channels to cascade down to the 
ground state, we expect the overall recombination rate to speed up, and that this 
will be noticeable if the rates of the extra forbidden transitions are significant. 



From Fig. 4.1, we can see that the recombination to Hei is discernibly faster 



in the new calculation. Fig. 4.2 shows the difference in Xc with and without 
the extra forbidden transitions. The dip at around z = 1800 corresponds to 
the recombination of He i and the one around z = 1200 is for H i. Overall, the 
addition of the forbidden transitions claimed by Dubrovich & Grachev (2005) 
leads to greater than 1% change in Xc over the redshift range where the CMS 
photons are last scattering. 

In the last Section, we found that the approximated two-photon rate given 
by Dubrovich & Grachev (2005) for H i with n = 3 was overestimated by more 
than a factor of 10. By considering only this extra two-photon transition, the 
approximate rate gives more than a per cent difference in Xo, while with the 
more accurate nu merical rates, the change in Xc is less than 0.1 per cent (as 



shown in Fig. [4.31) . Based on this result, we do not need to include this two- 
photon transition, as it brings much less than a per cent effect on Xc- For 
estimating the effect of the extra two-photon transitions for higher n, we use 
the scaled two-photon rate given by equation (4.6). The result is plotted in 



Fig. 4.4. The change in Xe with the scaled two-photon rates is no more than 0.4 
per cent, while the one with the Dubrovich & Grachev (2005) approximated 
rates brings about a 5 per cent change. For He i, Dubrovich & Grachev (2005) 
included the two-photon transitions from n — 6 to 40, since they claimed that 
the approximate formula (equation 4.7) is good for n>6. In our calculation, 
we use Ajj'g^^j-, from the approximate formula for the two-photon transitions of 
n = 3 to 10, since this is the best one can do for now (and the formula at least 
gives the right order of magnitude). The addition of the singlet-triplet 2'^Pi- 
l^So transition and the n^So-l^So and n^D2-l^So two-photon transitions with 
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Figure 4.1: The ionization fraction function of redshift z. The sohd line 

is calculated using the original multi- level code of Seager et al. (2000) |25), while 
the dashed line includes all the extra forbidden transitions discussed here. 
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Figure 4.2: The fractio nal d ifference ('new' minus 'old') in Xc between the two 
models plotted in Fig. O as a function of redshift z. The solid and dotted 
lines are the models with the two-photon rates for H i given by Dubrovich & 
Grachev (2005) and the scaled one given by equation (4.6), respectively. Both 
curves are calculated using all the He i forbidden transitions as discussed in the 
text. 
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Figure 4.3: Fractional change in Xc with the addition of the two-photon tran- 
sition from 3S and 3D to IS for Hi. The sohd hne is calculated with the 
approximate rate given by Dubrovich & Grachev (2005) while the dashed 
line is calculated with the numerical rates given by Cresser et al. (1986) 
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n = 3— 10 cause more than 1 per cent changes in (as shown in Fig. 2). The 
2^Pi-l^So transition has the biggest effect on Xe- 



Fig. 4.5 shows the fractional difference in Xc using different combinations 
of additional forbidden transitions. We can see that the 2^Pi-l-'^So transition 
alone causes more than a 1 per cent change in Xc, and the addition of each 
two-photon transition only gives about another 0.1 per cent change. The extra 
two-photon transitions from higher excited states (larger n) have a lower effect 
on Xc compared with that from small n, and we checked that this trend continues 
to higher n. However, the convergence is slow with increasing n. Therefore, one 
should also consider these two-photon transitions with ti > 10 for He I, and the 
precise result will require the use of accurate rates, rather than an approximate 
formula such as equation (4.7). For the 2'^Pi-l^So transition, Dubrovich & 
Grachev (2005) adopted an older and slightly larger rate, and this causes a 
larger change of the ionization fraction (about 0.5 per cent more compared with 



that calculated with our best rate), as shown in Fig. 4.6 



4.3.1 The importance of the forbidden transitions 

One might wonder why the semi-forbidden transitions are significant in recom- 
bination at all, since the spontaneous rate (or the Einstein A coefficient) of 
the semi- forbidden transitions are about 6 orders of magnitude (a factor of , 
where a is the fine-structure constant) smaller than those of the resonant tran- 
sitions. Let us take He i as an example for explaining the importance of the 
spin-forbidden 2'^Pi~l^So transition in recombination. The spontaneous rate 
is equal to 177.58 s^^ for this semi-forbidden transition, which is much smaller 
than 1.7989 x 10^ s~^ for the 2^Pi-l^So resonant transition. But when we cal- 
culate the net rate [see equation (4.7)], we also need to include the effect of 
absorption of the emitted photons by the surrounding neutral atoms, and we 
take this into account by multiplying the net bound-bound rate by the Sobolev 
escape probability pij |p5|| . If pij — 1 , the emitted line photons can escape to in- 
finity, while if Pij = the photons will all be reabsorbed and the line is optically 



thick. Fig. 4.7 shows that the escape probability of the 2-'^Pi-l^So resonant 



transition is about 7 orders of magnitude smaller than the spin-forbidden tran- 



sition. This makes the two net rates roughly comparable, as shown in Fig. 4.8 
From equation (4.11), we can see that the easier it is to emit a photon, the easier 
that photon can be re-absorbed, because the optical depth Tg is directly pro- 
portional to the Einstein A coefficient. So when radiative effects dominate, it is 
actually natural to expect that some forbidden transitions might be important 
(although this is not true in a regime where coUisonal rates dominate which is 
often the case in astrophysics). In fact for today's standard cosmological model, 
slightly more than half of all the hydrogen atoms in the Universe recombined 
via a forbidden transition |Q. Table 4.1 shows that this is also true for helium. 

In the previous multi-level calculation ||2^ , there was no direct transition 
between the singlet and triplet states. The only communication between them 
was via the continuum, through the photo-ionization and photo-recombination 
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transitions. Table 4.1 shows how many electrons cascade down through each 
channel from n — 2 states to the ground state. In the previous calculation, about 
70% of the electrons went down through the 2^Pi-l^So resonant transition. In 
the new calculation, including the spin-forbidden transition between the triplets 
and singlets, there are approximately the same fraction of electrons going from 
the 2^Pi and states to the ground state (actually slightly more going from 
2^Pi in the current cosmological model). This shows that we should certainly 
include this forbidden transition in future calculations. Our estimate is that 
only about 40% of helium atoms reach the ground state without going through 
a forbidden transition. 

How about the effect of other forbidden transitions in Hei recombination? 
We have included all the semi-forbidden electric-dipole transitions with n < 10 
and I < 7, and with oscillator strengths larger than 10~^ given by Drake & Mor- 
ton (2007) There is no significant change found in the ionization fraction. 
Besides the 2'^Pi-l^So transition, all the other extra semi-forbidden transitions 
are among the higher excited states where the resonant transitions dominate. 
This is because these transition lines are optically thin and the escape proba- 
bilities are close to 1. 



4.3.2 Effects on the anisotropy power spectrum 

The CMB anisotropy power spectrum Cg depends on the detailed profile of the 
evolution of the ionization fraction Xe- This determines the thickness of the 
photon last scattering surface, through the visibility function g{z) = e~'^ dr/dz, 
where r is the Thomson scattering optical depth {t ^ ca^ J nc{dt/dz) dz). The 
function Xc{z) sets the epoch when the tight coupling between baryons and 
photons breaks down, i.e. when the photon diffusion length becomes long, and 
the visibility function fixes the time when the fluctuations are effectively frozen 



in (see |ll|, 25 and references therein). The addition of the extra forbidden 
transitions speeds up both the recombination of H i and He i, and hence we 
expect that there will be changes in Ci. 

In order to perform the required calculation, we have used the code CMB- 
FAST and modified it to allow the input of an arbitrary recombination his- 



tory. Figs. 4.9 and 4.10 show the relative changes in the CMB temperature (TT) 
and polarization (EE) anisotropy spectra, respectively, with different combina- 
tions of extra forbidden transitions. The overall decrease of free electrons brings 
a suppression of Cg over a wide range of £. 

For He I, there is less at 2;~1400 — 2500, which leads to an earlier re- 
laxation of tight coupling. Therefore, both the photon mean free path and the 
diffusion length are longer. Moreover, the decrease of Xc in the high-z tail results 
in increased damping, since the effective damping scale is an average over the 
visibility function. This larger damping scale leads to suppression of the high-^ 



part of the power spectrum. From Figs. and [4.10| , we can see a decrease 
of Ci (for both TT and EE) toward high £ for Hei, with the maximum change 
being about 0.6 percent. 

For H I, the change of Cg is due to the decrease in Xc at z ~ 600 — 1400 (see 
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Table 4.1: The percentage of electrons cascading down in each chan- 
nel from n = 2 states to the I^Sq ground state for He I. 

2%^=ri% 2Tp7^7TTs;; 2^pr=ri%~ 

(two-photon) (resonant) (spin-forbidden) 
Seager et al. (2000) 30.9% 69.1% 

this work 17.3% 39.9% 42.8% 




Figure 4.4: Fractional change in with the addition of different forbidden 
transitions for H I. The long-dashed, dotted, dashed and solid lines include the 
two-photon transitions up to n = 10, 20, 30 and 40, respectively, using the 
approximation for the rates given by equation (4.5). The dot-dashed line is 
calculated with the scaled rate from equation (4.6). 
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Figure 4.5: Fractional change in Xe with the addition of different forbidden 
transitions for He I as a function of rcdshift. The sohd line corresponds to 
the calculation with only the 2''Pi l^So spin-forbidden transition. The short- 
dashed, dotted, long-dashed, dot-dashed and long dot-dashed lines include both 
the spin-forbidden transition and the two-photon (27) transition(s) up to n = 
3,4, 6, 8 and 10, respectively. 



Chapter 4. Forbidden transitions 



92 




Figure 4.6: Fractional change in Xc with only the Hei forbidden tran- 

sition. The soUd Hne is computed with our best value ^23Pi-iiSo — 177.58 s""'^ 
from Lach & Panchucki (2001) and the dashed line is calculated with the 
rate Aasp^^iig,, = 233 s^^ from Dubrovich & Grachev (2005) 0. 




Figure 4.7: Escape probability pij as a function of redshift. The solid line 
corresponds to the resonant transition between He I 2^Pi and l^So, while the 
dashed line refers to the spin-forbidden transition between Hei and I^Sq- 
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Figure 4.8: Not bound-bound rates for Hel as a funetion of rcdshift. The solid 
line is the resonant transition between 2^ Pi and I^Sq, the short-dashed line is 
the two-photon transition between 2^Si and I^Sq. And the long-dashed line is 
the spin- for bidden transition between 2^ Pi and I^Sq. 
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Fig. 4.2). There are two basic features in the curve of change in Cg (the dotted 
and dashed hues in Fig. |4j| ). Firstly, the power spectrum is suppressed with 
increasing i, due to the lower Xe in the high- 2; tail {z > 1000). Secondly, there are 
a series of wiggles, showing that the locations of the acoustic peaks are slightly 
shifted. This is due to the change in the time of generation of the C^s in the 
low-z tail. Cf^ actually shows an increase for £ < 1000 (see Fig. 4. 10] ); this is 
caused by the shift of the center of the visibility function to higher z, leading to 
a longer diffusion length. Polarization occurs when the anisotropic hot and cold 
photons are scattered by the electrons. The hot and cold photons can interact 
with each other through multiple scatterings within the diffusion length, and 
therefore, a longer diffusion length allows more scatterings and leads to a higher 
intensity of polarization at large scales. 

With the approximate rates used by Dubrovich & Grachev (2005) 0, the 
maximum relative change of Cf^ is about 4 percent and for Cf^ it is about 
6%. The overall change is thus more than 1% over a wide range of £. However, 
if we adopt the scaled two-photon rate given by equation (4.6), the relative 
changes of Cf^ and Cf^ are no more than 1 per cent. Note that we do not plot 
the temperature-polarization correlation power spectrum here, since there is no 
dramatically different change found (and relative differences are less meaningful 
since Cf-^ oscillates around zero). 



4.4 Discussion 



In our model we only consider the semi-forbidden transitions with n< 10 and 
I < 7 for He I and the two-photon transitions from the higher S and D states 
to the ground state for H and Hei. It would be desirable to perform a more 
detailed investigation of all the other forbidden transitions, which may provide 
more paths for the electrons to cascade down to the ground state and speed 
up the recombination process. In this paper we have tried to focus on the 
forbidden transitions which are likely to be the most significant. However we 
caution that, if the approximations used are inadequate, or other transitions 
prove to be important, then our results will not be accurate. 

There are several other approximations that we have adopted in order to per- 
form our calculations. For example, we consider the non-resonant two-photon 
rates for higher excited rates. The two-photon transitions from higher excited 
states {n > 3) to the ground state are more complicated than the 2S-1S transi- 
tion, because of the resonant intermediate states. For example, for the 3S-1S 
two-photon transition, the spectral distribution of the emitted photons shows 
infinities (resonance peaks) at the frequencies corresponding to the 3S-2P and 
2P-1S transitions!^. Here, we use only the non-resonant rates, by assum- 
ing a smooth spectral distribution of the emitted photons; this probably gives 
a lower limit on the change of a^c and Cg coming from these extra forbidden 
transitions. The correct way to treat this would be to consider the rates and 
feedback from medium using the full spectral distribution of the photons and 
radiative transfer; this will have to wait for a future study. 



Chapter 4. Forbidden transitions 



96 




multipole / 

Figure 4.9: Relative change in the temperature (TT) angular power spectrum 
due to the addition of the forbidden transitions. The solid line includes the 
spin-forbidden transition and also the two-photons transitions up to n = 10 for 
He I, the dotted line includes all the above transitions and also the two-photon 
transitions up to rt = 40 for Hi calculated with the approximate rates given 
by Dubrovich & Grachev (2005) §. The dashed line is computed with the 
same forbidden transitions as the dotted line, but with our scaled rates (and 
represents our best current estimate). 
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Figure 4.10: Relative change in the polarization {EE) angular power spectrum 
due to t he addition of the forbidden transitions, with the curves the same as in 
Fig. 1^9. 
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Besides the consideration of more forbidden transitions, there are many other 
improvements that could be made to the recombination calculation by the time 
when this work was published. In particular, Rubino-Marti'n et al. (2006) p2[ 
showed that a multi-level calculation of the recombination of H i with the in- 
clusion of separate Z-states can give more than 20 per cent difference in the 
population of some levels compared with the thermal equilibrium assumption 
for each n-shell. The latest calculation, considering up to 100 shells, is presented 
by Chluba et al. (2006) but does not include all the forbidden transitions 
studied here. A more complete calculation should be done by combining the 
forbidden transitions in a code with full angular momentum states, and we leave 
this to a future study. There are also other elaborations which could be included 
in future calculations, which we now describe briefly. 

The rate equation we use for all the two-photon transitions only includes 
the spontaneous term, assuming there is no interaction with the radiation back- 
ground (see equation (4.3)). Chluba & Sunyaev (2005) |H suggested that one 
should also consider the stimulated effect of the 2S-1S two-photon transition for 
H, due to photons in the low frequency tail of the CMB blackbody spectrum. 
Leung et al. (2004) |^ additionally argued that the change of the adiabatic 
index of the matter should also be included, arising due to the neutralization of 
the ionized gas. These two modifications have been studied only in an effective 
three-level atom model, and more than a percent change in Xe was claimed in 
each case (but see Chapter 5 34 for arguments against the effect claimed by 
Leung et al. 2004 [[| ). 

For the background radiation field J, we approximated it with a perfect 
blackbody Planck spectrum. This approximation is not completely correct for 
the recombination of Hi, since the He line distortion photons rcdshift into a 
frequency range that can in principle photo- ionize the neutral H j^, ^s], |33| . 
Althought we expect this secondary distortion effect to bring the smallest change 
on Xe among all the modifications suggested here, it is nevertheless important 
to carry out the calculation self-consistently, particularly for the spectral line 
distortions. In order to obtain an accurate recombination history, we therefore 
need to perform a full multi-level calculation with seperate Z-states and include 
at least all the improvements suggested above, which we plan to do in a future 
study. 

For completeness we also point out that the accuracy of the physical con- 
stants is important for recombination as well. The most uncertain physical quan- 
tity in the recombination calculation is the gravitational constant G. The value 
of G used previously in the recfast code is 6.67259 x 10^^^m'^kg~^s~^ and the 
latest value (e.g. from the Pariticle Data Group is 6.6742 x 10~^^m'^kg~^s~^. 
Another quantity we need to modify is the atomic mass ratio of *He and ^H, 
TO4fj(,/mifj, which was previously taken to be equal to 4 (as pointed out by 
Steigman 2006 |2^). By using the atomic masses given by Yao et al. (2006) |]3^ , 
the mass ratio is equal to 3.9715. The overall change in a;c is no more than 0.1 
per cent after updating these two constants in both recfast and multi-level 
code. 
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4.5 Conclusion 

In this paper, we have computed the cosmological recombination history by 
using a muhi-level code with the addition of the to l^So spin-forbidden 

transition for He i and the two-photon transitions from nS and nD states to the 
ground state for both H i and He i. With the approximate rates from Dubrovich 
& Grachec (2005) 0, we find that there is more than a per cent decrease in the 
ionization fraction, which agrees broadly with the result they claimed. However, 
the only available accurate numerical value of two-photon rate with n > 3 is 
for the 3S to IS and 3D to IS transitons for H. We found that the approximate 
rates from Dubrovich & Grachec (2005) were overestimated, and instead we 
considered a scaled rate in order to agree with the numerical n — 3 two-photon 
rate. With this scaled rate, the change in Xc is no more than 0.5 per cent. 

Including these extra forbidden transitions, the change in the CMB anisotropy 
power spectrum is more than 1 per cent, which will potentially affect the de- 
termination of cosmological parameters in future CMB experiments. Since one 
would like the level of theoretical uncertainty to be negligible, it is essential to in- 
clude these forbidden transitions in the recombination calculation. In addition, 
we still require accurate spontaneous rates to be calculated for the two-photon 
transitions and also a code which includes at least all the modifications sug- 
gested in Section 4.4, in order to obtain the C^s down to the 1 per cent level. 
Achieving sub-percent accuracy in the calculations is challenging! 

However, the stakes are high - the determination of the parameters which 
describe the entire Universe - and so further work will be necessary. Systematic 
deviations of the sort we have shown would potentially lead to incorrect values 
for the spectral tilt derived from Planck and even more ambitions future CMB 
experiments, and hence incorrect inferences about the physics which produced 
the density perturbations in the very early Universe. It is amusing that in order 
to understand physics at the 10^^ GeV energy scale we need to understand eV 
scale physics in exquisite detail! 



Chapter 4. Forbidden transitions 



100 



4.6 References 

[1] Bauman R. P., Porter R. L., Ferland G. J., MacAdam K. B. 2005, Astro- 
physical Journal, 628, 541 

[2] Chluba J., Sunyaev R. A. 2006, Astronomy and Astrophysics, 446, 39 

[3] Chluba J., Rubino-Martm J. A., Sunyaev R. A. 2006, Monthly Notices of 
the Royal Astronomical Society, 374, 1310 

[4] Cresser J. D., Tang A. Z., Salamo G. J., Chan F. T. 1986, Physical Review 

A, 33, 3, 1677 

[5] Dell'Antonio I. P., Rybicki G. B. 1993, in ASP Conf. Ser. 51, Observational 

Cosmology, cd. G. Chincarini ct al. (San Francisco: ASP), 548 

[6] Drake G. W. F., Morton, D. C. 2007, Astrophysical Journal Supplement, 
170, 251 

[7] Dubrovich V. K., Grachev S. I. 2005, Astronomy Letters, 31, 6, 35 

[8] Florescu V., Patrascu S., Stoican O. 1987, Physical Review A, 36, 2155 

[9] Florescu V., Schneider I., Mihailescu I. N. 1988, Physical Review A, 38, 4, 
2189 

[10] Goldman S. P. 1989, Physical Review A, 40, 1185 

[11] Hu W., Scott D., Sugiyama N., White M. 1995, Physical Review D, 52, 
5498 

[12] Hummer D. G., Storey P. J. 1998, Monthly Notices of the Royal Astronom- 
ical Society, 297, 1073 

[13] Labzowsky L. N., Shonin A. V., Solovyev D. A. 2005, Journal of Physics 

B, 38, 265 

[14] Lach G., Pachucki K. 2001, Physical Review A, 64, 042510 

[15] Leung P. K., Chan C. W., Chu M. C. 2004, Monthly Notices of the Royal 

Astronomical Society. 349, 2, 632 

[16] Lewis A., Weller J., Battye R. 2006, Monthly Notices of the Royal Astro- 
nomical Society, 373, 561 

[17] Lin C. D., Johnson W. R., Dalgarno A. 1977, Physical Review A, 15, 1, 
154 

[18] Morton D. C, Wu Q., Drake G. W. F. 2006, Canadian Journal of Physics, 
84, 83 

[19] Peebles P. J. E. 1968, Astrophysical Journal, 153, 1 



Chapter 4. Forbidden transitions 



101 



[20] Planck Collaboration 2006, ESA-SCI(2005)1, arXiv:astro-pli/0604069 

[21] Press W.H., Flannery B. P., Teukolsky S. A., Vetterling W. T. 1992, Nu- 
merical Recipes in C: The Art of Scientific Computing, Cambridge Univ. 

Press, Cambridge, UK 

[22] Rubiiio-Martfn J. A., Chluba J., Sunyaev R. A. 2006, Monthly Notices of 
the Royal Astronomical Society, 371, 1939 

[23] Santos J. P., Parente F., Indelicato P. 1998, European Physical Journal, 
D3, 43 

[24] Seager S., Sasselov D. D., Scott D. 1999, Astrophysical Journal, 523, LI 

[25] Seager S., Sasselov D. D., Scott D. 2000, Astrophysical Journal Supple- 
ment, 128, 407 

[26] Seljak U., Sugiyama N., White M., Zaldarriaga M. 2003, Physical Review 
D, 68, 083507 

[27] Seljak U., Zaldarriaga M. 1996, Astrophysical Journal, 463, 1 

[28] Spergel D. N. el al. 2006, Astrophysical Journal Supplement, 170, 377 

[29] Steigman G. 2006, Journal of Cosmology and Astro-Particle Physics, 10, 
16 

[30] Storey P. J., Hummer D. G. 1991, Computer Physics Communications, 66, 
12 

[31] Tung J. H., Ye X. M., Salamo G. J., Chan F. T. 1984, Physical Review A, 
30, 1175 

[32] Yao W.-M. et al., 2006, Journal of Physics G, 33, 1 

[33] Wong W. Y., Seager S., Scott D. 2006, Monthly Notices of the Royal As- 
tronomical Society, 367, 1666 

[34] Wong W. Y., Scott D. 2006, ArXiv e-prints, astro- ph/0612322 

[35] Zel'dovich Y. B., Kurt V. G., Sunyaev R. A. 1968, Zh. Eksp. Teor. Fiz., 
55, 278; English translation, 1969, Soviet Phys. JETP Lett., 28, 146 



102 



Chapter 5 

Matter temperature^ 



5.1 Introduction 

Detailed calculations of the process through which the early Universe ceased 
to be a plasma are increasingly important because of the growing precision of 
microwave anisotropy experiments. The standard way to calculate the evolution 
of the matter temperature during the process of cosmological recombination 
is to consider the expansion of radiation and matter separately, and include 



the relevant interactions, specifically Compton scattering (see Equation (3.9)) 
and photoionization cooling, as corrections |[ The matter is treated as 
a perfect gas which is assumed to envolve adiabatically. Recently, Leung et 
al. (2004) D suggested that we need to use a generalized adiabatic index, since 
the gas was initially ionized and so the number of species (ion + electron vs 
atom) changes in the recombination process. In their derivation, they considered 
the effect of photoionization, recombination and excitation on the matter, but 
assumed that the matter was undergoing an adiabatic process. However, an 
adiabatic approximation for only the ionized matter is not valid in this case, 
because the change of entropy of the matter is not zero. Moreover, the photons 
released from the recombination of atoms mostly escape as free radiation JlTI , 
instead of reheating the matter, since the heat capacity of the radiation is much 
larger than that of the matter (see, for example, |^ H). 

Here, wc try to study this problem in a consistent way by considering both 
the radiation and ionizing hydrogen as components in thermal equilibium and 
under adiabatic expansion. This is not exactly the way things happened during 
recombination, but this will give us the maximum effect of the heat if it is all 
shared by the radiation and matter. 



5.2 Discussion 

For simplicity, we consider that the matter consists only of hydrogen (including 
helium docs not change the physical picture). By assuming that the radiation 
field and the matter are in thermal equilibrium, the total internal energy per 

*A version of this chapter has been posted on the e-prints ArXiv: Wong W. Y. and 
Scott D. (2006) 'Comment on "Recombination induced softening and reheating of the cos- 
mic plasma"', ArXiv e-prints, arXiv:astro-ph/0612322. 
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unit mass of the system is 
1 



TlHTOH 

1 



07"* + x("-H + nc)k^T + npe^„ + ^ nf ef 



(5.1) 



where is the number density of neutral atoms in the zth state, rip is the 



number density of free protons, nn 



- nf is the total number density of 



neutral and ionized hydrogens, and is the number density of free electrons. 
Additionally ePjj and ej^ are the ionization energy for the ground state and the 
ith state of hydrogen, respectively, the xs are the fractional number densities 
normalized by nn, T is the temperature of the whole system, a is the radiation 
constant and fee is Boltzmann's constant. 

In equation (5.1), the first term is the radiation energy, the second term 



is the kinetic energy of the matter, and the last two terms are the excitation 
energy of the atoms. Here the energy of the ground state is set to be equal 
to zero No matter what energy reference is chosen, the change of energy 
should be the same, i.e. 



dE 



int 



1 



4aT3 



-dT 



?1H 



dnn + -(1 + Xc)k^dT 



--ksTdxc + e^ndxp 



(5.2) 



We know that the radiation and matter are not exactly in thermal equilibrium 
(the two temperatures are not precisely the same) during the cosmological re- 
combination of hydrogen, because the recombination rate is faster than the rate 
of expansion and cooling of the Universe. Nevertheless, the radiation back- 
ground and matter are tightly coupled and it is a good approximation to treat 
the two as if they were in thermal equilibrium (for example, Peebles 1971 
P.232). This simple approach allows us to estimate how much of the heat re- 
leased is shared with the radiation field and the matter during the recombination 
of hydrogen. In the expansion of the Universe, the whole system (radiation plus 
matter) is under an adiabatic process. However, this is not the case for the 
ionizing matter on its own, because the change of entropy of the matter is not 
zero. For an adiabatic process we have 



dE 



int 



p 



dp 



1 



{l + X,)kBT+- 

3 nn 



3dz 



(5.3) 
(5.4) 



where P and p are the pressure and mass density of the system and z is redshift. 
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By equating equations (^) and {pA\), we have 
1 + zdT 3{l + x,)kBT+^ 



1 + z 



T dz 



f(l + :.e)fcBT+i^ 



f(l+Xe)fcBT+itf^ 



(5.5) 



In order to see whether we can ignore the radiation field, we need to compare the 
two terms in the denominator, i.e. the radiation energy and the kinetic energy 
of the matter. If the matter energy were much greater than the radiation energy, 
then we would have 



l + zdT 
T dz 



2 - 



1 



f(l+Xe)fcBr 



3, rr.dXc 



H dxp 



,y.H 

dz 



(5.6) 



which is the result given by Leung et al. (2004) |]. However, in the current 
cosmological model with Tq = 2.725, Yp = 0.24, h = 0.73 and Ob = 0.04 (for 
example, 1^), we have 



matter 
TPint 
radiation 



1.6 X 10" 



(5.7) 



So, the radiation energy is much larger than both the kinetic energy of matter 
and also the total heat released during recombination. Hence, we definitely can- 



not ignore the radiation field. In such a case, the second term in equation (5.5) 
is much smaller than the first term, because 



1(1 



Xc)kT - 



(1 + z)nHeP„ 



10- 



(5.8) 



Hence the energy change due to the recombination process is taken up mostly by 
the radiation field, since there are many more photons than baryons. In other 
words, most of the extra photons (or heat) escape to the photon field, with just 
a very small portion (^ 10^^°) reheating the matter. Therefore, the change of 
the temperature of the system can be approximated as 



l + zdT 
T dz 



1±6, 



(5.9) 



where 5 < 10~^. This gives us back the usual formula for the radiation tempera- 
ture, which is consistent with the result that the matter temperature closely fol- 
lows the radiation temperatre (for example, |§, ^). Leung, Chan & Chu (2004) |^ 
assumed that the extra heat is shared by the matter only, and hence that the 
second term of equation (5.6) is significant, because 



(1 + ^)^Pon 
{l + Xe)kT 



6 X 10^ 



(5.10) 
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By comparing this and the ratio given in equation (5.8), we can see that the 
factor is about 10 orders of magnitude larger if we ignore the radiation field. An- 
other way to understand this overestimate is that the adiabatic approximation 
for matter only is not valid in their derivation, because the entropy of the matter 
is changing (i.e. dSma.tter/dz > 0). The Leung et al. (2004) paper ignored the 
last term (the sum of the excitation energy terms) in equation (5.6), which phys- 
ically means that there is a photon with energy equal to (^13.6 eV) emitted 
when a proton and electron recombine, and the energy of this distortion photon 
is used up to heat the matter. This is actually not true for the recombination 
of hydrogen, since there is no direct recombination to the ground state jj, [ij, || 
and there are about 5 photons per neutral hydrogen atom produced for each 
recombination 0. 

Note that what we calculate above is in the thermal equilibrium limit and it 
assumes that all the distortion photons are thermalized with the radiation back- 
ground and the matter. However, in the standard recombination calculation, 
most of these distortion photons escape to infinity with tiny energy loss to the 
matter through Compton scattering. The maximum fraction of energy loss by 
the distortion photons after multiple scatterings {AE^/E^) is very low [|l0|. An 
approximate estimate is 

A P 



Ej rric 



13 6 eV 

X 30 at z ~ 1500, when ~ 0.9 



511 keV 

8 X 10^^ 



where rric is the mass of electron, c is the speed of light and r is the optical 
depth. Therefore, the e^^dxp/dz term is in practise suppressed by at least 10~^ 
(since r decreases when more neutal hydrogen atoms form at lower redshift). 
Hence, although there is some heating of the matter, the ratio of the heat 
shared by the matter and the radiation is very small, and the effect claimed by 
Leung et al. (2004) |] is negligible for the recombination history and also for 
the microwave anisotropy power spectra. 



5.3 Conclusion 

By considering a simple model consisting of the radiation background and the 
ionizing gas under equilibrim adiabatic expansion, we show that the effect 
claimed by Leung et al. (2004) |^ is hugely overestimated. The appropriate 
method for calculating the matter temperature is to deal with Compton and 
Thomson scattering between the background photons, distortion photons and 
matter in detail. In general the Compton cooling time of the baryons off the 
CMB is very much shorter than the Hubble time until z ~ 200, hence it is 
extremely hard for any heating process to make the matter and radiation tem- 
peratures differ significantly at much earlier times. 
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How well do we understand 
recombination?^ 

6.1 Introduction 

Planck the third generation Cosmic Microwave Background (CMB) satelUte 
will be launched in 2008; it will measure the CMB temperature and polarization 
anisotropics Cg at multipoles ^ = 1 to ~ 2500 at much higher precision than 
has been possible before. In order to interpret these high fidelity experimental 
data, we need to have a correspondingly high precision theory. Understanding 
precise details of the recombination history is the major limiting factor in cal- 
culating the to better than 1 per cent accuracy. An assessment of the level 
of this uncertainty, in the context of the expected Planck capabilities, will be 
the subject of this chapter. 

The general physical picture of cosniological recombination was first given by 
Peebles (1968) |l|] and Zeldovich et al. (1968) They adopted a simple three- 
level atom model for hydrogen (H) , with a consideration of the Ly a and lowest 
order 2s-ls two-photon rates. Thirty years later, Seager et al. jl^ performed 
a detailed calculation by following all the resonant transitions and the lowest 
two-photon transition in multi-level atoms for both hydrogen and helium in a 
blackbody radiation background. Lewis et al. (2006) first discussed how 
the uncertainties in recombination might bias the constraints on cosmological 
parameters coming from Planck; this study was mainly motivated by the effect 
of including the semi- forbidden and high-order two-photon transitions Q , which 
had been ignored in earlier calculations. 

There have been many updates and improvements in the modelling of re- 
combination since then. Switzer & Hirata (2008) |]l8) presented a multi-level 
calculation for neutral helium (He i) recombination including evolution of the 
radiation field, which had been assumed to be a perfect blackbody in previ- 
ous studies. Other issues discussed recently include the continuum opacity due 
to neutral hydrogen (Hi) (see also [Q), the semi-forbidden transition 2'^p-l-'^s 
(the possible importance of which was first proposed by |^ ) , the feedback from 
spectral distortions between 2^p-l^s and 2'^p-l^s lines, and the radiative line 
transfer. In particular, continuum absorption of the 2^p-l^s line photons by 

version of this chapter has been pubhshed: Wong W. Y., Moss A. and Scott D. (2008) 
'How well do we understand cosmological recombination?', Monthly Notices of the Royal 
Astronomical Society, 386, 1023-1028. 



Chapter 6. How well do we understand cosmological recombination? 108 



neutral hydrogen causes helium recombination to end earlier than previously 
estimated (see Fig. 6.1). Hirata & Switzer (2008) |^ also found that the high 
order two-photon rates have a negligible effect on Hei, and the same conclu- 
sion was made by other groups for hydrogen as well ^ , largely because the 
approximate rates adopted by Dubrovich & Grachev (2005) Q had been over- 
estimated. The biggest remaining uncertainty in He i recombination is the rate 
of the 2^p-l^s transition, which causes a variation equal to about 0.1 per cent 
in the ionization fraction p9| | . 

For hydrogen, Chluba et al.(2006) |^ improved the multi- level calculation by 
considering seperate angular momentum £ states. This brings about a 0.6 per 
cent change in Xc at the peak of the visibility function, and about 1 per cent at 
redshifts z < 900. The effect of the induced 2s-ls two-photon rate due to the 
radiation background [Q is partially compensated by the feedback of the Ly a 
photons and the net maximum effect on Xc is only 0.55 per cent at z ~ 900. 
The high-order two photon transitions bring about a 0.4 per cent change in Xg 
at 2 ~ 1160 jj, |l|. There are also 0.22 per cent changes in Xc at z ~ 1050 when 
one considers the Lyman series feedback up to n = 30, and there is additionally 
the possibility of direct recombination, although this has only a roughly 10""* 
per cent effect Q . 

The list of suggested updates on Xe is certainly not complete yet, since some 
additional effects, such as the convergence of including higher excited states 
and the feedback-induced corrections due to the He I spectral distortions, may 
enhance or cancel other effects. In general we still need to develop a complete 
multi-level code for hydrogen, with detailed interactions between the atoms and 
the radiation field. However, what is really important here is establishing how 
these effects propagate into possible systematic uncertainties in the estimation 
of cosmological parameters. 

Since the uncertainties in cosmological recombination discussed in the Lewis 
et al. (2006) ||l^ paper have been reduced or updated, it is time to revisit the 
topic on how the new effects or remaining uncertainties might affect the con- 
straints on cosmological parameters in future experiments. The recent full ver- 
sion of the He i recombination calculation j^, |lj, 1^ takes too long to run to be 
included within the current Boltzmann codes for C{. So instead, in this paper, 
we try to reproduce the updated ionization history by modifying recfast |^ 
using a simple parametrization based on the fitting formulae provided by Kholu- 
penko et al. (2007) [1). We then use the COSMOMC jll| code to investigate how 
much this impacts the constraints on cosmological parameters for an experiment 
like Planck. 



6.2 Recombination model 

In this paper, we modify recfast based on the fitting formulae given by Kholu- 
penko et al. (2007) |^ for including the effect of the continuum opacity of neutral 
hydrogen for He i recombination. The basis set of rate equations of the ioniza- 
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Figure 6.1: Top panel: Ionization fraction function of redshift z. 

The dotted (red) line is calculated using the original RECFAST code. The 
solid (black) line is the numerical result from ||l^ , while the dashed (blue) and 
long-dashed (green) lines are both evaluated based on the modification given 
by - the dashed one has 6hc = 0.97 (the value used in the original paper) 
and the long-dashed one has 6hc — 0.86. Bottom panel: The visibility func- 
tion g(z) versus redshift z. The two curves calculated (dotted and long-dashed) 
correspond to the same recombination models in the upper panel. The cosmo- 
logical parameters used for these two graphs are, Sib = 0.04592, fl^ = 0.27, 
Qc ^ 0.22408, To = 2.728 K, Hq = 71 kms^^ Mpc^i and /ho = 0.079. 
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tion fraction of H i and He i used in recfast are: 
dxr 



H{z){l + z)-^ = (:re2:pnHaH-/3H(l-:^p)e-''"'^"^"/'^")CH, (6.1) 



H{Z){1 + ^)^^ = [xKcllXcnKaKel - /?Hel(/He " .THcIl)e-''^'^H=i.^^ = /'-'^^' ) CrcI 

+ fxHeiiXcnH«Li - ^^^^/3Hci(/He " a;Heii)e-'^^''H.i.23./fcTM ^ ^^^^^ (6.2) 

V .9HcL11s / 

where 

^ ^ 1 + KiiAiiniijl - Xp) 

" l + XH(AH + /3H)nH(l-Xp)' ^ ■ ^ 

'''' l + ifHcl(AHc + /3Hcl)nH(.fHc-XHeIl)e"-'^-/'=^-' 

^ l + ncl/3Ll^H(/Hc-:rHeIl)e^--/'=^"- ^^-^^ 

Note that Xe is defined as the ratio of free electons per H atom and so Xo > 1 
during He recombination. We foUow the exact notation used in Seager et 
al. (1999) [|l5) and we do not repeat the definitions of all symbols, except those 
that did not appear in that paper. The last term in equation ( |6.2| ) is added 
to the original dxncii/dz rate for the recombination of Hei through the triplets 
by including the semi- forbidden transition from the 2^p state to the l^s ground 
state. This additional term can be easily derived by considering an extra path 
for electrons to cascade down in He I by going from the continuum through 2'^p 
to ground state, and assuming that the rate of change of the population of the 
2^p state is negligibly small. The superscript 't' stands for triplets, so that, for 
example, a|jgj is the Case B Hei recombination coefficient for triplets. Based 
on the data given by Hummer & Storey (1998) Q, a^gj is fitted with the same 
functional form used for the q;hcI singlets (see equation (4), in jl^), with dif- 
ferent values for the parameters: p — 0.761; q = ^q^^^ '^'^^; Ti = 10^'^"'^"' K; and 
T2 = 3 K. This fit is accurate to better than 1 per cent for temperatures between 
10^'^ and 10^ K. Here /J^gj is the photoionization coefficient for the triplets and 
is calculated from a^gj by 

PHoI - "HoI 72 ^ ■ ' (.O.bj 

V "P / 5HcI,23s 



where (7hc+ ^-nd 5hcI,23s s-re the degeneracies of He 11 and of the He i atom with 
electron in the 2'^s state, and hpi'2is,c is the ionization energy of the 2'^s state. 

The correction factor ChoI accounts for the slow recombination due to the 
bottleneck of the Hei 2^p-l^s transition among singlets. We can also derive the 
corresponding correction factor Cjijgj for the triplets. The K-n, Kuci and K'^^^ 
quantities are the cosmological redshifting of the H Lya, Hci 2^p-l^s and Hei 
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2^p-l^s transition line photons, respectively. The factor K used in recfast is 
a good approximation when the line is optically thick (r ;2> 1) and the Sobolev 
escape probability's is roughly equal to l/r. In general, we can relate K and 
PS through the following equations (taking He I as an example) : 

.9HoI,11s 1 , „^ 

KhcI = and (6.7) 

3HcL2ip ^HoI,lis^2ip-lis^'S 



H 5HcI,11s 1 



KLi = ^^^^ ^ (6-8) 

where ^Hei,2ip-iis and ^Hoi,23p-iis are the Einstein A coefficients of the He I 
2^p-l^s and He I 2'^p-l-^s transitions, respectively. Note that AHei,23p-iis = 
5Hoi,23Pi/5Hci,23p X AhcL23Pi-i1s = 1/3 X 177.58 s"! @. For Hei 2ip-lis, 
we replace ps by the new escape probability Posc, to include the effect of the 
continuum opacity due to H, based on the approximate formula suggested by 
Kholupenko et al. (2007) §|. Exphcitly this is 

Pose = Ps + J3con,H I (6-9) 

where 

1 - e"^ 

PS = and (6.10) 



1 

1 + aHoT""" 



with 

7 - 



87r3/2CTH,ls(l^HoI,2ip)^'HcI,2ip^'^D,2ip(l " Xp) ' 

where crH,is(^'Hoi.2ip) is the Hi ionization cross-section at frequency 1^1101, 2ip and 
^^'D,2ip = ^'Hci,2ip\/2fcB7M/TOHoC^ is the thermal width of the Hei 2^p-l'^s line. 
The 7 factor in pcon,H is approximately the ratio of the Hei 2^p-l^s transition 
rate to the Hi photoionization rate. When 7^1, the effect of the continuum 
opacity due to neutral hydrogen on the Hei recombination is negligible. Here 
aHo and buc are fitting parameters, which arc equal to 0.36 and 0.97, based on 

the results from Kholupenko et al. (2007) j|] . 

We now try to reproduce these results with our modified recfast. Fig. |6.l| 
(upper panel) shows the numerical result of the ionization fraction Xq from differ- 
ent He I recombination calculations. The results from Kholupenko et al. (2007) ^ 
and Switzer & Hirata (2008) both demonstrate a significant speed up of Hei 
recombination compared with the original recfast. We do not expect these 
two curves to match each other, since Kholupenko et al. (2007) 1^ just included 
the effect of the continuum opacity due to hydrogen, which is only one of the 
main improvements stated in Switzer & Hirata (2008) ||l^. Nevertheless, we can 
regard the Kholupenko et al. (2007) study as giving a simple fitting model 
in a three-level atom to account for the speed-up of the He i recombination. 



Fig. 3.2 shows how the ionization history changes with different values of the 
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fitting parameter 6hc (with cho fixed to be 0.36). When ^ho is larger than 1.2, 
the effect of the neutral H is tiny and the fit returns to the situation with no con- 
tinuum opacity. However, if &ho is smaller than 1, the effect of the continuum 
opacity becomes more significant. Note that when 6hc ~ 0, both the escape 
probability Pesc and the correction factor Cuei are close to unity. This means 
that almost all the emitted photons can escape to infinity and so the ionization 
history returns to Saha equilibrium for He I recombination. 

This simple fitting formula can reproduce quite well the detailed numerical 
resu lt fo r the ionization history at the later stages of He i recombination. From 
Fig. 3.1, we can see that our model with 6ho ~ 0.86 matches with the numerical 
result at z < 2000 ||l^ . Although our fitting model does not agree so well with 
the numerical results for the earlier stages of Hei recombination, the effect on 
the Ci is neligible. This is because the visibility function g{z) = e^'^dr/dz, is 
very low at z > 2000 (at least 16 orders of magnitude smaller than the maximum 
value of g{z)), as shown in the lower panel of Fig. 6.1. Our fitting approach also 
appears to work well for other cosmological models (Switzer & Hirata, private 
communication) . 

In this paper, we employ the fudge factor Fn for H (which is the extra 
factor multiplying an) and the He I parameter 6hc in our model to represent 
the remaining uncertainties in recombination. For Hei, the factors cho and 
6hc in equation (S.ll) are highly correlated. We choose to fix ohc and use 
6hc as the free parameter in this paper; this is because it measures the power 
dependence of the ratio of the relevant rates 7 in the escape probability due to 
the continuum opacity Pcon,H- For hydrogen recombination, all the individual 
updates suggested recently give an overall change less than 0.5 per cent in Xe 
around the peak of the visibility function. Only the effect of considering the 
separate ^-states causes more than a 1 per cent change, and only for the final 
stages of hydrogen recombination (z < 900). Therefore, we think it is sufficient 
to represent this uncertainty with the usual fudge factor i'H, which basically 
controls the speed of the end of hydrogen recombination (see Fig. S.3). The 
best-fit to the current recombination calculation has Fn ~ 1.14. 



6.3 Forecast data 

We use the CosmoMC code to perform a Markov Chain Monte Carlo 
(MCMC) calculation for sampling the posterior distribution with given fore- 
cast data. The simulated Planck data and likelihood function are generated 
based on the settings suggested in Lewis et al. ||l2|. We use full polarization 
information for Planck by considering the temperature T and i?-type polariza- 
tion anisotropics for £ < 2400, and assume that they are statistically isotropic 
and Gaussian. The noise is also isotropic and is based on a simplified model 
with Nf^ ~ = 2 X 10~''/xK^, having a Gaussian beam of 7 arcmin- 

utes (Full Width Half Maximum, |Q). For our fiducial model, we adopt 
the best values of the six cosmological parameters in a ACDM model from 
the WMAP three- year result ijl^ . The six parameters are the baryon density 
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Figure 6.2: Ionization fraction function of redshift z calculated based on 

the modified He i recombination discussed here with different values of the he- 
lium fitting parameter 6hc- The curve with 6hc — (long-dashed, cyan) overlaps 
the line using Saha equilibrium recombination (solid, black). The cosmological 
parameters used in this graph are the same as for Fig. 6.1. 
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Figure 6.3: The ionization fraction function of redshift z calculated with 

different values of the hydrogen fudge fact or F a- The cosmological parameters 
used in this graph are the same as in Fig. 3.1. 
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r^b/i^ = 0.0223, the cold dark matter density fic/i^ = 0.104, the present Hub- 
ble parameter Hq = 73kms~^Mpc~^, the constant scalar adiabatic spectral in- 
dex ris = 0.951, the scalar amplitude (at k = 0.05 Mpc"!) 10^° = 3.02 and the 
optical depth due to reionization (based on a sharp transition) r = 0.09. For 
recombination, we calculate the ionization history using the original Recfast 
with the fudge factor for hydrogen recombination Fh set to 1.14 and the helium 
abundance equal to 0.24. 



6.4 Results 



Fig. |6.4| shows the parameter constraints from our forecast Planck likelihood 
function using the original recfast code with varying Fa and adopting differ- 
ent priors. For the Planck forecast data, can be well constrained away from 
zero (the same result as in |l^) and is bounded by a nearly Gaussian distribu- 
tion with (T approximately equal to 0.1. When we only vary Fn with different 
priors (compared with fixing it to 1.14), it basically does not change the size of 
the error bars on the cosmolo gica l parameters, except for the scalar adiabatic 
amplitude 10^° Ag. From Fig. |6.2| , we can see that the factor Fh controls the 
speed of the final stages of H i recombination, when most of the atoms and elec- 
trons have already recombined. Changing Fh affects the optical depth r from 
Thomson scattering, which determines the overall normalization amplitude of 
Ci (oc e"^"^) at angular scales below that subtended by the size of the horizon 
at last scattering {£ > 100). This is the reason why varying Fh affects the 
uncertainty in A^, since Ag also controls the overall amplitude of Ce (see the 



upper right panel in Fig. 6.6 for the marginalized distribution for Fh and Ag). 
The modified recombination model also changes the peak value (but not really 
the width) of the adiabatic spectral index Ug distribution, as one can see by 



comparing the dotted and dashed curves in Fig. 6.4. 

Based on all the suggested effects on Hi recombination, the uncertainty in Xo 
is at the level of a few per cent at z < 900, which corresponds to roughly a 1 per 



cent change in F^. In Fig. 6A, we have also tried to take this uncertainty into 
account by considering a prior on Fn with a = 0.01 (the long-dashed curves). 
We find that the result is almost the same as for the case using a = 0.1 for the 
-Fh prior. On the other hand, the error bar (measured using the 68 per cent 
confidence level, say) of Ag is increased by 40 and 16 per cent with a — 0.1 and 
0.01, respectively. 



Fig. 6.5 shows the comparison of the constraints in the original and modi- 
fied versions of recfast, with both Hi and He I parameterized. By comparing 
the solid and dotted curves in Fig. we can see that only the peaks of the 
spectra of the cosmological parameters are changed, but not the width of the 
distributions, when switching between the original and modified recfast codes. 
Allowing 6hc to float in the modified recombination model only leads to an in- 
crease in the error bar for spectral index Ug among all the parameters, including 
i^H- For the dashed curves, we used a very conservative prior for &ho, namely a 
flat spectrum from to 1.5 (i.e. from Saha recombination to the old Recfast 
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Figure 6.4: Marginalized posterior distributions for forecast Planck data varying 
the hydrogen recombination only. All the curves are generated using the original 
Recfast code. The solid (black) curve uses fixed F^, while the dotted (red) and 
dashed (green) allow for varying Fh with Gaussian distributions centred at 1.14, 
with a = 0.1 and 0.01, respectively. Note that using a flat prior (between and 
1.5) for Fu gives the same spectra as the case with a = 0.1 (the red dotted line). 
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Figure 6.5: Marginalized posterior distributions for forecast Planck data with 
hydrogen and helium phonomenological parameters both allowed to vary. The 
solid (black) curve shows the constraints using the original recfast code and 
allowing Fh to be a free parameter. The other curves also allow for the variation 
of Fh and use the fitting function for Hei recombination described in Section 
2: the dotted (green) line sets 6hc equal to 0.86; the dashed (red) one is with a 
flat prior for 6hc from to 1.5; and the long-dashed (blue) one is with a narrow 
prior for fcfjci consisting of a Gaussian centred at 0.86 and with a = 0.1. 
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behaviour). We can see that the value of 6hc is poorly constramed, because the 
CMB is only weakly sensitive to the details of He i recombination. Nevertheless, 
this variation allows for faster He i recombination than in the original Recfast 
code and this skews the distribution of rtg towards higher values (see also the 
upper left panel in Fig. 6.6). This is because a faster Hei recombination leads 
to fewer free electrons before H i recombination and this increases the diffusion 
length of the photons and baryons. This in turn decreases the damping scale 
of the acoustic oscillations at high £, which therefore gives a higher value of rig. 
In addition, this variation in 6hc increases the uncertainty (at the 68 per cent 
confidence level) of Ug by 11 per cent. 

Based on the comprehensive study of Switzer & Hirata (2008) [|l^, the domi- 
nant remaining uncertainty in Hei recombination is the 2'^p-l^s transition rate, 
which causes about a 0.1 per cent variation in Xc at 2; ~ 1900. For our fitting 
procedure this corresponds to about a 1 per cent change in foue- We try to 
take this uncertainty into account in our calculation by adopting a prior on bue 



which is peaked at 0.86 with width (sigma) liberally set to 0.1. From Fig. 6.5 



one can see that the error bar on Ug is then reduced to almost the same size as 
found when fixing &hc equal to 0.86 (the dotted and long-dashed curves). This 
means that, for the sensitivity expected from Planck, it is sufficient if we can 
determine 6ho to better than 10 per cent accuracy. 

As well as the individual marginalized uncertainties, we can also look at 
whether there are degeneracies among the parameters. From Fig. 3.6, we see 



that -Fh and 6hc are quite independent. This is because the two parameters 
govern recombination at very different times. As discussed before, ^hc controls 
the speed of Hei recombination, which affects the high-z tail of the visiblity 
function, while Fn controls the low-z part. 



6.5 Discussion and conclusions 

In this paper, we modify recfast by introducing one more parameter feuo (be- 
sides the hydrogen fudge factor ^h) to mimic the recent numerical results for 
the speed-up of He i recombination. By using the COSMOMC code with fore- 
cast Planck data, we examine the variation of these two factors to account for 
the remaining dominant uncertainties in the cosmological recombination calcu- 
lation. For Hei, the main uncertainty comes from the 2^p-l^s rate [|l^, which 
corresponds to about a 1 per cent change in 6ho- We find that this level of varia- 
tion has a negligible effect on the determination of the cosmological parameters. 
Therefore, based on this simple model, if the existing studies have properly con- 
sidered all the relevant physical radiative processes in order to provide Xc to 0.1 
per cent accuracy during He i recombination, then we already have numerical 
calculations which are accurate enough for Planck. 

For H, since there is still no comprehensive model which considers all the 
interactions between the atomic transitions and the radiation background, we 
consider the size of the updates as an indication of the existing level of uncer- 
tainty. We represent this uncertainty by varying the fudge factor Fn, because 
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Figure 6.6: Projected 2D likelihood for the four parameters rig, ^s, -Fh and ^He- 
Shading corresponds to the marginalized probabilities with contours at 68 per 
cent and 95 per cent confidence. 
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the largest update on Xc occurs at z < 900, and comes from a consideration 
of the separate angular momentum states We find that Fn needs to be 
determined to better than 1 per cent accuracy in order to have negligible effect 
on the determination of cosmological parameters with Planck. 

Hydrogen recombination is of course important for the formation of the 
CMB anisotropics C^, since it determines the detailed profile around the peak 
of the visibility function g{z). A comprehensive numerical calculation of the 
recombination of H i (similar to He i) to include at least all the recent suggestions 
for updates on the evolution of Xc is an urgent task. We need to determine 
that the phenomenological parameters Fh and 6hc arc fully understood at the 
< 1 per cent level before we can be confident that the uncertainties in the 
details of recombination will have no significant effect on the determination of 
cosmological parameters from Planck. 
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Chapter 7 

Summary and Future work 



7.1 Effects of distortion photons 

In this thesis, we have presented the detailed profile of the spectral distortion 
to the Cosmic Microwave Background (CMB) due to the Hi Lya and 2s-ls 
two-photon transitions, and the corresponding lines of Hei and Heii. The main 
peak of the distortion is from the Lya line and is located at A = 170 /xm in the 
standard cosmological ACDM model. Although the detection of these spectral 
distortions will be quite challenging due to the presence of the Cosmic Infrared 
Background (CIB), they would provide a direct probe for the detailed physical 
processes during the recombination epoch. These high energy distortion photons 
also have significant effects on the recombination of lithium pTf and formation of 
the primordial molecules in the cosmological 'dark ages' at redshift z < 500. 
Recently, Switzer & Hirata (2005) |^ showed that the distortion photons from 
H I strongly suppress and delay the formation of neutral lithium (Li i) . They 
found that neutral lithium is three orders of magnitude smaller than found in 
previous studies, which assumed a perfect blackbody radiation background (see 
|l6| for reviews) . This dramatically reduces the optical depth of Li i and 
makes the effects of Li i scattering on the CMB anisotropics unobservable |^ . 

Despite the effect of these spectral distortions reducing the strength of some 
potentially observable anisotropy effects, there may be other, related effects 
which are detectable. Basu et al. (2004) [|l| and Hernandez-Monteagudo & Sun- 
yaev (2005) jo) have shown that other sources of line scattering might lead to 
interesting signatures from the z ~3-25 universe. In a seperate study |l7| it 
was suggested that the spectral lines themselves, each with a different effective 
visibility function, could lead to anisotropy signatures which probe different 
epoches. Although all of these effects are relatively weak, as the sensitivity of 
experiments increases, it seems likely that these subtle effects, which are es- 
sentially mixed anisotropy and spectral signatures, will become of increasing 
importance. 

The primordial molecules (for example, II2, HD and LiH) arc important in 
the formation of the first stars, since molecular cooling plays a significant role 
in the first collapse of baryonic matter, when the amplitudes of structures grow 
non- linear and virialize ||^, 11 . With the addition of the distortion photons. 



the abundance of primordial H2 was found to be about 75% less compared 
with previous studies 0| . Note that the cooling of gas is more effective through 
the HD dipole radiation than through the quadrupole radiation from H2, and 
therefore understanding the formation of HD may be very important. Since the 
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main route for the formation of HD is H2 + D+^HD + H+, it will be worth 
performing a follow-up calculation for HD with the updated populations of H2 . 

7.2 A single numerical code for recombination 

From the above discussion, it is clear that the detailed spectrum of the distor- 
tion photons can have strong influence on the formation of primordial molecules. 
The distortion spectrum in turn depends strongly on details of the radiative 
processes in cosmological recombination. But the main motivation of improv- 
ing the recombination calculation is to obtain an accurate visibility function 
for CMB anisotropics. In anticipation of upcoming CMB experiments which 
push to smaller angular scales with higher sensitivity (for example, Planck , 
ACT @ and SPT (l^), it is crucial to understand all the relevant physical pro- 
cesses during recombination which may contribute more than (say) 0.1% to the 
ionization fraction a;o, in order not to bias the cosmological parameter extrac- 
tion. In this thesis, we studied the effect on recombination of the Hei 2^Pi-l^So 
spin-forbidden transition and also the higher order non-resonant two-photon 
transitions (nS-lS and nD-lS) of Hi and He I in a multi-level atom model. We 
found that more than 40% of electrons cascade down to the ground state through 
the 2^Pi-l^So spin-forbidden transition from the n = 2 state, and the inclusion 
of this transition brings more than a 1% change in Xe compared with previous 
studies. We also adopted improved two-photon rates for the transitions from 3S 
to IS and 3D to IS by including all the non-resonant poles through the high- 
lying intermediate nP states (n > 4) [|[ Q . Our best estimated H i non-resonance 
two-photon rates are lower than the ones from Dubrovich & Grachev (2005) |^ 
due to destructive interference in the matrix element; and so from this effect we 
found no more than a 0.5% change in Xo- 

Although in Chapter 4, we only considered the effect of some of these specific 
additional transitions, there have been many other recent updates on recombi- 
nation calculation, as discussed in Chapter 2 and the discussion sections in 
Chapter 3, 4 and 6. Most of the suggested improvements are concerned with 
consistently treating the radiative interactions between matter and the sur- 
rounding photons. We revisited one of the previous studies [p^ , which claimed 
that the matter was reheated by the distortion photons from recombination 
and that this delayed the H l recombination. We found that the energy transfer 
between the distortion photons and the matter (through Compton scattering) 
is very inefhcient, and the resulting effect on Xe is no more than 10~^. This 
is much lower than the previous estimate and hence this effect can be safely 
ignored. 

Many suggestions for improvements to recombination have been carried out 
in different independent numerical codes, and therefore the overall effect of all 
the modifications is still uncertain. Recently, there has been a comprehensive 
study of helium recombination ^ |2|, ^ , which includes most of the physical 
processes relevant of calculating Xc at the 0.1% level. Since 92% of the atoms 
in the Universe are hydrogen, it follows that Hi recombination is considerably 
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more important in determing the detailed profile of the last scattering surface 
for CMB photons. So a remaining task is to perform a similar systematic study 
for H I recombination, or even a full calculation combining the H and He cases. 

Once all the relevant corrections for the detailed numerical recombination 
calculation have been solidified, we need to incorporate a modified approximate 
version of these effects into a fast code similar to recfast for incorporating 
into the Boltzmann codes (for example, CMBFAST and CAMB [^) which are 
used for calculating the CMB anisotropics, C^s. This is because the current 
detailed numerical recombination calculations take far too long (typically more 
than a day) to yield results for a single cosmological model. In the previous 
chapter, we introduced an extra parameter 6ho in the current recfast to ap- 
proximately model the speed-up of Hei recombination due to the continuum 
opacity of Hi. This modified recfast can be considered as the first step in 
parametrizing the other recent result from the detailed numerical codes into a 
simple three-level atom calculation. 

We also studied how varying 6hc along with the existing hydrogen fudge 
factor .Fh might account for some of the remaining uncertainties in recombina- 
tion. Using the COSMOMC code with Planck forecast data {£ < 2500), we found 
that we need to determine the effective value of 6hc to better than 10% and 
Fh to better than 1%. The current Hei recombination studies seem to already 
calculate Xc accurately enough for Planck, but we still require a comprehen- 
sive study for Hi to reach the same level of accuracy. Note that these two 
phenomenological parameters mainly affect the determination of the scalar am- 
plitude As and the spectral index rig of the primordial perturbation spectrum. 
There are other CMB experiments, such as the Atacama Cosmology Telescope 
(ACT) which will be able to measure CiS over a wide range of angular scales 
(1000 < ^ <10000); such measurements can put tight constraints on the tilt of 
the temperature power spectrum, which is characterized by the primodial spec- 
tral index Ug. For these and even better future experiments, we may need to 
determine these two phenomemological parameters (-Fh and 6hc) to better than 
the 1% level in order to obtain the correct inferences about inflationary models. 
Alternatively, we should systematically account for all the relevant updates on 
recombination, in additional to the one recent correction which we included in 
the modified recfast code. There is still much work to be done! 
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